Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,131,261	A→G	T82A (ACC→GCC)	flgB →	flagellar basal‑body rod protein FlgB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,131,261	0	A	G	100.0%	110.4 / NA	32	T82A (ACC→GCC)	flgB	flagellar basal‑body rod protein FlgB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/15);  total (17/15)

TTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAA  >  NC_000913/1131128‑1131395                                                                                                                                      |                                                                                                                                       ttATCAGGCGCGCGATATCGATGTTGCCAGTGAACTTTAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAg                                                                                                                                  <  1:291514/140‑1 (MQ=255)    tCAGGCGCGCGATATCGATTTTCCCAGTGAAATTACAAAGGTCATGCAGCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAAc                                                                                                                               <  2:216179/140‑1 (MQ=255)          gcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTACGGCGGAGGCGCTGACCCCTCCTGCCGCAGAACTGCAAt                                                                                                                         >  2:330733/1‑140 (MQ=255)                  tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGTTGACGCCTCCTGCCGCAGAACTGCAATACCGTAtt                                                                                                                 <  2:691403/140‑1 (MQ=255)                  tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGTTGACGCCTCCTGCCGCAGAACTGCAATACCGTAtt                                                                                                                 <  2:324721/140‑1 (MQ=255)                  tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTAtt                                                                                                                 >  1:31419/1‑140 (MQ=255)                          gCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCa                                                                                                                           <  1:525491/122‑1 (MQ=255)                          gCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCa                                                                                                                           >  2:525491/1‑122 (MQ=255)                                    ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCtt                                                                                               <  1:658954/140‑1 (MQ=255)                                               aTGCAACATGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                                                                    <  2:677465/140‑1 (MQ=255)                                                tGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAcc                                                                                   >  2:173523/1‑140 (MQ=255)                                                 gCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCg                                                                                  <  2:93130/140‑1 (MQ=255)                                                 gCAACGTGGACGGGATGAAACCAGTGTGGTTACACTGACAATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCg                                                                                  <  2:234241/140‑1 (MQ=255)                                                      gTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGat                                                                             >  2:638170/1‑140 (MQ=255)                                                      gTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTACAATACCGTATTACGGACCAGCCTTCGCTTGACGGTAAAACCGTCGat                                                                             >  2:257084/1‑140 (MQ=255)                                                                 gCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGa                                                                  >  1:609062/1‑140 (MQ=255)                                                                    aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGAAGGTAATACCGTCGATATGGATCGCGAACg                                                               >  1:243082/1‑140 (MQ=255)                                                                        gtgtGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGAGCAGCCTTCACTTGACGGTAATACCGGCGATTTGGATCGCGAACGCAcc                                                           >  2:203945/1‑140 (MQ=255)                                                                        gtgtGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc                                                           >  1:564101/1‑140 (MQ=255)                                                                         tgtgGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa                                                             >  2:134877/1‑137 (MQ=255)                                                                         tgtgGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa                                                             <  1:134877/137‑1 (MQ=255)                                                                              ttGCACTGACGATGACCCCAACTCAACACATTCCGCCGCCGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt                                                     <  2:187533/140‑1 (MQ=255)                                                                                gCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGc                                                   <  1:173523/140‑1 (MQ=255)                                                                                        gATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACa                                           >  1:590491/1‑140 (MQ=255)                                                                                            aCCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGTGAACGCACCCAGTTTGCCGATAACAGCCt                                       >  1:269393/1‑140 (MQ=255)                                                                                                cAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                   >  2:31658/1‑140 (MQ=255)                                                                                                cAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                   >  1:668895/1‑140 (MQ=255)                                                                                                  aCGCAACACATTCCGGCCCAGGAGCTGACGCTTCCTGCCTCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATa                                 <  2:269393/140‑1 (MQ=255)                                                                                                      aaCACATTCCGGCGCAGGCGCTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAg                             >  1:184355/1‑140 (MQ=255)                                                                                                                         gcTGACGCCTCCTGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGa          <  2:609062/140‑1 (MQ=255)                                                                                                                                cctcctGCCGCAGAATTGCAATTCCGTATTCCGGACCTGCCTTCGCTTGACGGTAAAACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCa   <  1:257084/140‑1 (MQ=255)                                                                                                                                cctcctGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCa   <  1:638170/140‑1 (MQ=255)                                                                                                                                 ctcctGCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCaa  >  2:148620/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       TTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAA  >  NC_000913/1131128‑1131395

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