Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	223,514	A→G	intergenic (+106/‑257)	gmhB → / → rrsH	D‑glycero‑beta‑D‑manno‑heptose‑1,7‑bisphosphate 7‑phosphatase/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	223,514	0	A	G	100.0%	139.7 / NA	39	intergenic (+106/‑257)	gmhB/rrsH	D‑glycero‑beta‑D‑manno‑heptose‑1,7‑bisphosphate 7‑phosphatase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/18);  total (21/18)

AGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAG  >  NC_000913/223386‑223647                                                                                                                                 |                                                                                                                                      agcagcAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCcgg                                                                                                                            <  1:392384/140‑1 (MQ=255) agcagcAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCcgg                                                                                                                            <  1:436416/140‑1 (MQ=255)   cagcaAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCcggtt                                                                                                                          <  1:287856/140‑1 (MQ=255)     gcaAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCcggttcg                                                                                                                        >  1:537092/1‑140 (MQ=255)                              gATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc                                                                                               <  2:634133/140‑1 (MQ=255)                                  aGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAggg                                                                                           >  2:629122/1‑140 (MQ=255)                                  aGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAggg                                                                                           >  1:94302/1‑140 (MQ=255)                                      gTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt                                                                                                 <  1:137679/130‑1 (MQ=255)                                      gTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt                                                                                                 >  2:137679/1‑130 (MQ=255)                                       ttGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGAc                                                                                      <  2:424909/140‑1 (MQ=255)                                       ttGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGAc                                                                                      <  2:557478/140‑1 (MQ=255)                                         gAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGActc                                                                                    >  1:375229/1‑140 (MQ=255)                                         gAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGActc                                                                                    <  1:629122/140‑1 (MQ=255)                                                      aTTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGaaa                                                                               >  2:707657/1‑132 (MQ=255)                                                      aTTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGaaa                                                                               <  1:707657/132‑1 (MQ=255)                                                       tttttCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGc                                                                      <  1:571491/140‑1 (MQ=255)                                                        ttttCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCg                                                                     <  1:666310/140‑1 (MQ=255)                                                        ttttCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCg                                                                     <  2:375229/140‑1 (MQ=255)                                                               cTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAAtata                                                              >  1:689786/1‑140 (MQ=255)                                                                    cTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCa                                                         <  2:171597/140‑1 (MQ=255)                                                                         tGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTcg                                                    <  2:537092/140‑1 (MQ=25)                                                                              cGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAAGCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACa                                               <  1:664779/140‑1 (MQ=14)                                                                               gACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAg                                              >  2:480924/1‑140 (MQ=11)                                                                                  tCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGa                                           >  2:69861/1‑140 (MQ=11)                                                                                           tGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGc                                  >  2:77184/1‑140 (MQ=37)                                                                                            gcgcCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGcc                                 >  2:389261/1‑140 (MQ=38)                                                                                              gcCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgc                               >  1:576845/1‑140 (MQ=255)                                                                                              gcCTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgc                               >  1:571514/1‑140 (MQ=255)                                                                                               ccgccATCGAAACGGCGGATGGGAATCACTTCACGCAAAAAGCCGGCGCGGTTGAATAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCAACTCGCGACAGTGAGCTGAAAGCcgcg                              >  2:684373/4‑140 (MQ=14)                                                                                               ccTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgcg                              >  1:457546/1‑140 (MQ=255)                                                                                               ccTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgcg                              >  2:580243/1‑140 (MQ=255)                                                                                                cTCCATCGACACGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGt                             >  1:345559/1‑140 (MQ=255)                                                                                                        acacGGCGGATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                     <  2:689786/140‑1 (MQ=255)                                                                                                            ggcggATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctct                 >  2:614544/1‑140 (MQ=255)                                                                                                               ggATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                     <  2:196257/133‑1 (MQ=255)                                                                                                               ggATGTGAATCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                     >  1:196257/1‑133 (MQ=255)                                                                                                                      aaTCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAAttt       >  1:33608/1‑140 (MQ=255)                                                                                                                        tCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                       >  1:323517/1‑122 (MQ=255)                                                                                                                        tCACTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                       <  2:323517/122‑1 (MQ=255)                                                                                                                           cTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                   <  2:170571/123‑1 (MQ=255)                                                                                                                           cTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                   >  1:170571/1‑123 (MQ=255)                                                                                                                           cTTCACGCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAg  >  2:682291/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                      AGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAG  >  NC_000913/223386‑223647

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