Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,131,195	A→G	T60A (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,195	0	A	G	100.0%	98.6 / NA	29	T60A (ACC→GCC)	flgB	flagellar basal‑body rod protein FlgB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/14);  total (15/14)

AGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCG  >  NC_000913/1131058‑1131328                                                                                                                                          |                                                                                                                                      agGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGcc                                                                                                                                     >  2:423630/1‑140 (MQ=255)   gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAg                                                                                                                                   <  2:261029/140‑1 (MQ=255)   gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAg                                                                                                                                   >  1:827960/1‑140 (MQ=255)     gcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAgtg                                                                                                                                 >  2:265484/1‑140 (MQ=255)                 cgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGa                                                                                                                     >  1:760337/1‑140 (MQ=255)                                    cTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                                                                                  >  1:729479/1‑140 (MQ=255)                                        cagcaAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCt                                                                                                            <  1:202324/126‑1 (MQ=255)                                        cagcaAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCt                                                                                                            >  2:202324/1‑126 (MQ=255)                                                     cAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTAATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                                                                                 >  2:364027/1‑140 (MQ=255)                                                       aTGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                               >  2:381445/1‑140 (MQ=255)                                                               aCCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcc                                                                       <  2:760337/140‑1 (MQ=255)                                                                   cTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGa                                                                              <  2:390653/129‑1 (MQ=255)                                                                   cTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGa                                                                              >  1:390653/1‑129 (MQ=255)                                                                     ggTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                               <  2:44968/126‑1 (MQ=255)                                                                     ggTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                               >  1:44968/1‑126 (MQ=255)                                                                       ttATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAg                                                               <  1:381445/140‑1 (MQ=255)                                                                           cAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAg                                                               <  1:374086/136‑1 (MQ=255)                                                                           cAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAg                                                               >  2:374086/1‑136 (MQ=255)                                                                               cgcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAATCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCaa                                                       <  1:282720/140‑1 (MQ=255)                                                                                   cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATAcc                                                   >  1:769913/1‑140 (MQ=255)                                                                                             tttGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                         <  1:539051/140‑1 (MQ=255)                                                                                                gCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCa                                      >  1:495361/1‑140 (MQ=255)                                                                                                  cAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                    >  2:83475/1‑140 (MQ=255)                                                                                                                 aGTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTa                     >  1:691262/1‑140 (MQ=255)                                                                                                                  gTCATGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCATGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTaa                    <  2:729479/140‑1 (MQ=255)                                                                                                                     aTGCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                 <  2:769913/140‑1 (MQ=255)                                                                                                                       gCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCg               <  1:83475/140‑1 (MQ=255)                                                                                                                       gCAACGTGGACGGGATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCg               <  2:495361/140‑1 (MQ=255)                                                                                                                                    gATGCAGCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATcg  <  2:642747/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      AGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCG  >  NC_000913/1131058‑1131328

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