Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,684,366	C→A	R216S (CGT→AGT)	ECB_03438 →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,684,366	0	C	A	100.0%	270.7 / NA	78	R216S (CGT→AGT)	ECB_03438	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (44/34);  total (44/34)

ACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGG  >  REL606/3684318‑3684413                                                 |                                                aCTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgt                                               <  1:2575655/51‑1 (MQ=255) aCTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgt                                               >  1:4193812/1‑51 (MQ=255) aCTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgt                                               <  1:4631988/51‑1 (MQ=255)   tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtgg                                             >  1:3126304/1‑51 (MQ=255)   tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtgg                                             <  1:2152178/51‑1 (MQ=255)   tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtgg                                             >  1:2509138/1‑51 (MQ=255)    gAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggt                                            >  1:2326478/1‑51 (MQ=255)    gAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggt                                            >  1:6734160/1‑51 (MQ=255)    gAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggt                                            >  1:5769185/1‑51 (MQ=255)    gAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggt                                            >  1:3318208/1‑51 (MQ=255)     aaGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggtg                                           >  1:4967872/1‑51 (MQ=255)       gagCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgcg                                         <  1:4047579/51‑1 (MQ=255)        agCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgcgc                                        <  1:5052788/51‑1 (MQ=255)         gCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCa                                       >  1:1406398/1‑51 (MQ=255)          ccGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCaa                                      >  1:4415060/1‑51 (MQ=255)          ccGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCaa                                      >  1:5195985/1‑51 (MQ=255)           cgcgCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgc                                          <  1:4524136/46‑1 (MQ=255)           cgcgCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAAc                                     >  1:4321448/1‑51 (MQ=255)            gcgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAAcc                                    <  1:4657676/51‑1 (MQ=255)            gcgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAAcc                                    >  1:5988254/1‑51 (MQ=255)             cgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCg                                   <  1:891189/51‑1 (MQ=255)             cgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCg                                   >  1:4662497/1‑51 (MQ=255)              gcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGc                                  >  1:4997231/1‑51 (MQ=255)              gcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGc                                  >  1:3182075/1‑51 (MQ=255)              gcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGc                                  >  1:5490/1‑51 (MQ=255)               cTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCa                                 <  1:4769457/51‑1 (MQ=255)                tACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCAc                                >  1:5654245/1‑51 (MQ=255)                tACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCAc                                >  1:1962488/1‑51 (MQ=255)                tACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCAc                                >  1:4823765/1‑51 (MQ=255)                tACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCAc                                >  1:665723/1‑51 (MQ=255)                 aCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACt                               >  1:3869315/1‑51 (MQ=255)                 aCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACt                               >  1:347797/1‑51 (MQ=255)                   cTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAt                             <  1:5017084/51‑1 (MQ=255)                   cTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAt                             <  1:6235125/51‑1 (MQ=255)                    tGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAtt                            <  1:4690899/51‑1 (MQ=255)                    tGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAtt                            <  1:3765743/51‑1 (MQ=255)                     ggCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTa                           >  1:3823840/1‑51 (MQ=255)                     ggCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTa                           <  1:4352963/51‑1 (MQ=255)                     ggCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTa                           >  1:6229280/1‑51 (MQ=255)                     ggCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTa                           <  1:521381/51‑1 (MQ=255)                      gcgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTAc                          <  1:6841752/51‑1 (MQ=255)                       cgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCTCACTATTACg                         <  1:5024383/51‑1 (MQ=255)                       cgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg                         <  1:6911968/51‑1 (MQ=255)                       cgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg                         <  1:1936854/51‑1 (MQ=255)                       cgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg                         <  1:1868821/51‑1 (MQ=255)                         cTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAc                       >  1:665472/1‑51 (MQ=255)                          tGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAcc                      >  1:6308743/1‑51 (MQ=255)                           gACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCa                     >  1:2212751/1‑51 (MQ=255)                           gACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCa                     >  1:6208241/1‑51 (MQ=255)                           gACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCa                     <  1:435686/51‑1 (MQ=255)                           gACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCa                     >  1:6950172/1‑51 (MQ=255)                            aCGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCaa                    >  1:4525682/1‑51 (MQ=255)                            aCGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCaa                    <  1:211692/51‑1 (MQ=255)                             cGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAg                   <  1:6934808/51‑1 (MQ=255)                             cGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAg                   >  1:1929677/1‑51 (MQ=255)                             cGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAg                   <  1:2171688/51‑1 (MQ=255)                              gAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGa                  >  1:826854/1‑51 (MQ=255)                               aaCTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGaa                 >  1:1755240/1‑51 (MQ=255)                                aCTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAat                <  1:5247957/51‑1 (MQ=255)                                 cTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAata               >  1:1047879/1‑51 (MQ=255)                                 cTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAata               <  1:2349517/51‑1 (MQ=255)                                 cTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAata               >  1:26901/1‑51 (MQ=255)                                  tATTTTGTTGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAatat              <  1:4578621/51‑1 (MQ=255)                                  tATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAatat              <  1:4403494/51‑1 (MQ=255)                                  tATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAatat              <  1:2069183/51‑1 (MQ=255)                                  tATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAatat              >  1:525515/1‑51 (MQ=255)                                    ttttGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGa            >  1:3800943/1‑51 (MQ=255)                                      ttGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaaa          <  1:4077315/51‑1 (MQ=255)                                       tGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaaaa         <  1:5301734/51‑1 (MQ=255)                                         tCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAgc       <  1:1103212/51‑1 (MQ=255)                                         tCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAgc       >  1:7078060/1‑51 (MQ=255)                                          cGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAgcg      >  1:1939244/1‑51 (MQ=255)                                           gAACAGAGTGGTGTGCAACCGCACTATTACGACCAAGAATATGAAAAgcgc     >  1:1805237/1‑51 (MQ=255)                                            aaCAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGcgcg    <  1:5589974/51‑1 (MQ=255)                                            aaCAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGcgcg    >  1:6278343/1‑51 (MQ=255)                                            aaCAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGcgcg    <  1:2170652/51‑1 (MQ=255)                                             aCAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCgg   >  1:6288098/1‑51 (MQ=255)                                              cAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCggg  >  1:698012/1‑51 (MQ=255)                                                 |                                                ACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGG  >  REL606/3684318‑3684413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑