Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,153,326	G→A	noncoding (16/144 nt)	sibA →	small RNA SibA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,153,326	0	G	A	100.0%	23.0 / NA	11	noncoding (16/144 nt)	sibA	small RNA SibA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/5);  total (6/5)

GGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAGGTCCGGTATCATACCGACTGCAATAATGAAATAAGAGCAAACGAAGAAACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCC  >  NC_000913/2153105‑2153539                                                                                                                                                                                                                              |                                                                                                                                                                                                                      ggCAAACGTGCTTTGTTCGCTCATTAGAAATCCTTTCTTAGAGAATTATCCAGGTCCGGTATCATACCGACTGCAATAATGCAATAAGAGCAACCGAAGAAACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGAGACATAATTTCAGGAGGGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTa                                                                                                                                                                                           >  1:38913/1‑250 (MQ=255)                                             ntantCGNGNNCNNGTANCATACCGACTGCAATAATGAANNAAGAGCAAANGAAGAAACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTNGCAAGGACGNNGANCCGTCNANNNANNNTANGGGCGGAGTGACATAAATTNAGNNGTGAGGNTTAGGGAGAGATTtccccctcccccc                                                                                                                                                                                                        <  2:40558/189‑2 (MQ=255)                                                           gTATCATACCGACTGCAATAATGAAATAAGAGCAAACGAAGAAACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAAccn                                                                                                                                <  1:103097/251‑2 (MQ=255)                                                                                        gagCAAACGAAGAAACTTTCAGAATTAGGCCCGGCGGAGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCAGTCGATTGACACTAAGGGCGGAGTGACATAATTTAAGGAGTGAGGGTTAGGGAGAGATTTCCCAATCCCCCTGGTGTTATTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTatgat                                                                                                                                                    >  1:76265/1‑202 (MQ=255)                                                                                                 nagaaACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGGGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGAATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCttt                                                                                          >  1:42208/2‑251 (MQ=255)                                                                                                                    gCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCACGACCGGAGAGACCTCGATTGACAATAAGGGTGGAGTGCCCTAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGAGTTCTTAGTAAGCCTGGAAGCTACTCACTAAGAGTATCACCAGTATGCTGACGTGCTTCATCATAACCCTTTCCTTATTAACAGCCCTCTTCTCCGGTAGACGCTTTCCAGTTTCAGCGTCCCGCt                                                                       >  1:84989/1‑251 (MQ=255)                                                                                                                            ncgcAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGaa                                                                                                                                                                                  >  1:49056/2‑136 (MQ=255)                                                                                                                                                                        cGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTGCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCGGTCGCCCTTAAGGCAGTCTATCGCTAAcgcg                   <  1:145973/251‑1 (MQ=255)                                                                                                                                                                                        nggAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTc   >  1:43855/2‑251 (MQ=255)                                                                                                                                                                                         ggAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTAGTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCTCTAACGCGTGAATAAAATAGTTTcc  <  1:52473/251‑1 (MQ=255)                                                                                                                                                                                                             gTGAGGGTTAGGGAGAGATTTCCCCCTCCCCCTGGTGTTCTTATTAAGCCTGGAAGCGAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCgg                                                            <  1:140509/173‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                      GGCAAACGTGATTTGTTCGCTCATAAGAAATCCTTTCTTAGAGAATTATCGAGGTCCGGTATCATACCGACTGCAATAATGAAATAAGAGCAAACGAAGAAACTTTCAGAATTAGGCCCGGCGGCGCAGAAAATTTTGCTATTCGCCTCGCAAGAACGGAGAGCCGTCGATTGACACTAAGGGCGGAGTGACATAATTTCAGGAGTGAGGGTTAGGGAGAGGTTTCCCCCTCCCCCTGGTGTTCTTAGTAAGCCTGGAAGCTAATCACTAAGAGTATCACCAGTATGATGACGTGCTTCATCATAACCCTTTCCTTATTAAAAGCCCTCTTCTCCGGGAGAGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAATGCAGTCTATCGCTAACGCGTGAATAAAATAGTTTCC  >  NC_000913/2153105‑2153539

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