Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,131,098	A→G	A27A (GCA→GCG)	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,098	0	A	G	100.0%	131.2 / NA	36	A27A (GCA→GCG)	flgB	flagellar basal‑body rod protein FlgB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/20);  total (16/20)

CGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACA  >  NC_000913/1130968‑1131232                                                                                                                                   |                                                                                                                                       cgTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATc                                                                                                                               >  2:339408/1‑140 (MQ=255)    aaaGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGcc                                                                                                                            <  2:233464/140‑1 (MQ=255)          aTTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGcc                                                                                                                      >  1:397503/1‑140 (MQ=255)            ttAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGa                                                                                                                    <  1:443759/140‑1 (MQ=255)               aGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATAc                                                                                                                 <  2:397503/140‑1 (MQ=255)                gCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATAcc                                                                                                                <  1:817155/140‑1 (MQ=255)                  tGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATAcccc                                                                                                              <  2:40182/140‑1 (MQ=255)                  tGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATAcccc                                                                                                              <  2:692555/140‑1 (MQ=255)                    aTGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTg                                                                                                            >  1:688859/1‑140 (MQ=255)                            aTTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAg                                                                                                    <  2:688859/140‑1 (MQ=255)                                                   atGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAg                                                                             >  2:85661/1‑140 (MQ=255)                                                        cGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAAc                                                                        >  1:790479/1‑140 (MQ=255)                                                          aTAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCA‑GCGCGCGATATCGATTTTGCCAGTGAACTTa                                                                     >  1:525276/1‑140 (MQ=255)                                                           tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTa                                                                     >  2:425528/1‑140 (MQ=255)                                                           tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTa                                                                     >  1:622267/1‑140 (MQ=255)                                                           tAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCA‑GCGCGCGATATCGATTTTGCCAGTGAACTTaa                                                                    <  2:525276/140‑1 (MQ=255)                                                                        gccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCa                                                        <  1:164317/140‑1 (MQ=255)                                                                         ccTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCaa                                                       >  1:202600/1‑140 (MQ=255)                                                                                ttttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGAc                                                <  1:618497/140‑1 (MQ=255)                                                                                ttttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGAc                                                <  1:124370/140‑1 (MQ=255)                                                                                 tttCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACg                                               >  1:71020/1‑140 (MQ=255)                                                                                   tCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATAGCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACggg                                             >  2:601432/1‑140 (MQ=255)                                                                                    caacaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGa                                            <  1:426409/140‑1 (MQ=255)                                                                                    caacaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGa                                            <  1:425528/140‑1 (MQ=255)                                                                                          gagGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAAc                                      <  2:790479/140‑1 (MQ=255)                                                                                           agGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAAcc                                     <  1:339408/140‑1 (MQ=255)                                                                                            ggCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCa                                    >  1:135410/1‑140 (MQ=255)                                                                                             gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCa                                        <  1:93427/135‑1 (MQ=255)                                                                                             gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCa                                        >  2:93427/1‑135 (MQ=255)                                                                                             gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAg                                   >  1:455710/1‑140 (MQ=255)                                                                                                   aaTCTGCGCGCCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGt                             <  2:71020/140‑1 (MQ=255)                                                                                                        gcgcgcCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCAc                        <  2:622267/140‑1 (MQ=255)                                                                                                          gcgcCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTg                      <  1:42428/140‑1 (MQ=255)                                                                                                           cgcCCAGCGTCAGGAAGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGa                     >  2:500065/1‑140 (MQ=255)                                                                                                                         aaGTGCTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACg       <  1:539235/140‑1 (MQ=255)                                                                                                                              cTGGCGGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca  >  2:696348/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                       CGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACA  >  NC_000913/1130968‑1131232

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: ‑