Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,731,312	G→A	intergenic (‑155/+288)	rrsG ← / ← clpB	16S ribosomal RNA of rrnG operon/protein disaggregation chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,731,312	0	G	A	100.0%	159.5 / NA	47	intergenic (‑155/+288)	rrsG/clpB	16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/28);  total (19/28)

GATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC  >  NC_000913/2731275‑2731406                                      |                                                                                               gATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                  >  2:346034/1‑100 (MQ=14) gATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                  >  1:171159/1‑100 (MQ=14)   ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                <  1:657946/100‑1 (MQ=17)   ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                <  2:1932785/100‑1 (MQ=17)    tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                               <  1:653396/100‑1 (MQ=18)    tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                               <  1:1472399/100‑1 (MQ=18)     gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                              <  1:966678/100‑1 (MQ=25)       cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  2:782831/1‑100 (MQ=255)        tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                           <  2:2459397/100‑1 (MQ=255)        tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                           >  1:2634398/1‑100 (MQ=255)          aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                         <  1:2347835/100‑1 (MQ=255)          aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                         <  1:798229/100‑1 (MQ=255)          aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                         <  1:270616/100‑1 (MQ=255)           tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                        <  1:838377/100‑1 (MQ=255)             aaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                      >  1:1562287/1‑100 (MQ=255)               ttGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt                    >  1:2139272/1‑100 (MQ=255)                tGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt                   <  2:1128261/100‑1 (MQ=255)                 gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  1:297551/100‑1 (MQ=255)                 gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  1:540552/100‑1 (MQ=255)                 gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:989824/1‑100 (MQ=255)                 gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:197679/100‑1 (MQ=255)                  ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                 <  2:665189/100‑1 (MQ=255)                  ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                 >  1:1070642/1‑100 (MQ=255)                  ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                 >  1:1574581/1‑100 (MQ=255)                   tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGTCGCGGGGGGTGCATAATACGCCTTCCCGCTACAGAGGCAAGCATTTCTTTGCTTTTCTCTGTTgag                >  2:124504/1‑100 (MQ=17)                   tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                <  2:1400542/100‑1 (MQ=255)                   tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                <  1:356123/100‑1 (MQ=255)                   tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                >  1:2100440/1‑100 (MQ=255)                   tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGGCAAGCATTTCTTTGCTTTTCTCTGTTgag                >  2:2629152/1‑100 (MQ=255)                    aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgagn               <  1:965736/100‑2 (MQ=255)                    aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga               <  1:965376/100‑1 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              >  1:927756/1‑100 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              <  1:2767340/100‑1 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              >  2:1683718/1‑100 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              >  1:1373051/1‑100 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              >  1:1418962/1‑100 (MQ=255)                     aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTCTCTCTGTTGAGAt              <  1:1856522/100‑1 (MQ=255)                       gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc            >  1:2501043/1‑100 (MQ=255)                        agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct           <  1:222302/100‑1 (MQ=255)                         gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc          <  1:1154618/100‑1 (MQ=255)                         gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc          <  2:213744/100‑1 (MQ=255)                          cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa         <  2:489260/100‑1 (MQ=255)                           cgatctgcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg        <  1:415210/94‑1 (MQ=255)                           aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg        <  1:103589/100‑1 (MQ=255)                           aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg        >  1:446231/1‑100 (MQ=255)                            gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTTTTTGCTTTTCTCTGTTGAGATTCTCAg        >  1:1798081/1‑99 (MQ=255)                            gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       <  1:358934/100‑1 (MQ=255)                             tGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga      <  1:2028905/100‑1 (MQ=255)                                 gcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgagatcg  >  2:415206/1‑94 (MQ=255)                                      |                                                                                               GATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC  >  NC_000913/2731275‑2731406

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