Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,313,176 (T)8→7 noncoding (93/93 nt) micF → small regulatory RNA MicF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,313,1690T.75.0% 30.2 / 2.9 21noncoding (86/93 nt)micFsmall regulatory RNA MicF
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (6/9);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

AAATGCAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAATAAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGCAATTATTTTTGTTGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGG  >  NC_000913/2312926‑2313386
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                               
caaaggggaaacttttttaAAGATTATTGCGGAATGGCGAAATAAAGCACCTAACATCAAGCAATAATAATTCAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACta                                                                                                                                                                      <  1:1263962/286‑1 (MQ=255)
     cAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAATAAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAACTATAACGCTTGTTAAc                                                                                                                                                        <  1:1263961/306‑1 (MQ=255)
                                     cGAAATAAAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaa                                                                                                                                                                            <  1:1263966/257‑1 (MQ=255)
                                                          aaGCAATAATAATTCAAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCa                                                                                                                  <  1:1263968/293‑1 (MQ=255)
                                                              caatataatTCAAGGTTTAAAAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCTCATTACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTCCCCTATTTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACtt                                                                                                                                                                                                      >  1:415792/5‑209 (MQ=255)
                                                                   taatTCAAGGTT‑‑‑AAAATCAATTAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                                                                                                                                           <  1:1263969/244‑1 (MQ=255)
                                                                                  taaaTCAATAACTTATTCTAAGTATTTGACAGCACTGAATGTCAAAACAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGACTGACTTCGATT‑‑GTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGtt                                                                                                         <  1:1263970/278‑1 (MQ=255)
                                                                                      tCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTATCCCTTCTTTACACACTTTTCATTATTCt                                                                                                                                                                                           >  1:415795/1‑197 (MQ=255)
                                                                                                     ttAAGTATTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGA‑‑TGCCTCGCATTCGGttt                                                                                                                                                                                                                               >  1:1263972/1‑148 (MQ=255)
                                                                                                                          aTGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATaa                                                                                                                           >  1:1263974/1‑225 (MQ=255)
                                                                                                                               aaaaacaaaaaCCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCtgtg                                                                                                                                                                                        >  1:415798/2‑162 (MQ=255)
                                                                                                                                        aacCTTCACTCGCAACTAGAATACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGAGCTGACTTCGGATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCt                                                                                                                                                                                      <  1:415796/156‑1 (MQ=255)
                                                                                                                                                   gCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGG‑‑TTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaaactat                                                                                                                                                                       <  1:1263973/158‑5 (MQ=255)
                                                                                                                                                             tAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTATCCCTTCTTTACACACTTTTCa                                                                                                                                                                                                  >  1:415800/1‑123 (MQ=255)
                                                                                                                                                             tAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGCAATTATTTTTGTTGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACgg  >  1:1263976/1‑310 (MQ=255)
                                                                                                                                                                                     ttatttatTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTAATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGCAATTATTTTTGTTGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGgatga              >  1:1263978/1‑275 (MQ=255)
                                                                                                                                                                                           attaCCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCa                                                                                                                     <  1:1263977/166‑1 (MQ=255)
                                                                                                                                                                                                      tcattTCTGAATGTCTGTTTACCCCTATTTCAACCGGGA‑‑TGCCTCGCATTCGGTTTGTTTTCACCCTTCTTTACACACTTTTCATATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCttgt                                                                                                           <  1:415802/169‑3 (MQ=255)
                                                                                                                                                                                                                      tGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTATCCCTTCTTTACACACTTTTCa                                                                                                                                                                                                  >  1:415803/1‑67 (MQ=255)
                                                                                                                                                                                                                         ttACCCCTATTTCAACCGGA‑‑TGCCTCGCATTTCGGTTTTTTTACCCTTCTTTACACACttt                                                                                                                                                                                                     >  1:415804/1‑61 (MQ=255)
                                                                                                                                                                                                                                     tCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTACCCTTCtttt                                                                                                                                                                                                              >  1:415805/1‑40 (MQ=255)
                                                                                                                                                                                                                                      aaaCCGGA‑‑TGCCTCGCATTCGGTTTTTTTTTACCCTTCTTTACACACTTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGTCAATTATTTTTGTTGGGCTTTTTGTTAGGTTAtttt                                                              >  1:1263980/2‑184 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                               
AAATGCAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAATAAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAATAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAACAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCTGAATGTCTGTTTACCCCTATTTCAACCGGA‑‑TGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGCAATTATTTTTGTTGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGG  >  NC_000913/2312926‑2313386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: