Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,731,312	G→A	100%	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%	146.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 (23/24);  total (23/24)

TGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCC  >  NC_000913/2731278‑2731408                                   |                                                                                                 tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                 <  2:524911/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                 <  1:222478/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                 <  1:116903/100‑1 (MQ=18)  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                <  2:437334/100‑1 (MQ=25)  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                >  2:1419692/1‑100 (MQ=25)  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                <  2:1133392/100‑1 (MQ=25)  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                <  1:876053/100‑1 (MQ=25)    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  1:1208429/1‑100 (MQ=255)    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              <  1:2040946/100‑1 (MQ=255)    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  1:466477/1‑100 (MQ=255)    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGAGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  2:1320180/1‑100 (MQ=255)     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTCCAGAGGCAAGCATTTCTTTGCttt                             >  2:529623/1‑100 (MQ=255)     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGAGGGGTGTGCATAATACGCCTTCCCGCTACAGAGCAAAGGATTTCTTTGCttt                             >  2:529580/1‑100 (MQ=14)      gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTTTTTTGGTGTGCATAATACGCCTTCCTGCTACTGAGTCAAGCATTTCTTTGCtttt                            >  2:1033323/1‑100 (MQ=9)      gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                            >  2:359868/1‑100 (MQ=255)       ataaacttgtAAAGAGCAATCCCGTTTCACTCTTTCGCAGCGACCCGAGGTGTATATAATACACCTTCCCGCTACAGAGCCAAGCATTTCTCTCCTTTtc                           <  2:1728497/91‑1 (MQ=2)       aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                           <  1:565583/100‑1 (MQ=255)       aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                           <  2:2330537/100‑1 (MQ=255)        tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                          >  1:594190/1‑100 (MQ=255)         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         >  2:582397/1‑100 (MQ=255)         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         >  2:2158948/1‑100 (MQ=255)          aaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         <  2:1498667/99‑1 (MQ=255)          aaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                        >  1:34244/1‑100 (MQ=255)          aaTGGTTAAAGAGCAGTGCCGATTCAATGTTTCGAAGCGGCACGTTGTGTGCATAATACGCCTTCGAGCAACAAGAACACGCATTCCTTTGAGTTTctct                        <  2:2024253/100‑1 (MQ=2)            gatctAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt                      <  2:1788867/96‑1 (MQ=255)              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCTGCGCGGGGTGTGGATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtn                     <  2:811479/99‑2 (MQ=255)              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    >  1:75077/1‑100 (MQ=255)              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    <  2:1709230/100‑1 (MQ=255)              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    >  1:884796/1‑100 (MQ=255)              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    <  1:412207/100‑1 (MQ=255)               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                   <  1:286216/100‑1 (MQ=255)               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                   >  1:674357/1‑100 (MQ=255)               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                   >  1:468361/1‑100 (MQ=255)               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGAGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                   >  2:912479/1‑100 (MQ=255)                tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATATGCCTTCCCGATACAGGGTCAAGCAGTTCTTTGCTTTTCTCTGTTgag                  >  2:1611535/1‑100 (MQ=17)                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtagatc                 >  1:1788882/1‑95 (MQ=255)                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgagn                 <  1:748556/100‑2 (MQ=255)                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                 <  1:761388/100‑1 (MQ=255)                   agagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAtt               >  1:405090/1‑100 (MQ=255)                    gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc              <  2:2158890/100‑1 (MQ=255)                    gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc              >  2:270412/1‑100 (MQ=255)                    gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc              <  1:2228266/100‑1 (MQ=255)                        aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg          <  2:430397/100‑1 (MQ=255)                        aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg          <  1:1775194/100‑1 (MQ=255)                        aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg          <  1:279765/100‑1 (MQ=255)                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg         >  1:842047/1‑100 (MQ=255)                           gCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag       >  1:1846036/1‑100 (MQ=255)                             cggttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg         <  2:2196021/93‑1 (MQ=255)                                ttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAccc  <  2:2234214/100‑1 (MQ=255)                                   |                                                                                                 TGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCC  >  NC_000913/2731278‑2731408

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