Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	4,029,087	A→T	noncoding (197/1542 nt)	rrsA →	16S ribosomal RNA of rrnA operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,029,087	0	A	T	100.0%	63.7 / NA	34	noncoding (197/1542 nt)	rrsA	16S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (12/22);  total (12/22)

AGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGAAGAAGGCCTTCGGG  >  CP009273/4028863‑4029315                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     aGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCgg                                                                                                                                                                                                                       >  1:104417/1‑240 (MQ=14)                   tCAAACTTTTAAATTGAAGAGTTTTATCATGGCTCAGATTGAACGCTGGCTGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGAt                                                                                                                                                                                                      <  1:435861/239‑1 (MQ=14)                       aCTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGAt                                                                                                                                                                                                      <  2:319727/235‑1 (MQ=14)                           ttAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                                                                                                                                                                             <  1:345210/240‑1 (MQ=14)                                                      aGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                                                                                                                                                                             <  2:77439/213‑1 (MQ=14)                                                         ttGAACGCTGGCGGCAGGCCTAACACATGCAAATCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCt                                                                                                                                                               <  1:308955/240‑1 (MQ=14)                                                                  ggcggcAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAg                                                                                                                                                                                    >  2:417673/1‑210 (MQ=14)                                                                  ggcggcAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAg                                                                                                                                                        <  1:146749/238‑1 (MQ=14)                                                                            cTAACACAGGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTc                                                                                                                                                              <  2:316579/222‑1 (MQ=14)                                                                                    tGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAg                                                                                                                                                        >  1:452135/1‑220 (MQ=14)                                                                                                    aCAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                                                                                                    <  1:324409/240‑1 (MQ=14)                                                                                                        gacggagCTTGCTTCTTTGCTGACGAGTGGCGGACGTGTGAGTAATGTCTGGGAAACTGCCTGATTGAGGGGGATAATTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCATCGGGCCTCTTGCCAt                                                                                                                                                                                                           <  2:136518/144‑1 (MQ=11)                                                                                                        gaagaaGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGGGGGGTAACGGCTCACCTAGGc                                                                                                                                                      >  1:263762/1‑202 (MQ=14)                                                                                                                       tttGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                                  <  1:288834/239‑1 (MQ=14)                                                                                                                       tgtGCTGACGAGGGGCGGACGGGTGAGTAAGGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                                                                                                   <  1:134124/236‑1 (MQ=14)                                                                                                                               cGAGTGGCGGACGGGTGAGTAAGTTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTgg                                                                                                            <  2:285604/221‑1 (MQ=14)                                                                                                                                      cGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAGAACTACTGGAACCGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCt                                                                                                                                   <  2:330044/191‑1 (MQ=14)                                                                                                                                        gACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAAAACGGCATAACGTCGCAAGACTAATGAGGGGTACATTCGGGACTCTTGCCATCGGATGTGCCCAGATGGGAt                                                                                                                                                                                       <  1:310697/137‑1 (MQ=11)                                                                                                                                               gAGTAATGGCGGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCa                                                                                           <  1:324385/222‑1 (MQ=14)                                                                                                                                                aGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCGGAGAGGATGACCAGCCACAATGGAACTGAGACACg                                                                                                >  1:309681/1‑216 (MQ=14)                                                                                                                                                      gTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg                                                                  <  1:236500/240‑1 (MQ=14)                                                                                                                                                                                   aaCTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGa                                      <  1:481749/239‑1 (MQ=14)                                                                                                                                                                                    actactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGa                                      <  1:146251/238‑1 (MQ=14)                                                                                                                                                                                             aaCGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGact                                                              >  2:71466/1‑203 (MQ=14)                                                                                                                                                                                                ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGcc                                                                                                                                                                                                 <  2:94315/71‑1 (MQ=255)                                                                                                                                                                                                ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGcc                                                                                                                                                                                                 >  1:94315/1‑71 (MQ=25)                                                                                                                                                                                                ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGc                                                                                                                                         >  2:291671/1‑127 (MQ=18)                                                                                                                                                                                                ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGc                                                                                                                                         <  1:291671/127‑1 (MQ=17)                                                                                                                                                                                                ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTAc                                                                                  >  1:77110/1‑182 (MQ=14)                                                                                                                                                                                                     cTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTa                   <  1:121289/240‑1 (MQ=14)                                                                                                                                                                                                           ccGCATACCGTCCCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGGGCCCAGATGGGATTa                                                                                                                                                                                     <  2:172905/72‑1 (MQ=11)                                                                                                                                                                                                              cATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATgaagca            >  1:367205/1‑236 (MQ=14)                                                                                                                                                                                                              cATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGCAGGTGGGGTAATGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACgg                                                                                >  2:208600/1‑170 (MQ=12)                                                                                                                                                                                                                      cGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCCCCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTCTGAAGAAGGCCTTCggg  >  1:327759/1‑240 (MQ=14)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     AGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGAAGAAGGCCTTCGGG  >  CP009273/4028863‑4029315

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