Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	961,913	A→G	intergenic (+29/‑82)	cmk → / → rpsA	cytidylate kinase/30S ribosomal subunit protein S1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,913	0	A	G	100.0%	132.7 / NA	39	intergenic (+29/‑82)	cmk/rpsA	cytidylate kinase/30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/20);  total (19/20)

CGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA  >  NC_000913/961784‑962046                                                                                                                                  |                                                                                                                                      cGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGTATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGAc                                                                                                                             <  2:811938/140‑1 (MQ=255) cGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGAc                                                                                                                             <  2:29713/140‑1 (MQ=255)   gCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGGATTTGCAGTACCCCCGTTGCGATGGAATGACAg                                                                                                                           <  2:137771/140‑1 (MQ=255)      gCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCgg                                                                                                                        <  2:41149/140‑1 (MQ=255)       ccGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGAc                                                                                                                             >  2:490406/1‑134 (MQ=255)       ccGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGAc                                                                                                                             <  1:490406/134‑1 (MQ=255)                     ttGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAAcc                                                                                                         <  1:689051/140‑1 (MQ=255)                             caccacCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCgg                                                                                                 >  2:809909/1‑140 (MQ=255)                                    ttAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAg                                                                                          <  1:156026/140‑1 (MQ=255)                                             gAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGa                                                                                 >  1:786515/1‑140 (MQ=255)                                                 aaGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGtt                                                                             <  1:192629/140‑1 (MQ=255)                                                   gTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTaa                                                                           >  2:581521/1‑140 (MQ=255)                                                   gTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTaa                                                                           >  1:395840/1‑140 (MQ=255)                                                             aaGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCt                                                                 <  1:674137/140‑1 (MQ=255)                                                             aaGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCt                                                                 >  2:674137/1‑140 (MQ=255)                                                                  cTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGa                                                            >  1:380687/1‑140 (MQ=255)                                                                    aCAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGAtt                                                          >  2:825283/1‑140 (MQ=255)                                                                    aCAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGAtt                                                          <  2:332815/140‑1 (MQ=255)                                                                       aTACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTaaa                                                       <  2:480795/140‑1 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTg                                                                <  1:80258/123‑1 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTg                                                                >  2:80258/1‑123 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAAc                                                      <  1:67062/133‑1 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAAc                                                      >  2:67062/1‑133 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTg                                               >  2:305759/1‑140 (MQ=255)                                                                               gcCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTg                                               >  2:776104/1‑140 (MQ=255)                                                                                ccAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGa                                              <  2:796862/140‑1 (MQ=255)                                                                                ccAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGa                                              >  1:277716/1‑140 (MQ=255)                                                                                     aaTTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCtt                                         <  1:46611/140‑1 (MQ=255)                                                                                            tctcGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCa                                  <  2:786515/140‑1 (MQ=255)                                                                                              tcGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAAc                                >  1:649884/1‑140 (MQ=255)                                                                                               cGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAAct                               >  2:233535/1‑140 (MQ=255)                                                                                                      gCGACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGaa                        <  2:417091/140‑1 (MQ=255)                                                                                                        gACCGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAaga                      <  2:10297/140‑1 (MQ=255)                                                                                                          ccGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGt                    >  1:396855/1‑140 (MQ=255)                                                                                                          ccGAATTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGt                    >  2:829592/1‑140 (MQ=255)                                                                                                             aaTTTGCAGTACCCCCGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCt                 <  2:649884/140‑1 (MQ=255)                                                                                                                        cccccGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCTTAAAAGAAATc      >  1:647261/1‑140 (MQ=255)                                                                                                                        cccccGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCTTAAAAGAAATc      >  2:808136/1‑140 (MQ=255)                                                                                                                            cGTTGCGATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGaaa  <  2:277716/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      CGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACTCTTTGAAGAGTCCTTAAAAGAAATCGAAA  >  NC_000913/961784‑962046

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