Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,874	C→T	noncoding (67/1542 nt)	rrsC →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,874	0	C	T	100.0%	65.7 / NA	29	noncoding (67/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/17);  total (12/17)

ACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACC  >  NC_000913/3941631‑3942115                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  acCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGGTCTTTAACAATTTATCAGACAATCTGTGAGGGCTCTCGAAGATACGGAGTCTTAACGGCGCAAGTCGAAAAAGGAATACCAAGTCGCAAGAGTGAAGACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATGGAAGAGTTTGATCATGGCGCATATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACa                                                                                                                                                                                                                                            <  1:521967/251‑1 (MQ=11)       cgccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAAc                                                                                                                                                                                                                                      <  2:275030/251‑1 (MQ=17)          cgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGc                                                                                                                                                                                                                                   >  1:411298/1‑251 (MQ=18)                 aaaaaGCAGAGCGGCACTGCTCTGTAACTATTTATCAGACAATCTGTGTGGTCACTCGAAGATACGGAATCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAATAGTTAACACGTAATTCATTACGAAGGTTAATTCGTGGAGCATCAAACGTTTAAATTGAATACATTGATCATGGCTCAGATTGAACGCTTTCGGCTTGCCTAACACATGTAAGTCGAATGGTAATACGAAACTGCTTGCTGt                                                                                                                                                                                                                            <  2:283803/251‑1 (MQ=11)                      gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGc                                                                                                                                                                                                                       <  1:68761/251‑1 (MQ=18)                       caaagcGGCACTGCTCTTTAACAATTTATCAGACAAGCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAGCTTTGAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGTCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGGTTCGCt                                                                                                                                                                                                                      <  2:268247/251‑1 (MQ=11)                       aaaagcGGCTCTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACTGATTCTTAACGGCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGTAGTTTAATTCTTTGAGCATCAAACTTTTAAATTTAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAATTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCt                                                                                                                                                                                                                      <  1:354859/250‑1 (MQ=11)                          agcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAAAGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGCAACAGCTTGCTGTTTCGCTGAc                                                                                                                                                                                                                   >  1:42158/1‑251 (MQ=17)                                                                                                                                         tAATTCATATCGAAGGTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATTGCTCAGATTGAACTCTGGAGGCAGGCCGTACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGTACTGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACGACTGGAAACGGTAGCTAATTCCGCATAACGTCGCAAGACCAAAGAGGGGGCCCtt                                                                                                     <  2:312332/250‑1 (MQ=11)                                                                                                                                              cATTACGAAGTTTAATTCTTTGAGCATCAAACGTTTAAATGGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGGCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCAGAACGTCGCAAGGCCAAAGAGGGGGACCTTCGGGc                                                                                                <  2:507767/250‑1 (MQ=255)                                                                                                                                                                      cATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCTCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGAt                                                                       <  1:363796/251‑1 (MQ=18)                                                                                                                                                                           aaCTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGTGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGTGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGAGGGGAt                                                                  <  1:135712/251‑1 (MQ=16)                                                                                                                                                                               tttAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGACACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTATTGCCATCAGATGTGCCCAGATGGGATTAGc                                                              >  1:194175/1‑251 (MQ=18)                                                                                                                                                                                 tAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCt                                                             >  1:202286/1‑250 (MQ=18)                                                                                                                                                                                 gaaaTTGAAGAGTTTGATCATGGCTCAGATTGAACGCGGGCGGCAGTCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCGTGCTGTTTCGCTGACGATTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGTGGGACCTTCGGGCCTCTTGCCATCAGATGTGTCCAGATGGGATTAGCta                                                            <  1:480491/250‑1 (MQ=11)                                                                                                                                                                                         agagTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTCCTGGAAACGGTAGCTAATACCGCCTAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCCGATGGGCCCAGATGGGGTTAGCTCGTAggggg                                                    >  2:75277/1‑248 (MQ=11)                                                                                                                                                                                               tGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTaa                                               <  1:245405/250‑1 (MQ=37)                                                                                                                                                                                                 aTCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATCACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGGTTAGCTAGTAGGGGGGGTAACgg                                            >  2:182925/1‑251 (MQ=255)                                                                                                                                                                                                 aTCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGTGGAACgg                                            <  2:518764/251‑1 (MQ=255)                                                                                                                                                                                                  tCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCCTCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGc                                           >  2:421766/1‑251 (MQ=37)                                                                                                                                                                                                   cATGGATCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGCTGGGATTAGCTCGTAGGTGGGGTAACGGcc                                          >  2:577486/1‑250 (MQ=255)                                                                                                                                                                                                      ggCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGTGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCTGATGGGCCCAGGTGGGATTAGCTAGTAGGTGGGGTTACGGCTCAc                                       <  2:299779/251‑1 (MQ=255)                                                                                                                                                                                                      ggCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGGCGAATGGGAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGGCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAc                                       <  2:36035/251‑1 (MQ=255)                                                                                                                                                                                                         tCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAAGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAAATACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGGAGGAGGGGTAACGGCTCACCtc                                    >  1:247534/1‑250 (MQ=255)                                                                                                                                                                                                             aTTGAACGATGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACAACAGAGGGGGACCTTCGGGCCTCTTGCAATCAGCTGTGCCCAGATGGGATTAGATAGTAGGTGGGGGAACGGCTCACCttgg                                  >  1:401300/1‑246 (MQ=11)                                                                                                                                                                                                                gAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCGGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgac                              <  1:74998/250‑1 (MQ=255)                                                                                                                                                                                                                                      acaatcAAGTTGAATGGTAACAGGAATCTGCTTGCTTTTTCTCTGACGAGTGGCGGACGGGTTAGTAATGTCTGGGAAACTTCCTGATGGAGGGGGATAACTACTGGAATCGGTATCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCCCCTAGGCGACGATCCCTCTCTGGTCtcatagga       <  2:466063/246‑8 (MQ=11)                                                                                                                                                                                                                                           cAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGAcc  >  1:28055/1‑251 (MQ=38)                                                                                                                                                                                                                                             aGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGAcc  >  1:504146/1‑249 (MQ=39)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  ACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACC  >  NC_000913/3941631‑3942115

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