Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,997,033	G→T	noncoding (113/137 nt)	sdsN →	small regulatory RNA SdsN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,997,033	0	G	T	100.0%	122.7 / NA	36	noncoding (113/137 nt)	sdsN	small regulatory RNA SdsN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/13);  total (23/13)

TATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCA  >  NC_000913/1996915‑1997165                                                                                                                       |                                                                                                                                     tataTTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGAttt                                                                                                                 >  1:401206/1‑140 (MQ=255)         tCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCttat                                                                                                         >  2:944403/1‑140 (MQ=255)         tCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCttat                                                                                                         >  2:744092/1‑140 (MQ=255)         tCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCttat                                                                                                         >  2:662798/1‑140 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACTTGTGAATCGTCGGGGCTGAt                                                                                                                   <  2:275694/126‑1 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACTTGTGAATCGTCGGGGCTGAt                                                                                                                   >  1:275694/1‑126 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGAt                                                                                                                   <  2:904860/126‑1 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGAt                                                                                                                   >  1:904860/1‑126 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTattatt                                                                                                     >  2:832950/1‑140 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTattatt                                                                                                     >  1:429404/1‑140 (MQ=255)             gCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTattatt                                                                                                     >  2:509164/1‑140 (MQ=255)                       aCTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGaaa                                                                                           >  1:596693/1‑140 (MQ=255)                        cTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAAc                                                                                          <  2:401206/140‑1 (MQ=255)                           ttGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTtc                                                                                       >  1:91308/1‑140 (MQ=255)                              cAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCg                                                                                    >  2:935602/1‑140 (MQ=255)                              cAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCg                                                                                    <  2:429404/140‑1 (MQ=255)                                    aTCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAAc                                                                              <  2:91308/140‑1 (MQ=255)                                             tctTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTg                                                                     <  1:509164/140‑1 (MQ=255)                                                tATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCgg                                                                  >  1:226276/1‑139 (MQ=255)                                                  tGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAgg                                                                <  1:662798/140‑1 (MQ=255)                                                       ttaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGc                                                                          >  1:422707/1‑125 (MQ=255)                                                       ttaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGc                                                                          <  2:422707/125‑1 (MQ=255)                                                       ttaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGATCGGGgccgct                                                           >  1:965149/1‑140 (MQ=255)                                                                       aTTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGGGAGGGTGCTCGGGGGCGGCAAATaa                                           >  1:462256/1‑140 (MQ=255)                                                                       aTTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGgctgc                                                            <  2:305749/123‑1 (MQ=255)                                                                       aTTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGgctgc                                                            >  1:305749/1‑123 (MQ=255)                                                                            ttAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGGAAAAAAAGCtt                                      >  2:362770/1‑140 (MQ=255)                                                                            ttAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGGTGCTCGGGGGCGGGAAAAAAAGCtt                                      >  2:899750/1‑140 (MQ=255)                                                                              aaGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCtttt                                    >  2:788175/1‑140 (MQ=255)                                                                                      ttATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGGTGCTCGGGGGCGGCAAATAAAGCTTTTGCGGcccc                            >  2:695229/1‑140 (MQ=255)                                                                                         tCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGGGGCAAATAACGCTTTTTCGGccccccc                         >  2:863931/1‑138 (MQ=255)                                                                                              cTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGa                    >  1:661672/1‑140 (MQ=255)                                                                                                           ccAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTg       <  1:899750/140‑1 (MQ=255)                                                                                                            cAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTgg      <  1:863931/140‑1 (MQ=255)                                                                                                             gagacagACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGt     <  1:401862/134‑1 (MQ=255)                                                                                                                cccGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCa  <  1:362770/140‑1 (MQ=255)                                                                                                                       |                                                                                                                                     TATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCA  >  NC_000913/1996915‑1997165

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