Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,986,329	C→T	intergenic (‑322/‑103)	aslA ← / → glmZ	putative sulfatase AslA/small regulatory RNA GlmZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,986,329	0	C	T	100.0%	56.4 / NA	21	intergenic (‑322/‑103)	aslA/glmZ	putative sulfatase AslA/small regulatory RNA GlmZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/13);  total (8/13)

AAATCAGGAGAGCGTTTTCAATCCTACCTCTGGCGCAGTTGATATGTAAGGCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTGCCTCATAAACTCCGGAATGACGCAGAGCCGTTTACGGTGCTTATCGTCCACTGACAGATGTCGCTTATGCCTCATCAGACACCATGGACACAACGTTG  >  NC_000913/3986147‑3986566                                                                                                                                                                                       |                                                                                                                                                                                                                                              aaaTCAGGAGAGCGTTTGCAATCCTACCTCTGGCGCATTTGATATGTTAGGCAGGTTTATTATATCGTGTTGAGTATTGATGCTGTTTTTAGTTTTAACGGCAATTCATATATGTGTTAGTAATTTAATGAATTTTATCATGCATAATAAATATGTGTAAGATCAAGCTCAGGTTAAAGATTTTCTCAGGAAGATATGACACAGGAAGCAAAGAGGATTACAGAATTATCGCATAACAAGTGTTAAGGGAt                                                                                                                                                                           <  1:505135/251‑1 (MQ=255)          gagCGTTTTCAATCCTACCTCTGGCGCAGTTGATATGTAAGGCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATTTGTTATTAATTGAATGAATTTTATCATTCATAATAATTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAg                                                                                                                                                                                                                                         <  1:448346/180‑1 (MQ=255)                                ggCGCAGTTGATATGTAAGGCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTAc                                                                                                                                                                                                                              <  1:398260/169‑1 (MQ=255)                                             tGTAAGGCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCtg                                                                                                                                                        <  1:260108/226‑1 (MQ=255)                                                   gCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTCCTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCTGCTTCTCTGTGGCATAATACACGAGTAGATGCTCATTcc                                                                                                                          >  1:456583/1‑250 (MQ=255)                                                    cAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTgg                                                                                                                                                     >  1:237736/1‑222 (MQ=255)                                                                   gcgTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTACTATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCAtctc                                                                                                                     >  1:461877/1‑239 (MQ=255)                                                                                 tGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGNgg                                                                                                                                                                                                              >  2:312134/1‑136 (MQ=255)                                                                                     gTTTTTAGGTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGGATGTGTAGGATCAAGCTCAGGTTAAATATTTACTAAGGAAGTTATTACTCAGGAAGCAAAAAGGATTATAGAATTATCTCATAACAATAGTTAAGGTATGTTATTTACCGATTCTCTGTTGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCTTTAGTGCTTAATCAACTcc                                                                                       >  1:108796/1‑251 (MQ=255)                                                                                       ttttAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTGATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTAGTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGcg                                                                                                                                <  1:145048/208‑2 (MQ=255)                                                                                              tttAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTGCGCCTTAGTGCCTCATAAACTCCGGAATGan                                                                               <  1:312134/250‑2 (MQ=255)                                                                                                   cGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTGCCTCATAAACTCCGGAATGACGCaga                                                                          <  1:178326/250‑1 (MQ=255)                                                                                                      cAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGAttn                                                                                                                                                            <  1:101023/165‑2 (MQ=255)                                                                                                                             tgaatgaaTTTTATCATTCATAATAAGTATGGTTAGGATCAAGCTCAGGTTAAATATTTACTCAGGAAGTTCTTAATCAGGaa                                                                                                                                                                                                                       >  1:292838/1‑83 (MQ=255)                                                                                                                              gaatgaatTTTATCATTCATAATAAGTATGTGTAGGATCAAGCGCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCGCGTATGTTCGCTTTAGTGCCTCATAAAc                                                                                          <  1:184210/207‑1 (MQ=255)                                                                                                                                                           gtgtAGGATCAAGCTCAGGTTAAAGATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATg                                                                                                                                                                          <  1:218011/98‑1 (MQ=255)                                                                                                                                                                        cTCAGGTTAAATATTTACTCAGGTAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTCTCTCATAACAAGTGTTAAGGGATGTTATTTCCAGCTTCTCTGTGGCATAATACACGAGTAGATGCTCATTCCATATCTTATGTTCGCTTTAGTGCCTCATAAACTCCGGAATGACGCAGAGCCGTTTACGGTGCTTATCGTCCACTGACAGATGTCGCTTATg                               >  1:503892/1‑224 (MQ=255)                                                                                                                                                                        cTCAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTCTCTCATAACAAGTGTTCAGGGATGTTATTTCCAGATTCTCTGTGGCCTAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTg                                                                                                    >  1:336812/1‑155 (MQ=255)                                                                                                                                                                          cAGGTTAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTGCCTCATAAACTCCGGAATGACGCAGAGCCGTTTACGGTGCTTATCGTCCACTGACAGATGTCGCTTATGCCTCATCAGACACCATGGACACAACGTTg  <  1:397613/251‑1 (MQ=255)                                                                                                                                                                               tAAATATTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTGCCTCATAAACTCCGGa                                                                                    <  1:232982/164‑1 (MQ=255)                                                                                                                                                                                   tatTTACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTATCTGTgg                                                                                                                                                     <  1:145217/95‑1 (MQ=255)                                                                                                                                                                                       |                                                                                                                                                                                                                                              AAATCAGGAGAGCGTTTTCAATCCTACCTCTGGCGCAGTTGATATGTAAGGCAGGTTTATTATATCGCGTTGATTATTGATGCTGTTTTTAGTTTTAACGGCAATTAATATATGTGTTATTAATTGAATGAATTTTATCATTCATAATAAGTATGTGTAGGATCAAGCTCAGGTTAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTGTTAAGGGATGTTATTTCCCGATTCTCTGTGGCATAATAAACGAGTAGATGCTCATTCCATCTCTTATGTTCGCCTTAGTGCCTCATAAACTCCGGAATGACGCAGAGCCGTTTACGGTGCTTATCGTCCACTGACAGATGTCGCTTATGCCTCATCAGACACCATGGACACAACGTTG  >  NC_000913/3986147‑3986566

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