Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,861 A→G intergenic (+79/+6) C0293 → / ← ymfD small regulatory RNA C0293/putative SAM‑dependent methyltransferase YmfD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,8610AG96.2% 76.7 / NA 26intergenic (+79/+6)C0293/ymfDsmall regulatory RNA C0293/putative SAM‑dependent methyltransferase YmfD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (8/17);  minor base T (1/0);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CCAAAGTTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAG  >  NC_000913/1196614‑1197091
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
ccAAAGGTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCTCCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGTGGCCCCTTCACATACATCTTTAGTACTGAGGCTGTTTAACCTAGGGAATTATTGTCg                                                                                                                                                                                                                                     <  2:82474/251‑1 (MQ=255)
   aaGTTAGTTATTTATTCACTGTTATTAGGACAATTATGAATTCCCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAc                                                                                                                                                                                                                                  <  1:501212‑M1/251‑2 (MQ=255)
     gttagttaTGTATGCACTGTTATTAGGACAATTATGAATGACCACAACTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAATGCCGCGAGCGCGGCCCCTTCACATACATCTATAGTACTGAGACTGATTAACCTAGGGAATTATTGTCGtaac                                                                                                                                                                                                                                 >  2:782/1‑248 (MQ=255)
        agttaTGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatc                                                                                                                                                                                                                              <  2:322896‑M1/250‑6 (MQ=255)
                                 cAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTACCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcag                                                                                                                                                                                                    <  1:157825‑M1/251‑32 (MQ=255)
                                  aaTTATGAATTACCACTCCTTACACCTGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTCTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagt                                                                                                                                                                                                   >  1:163010‑M1/1‑219 (MQ=255)
                                                     ttACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGCATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttg                                                                                                                                                                                >  2:517813‑M1/1‑200 (MQ=255)
                                                         acCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaa                                                                                                                                                                              >  1:465925‑M1/1‑196 (MQ=255)
                                                            cGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaaggga                                                                                                                                                                         <  1:289250‑M1/251‑59 (MQ=255)
                                                                 caatcttgttaaATTGCCAGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGGATGAAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTTTCGTAccatcaacgattaaactgcaattagacgaagtttacgtttacaggatgttgaaagggaaaatt                                                                                                                                                                    >  2:487692‑M1/6‑188 (MQ=255)
                                                                 aaaTATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACCTGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcg                                                                                                                                                                                                       >  2:355401‑M1/1‑188 (MQ=255)
                                                                 aaaTATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcg                                                                                                                                                                                                       <  1:355401‑M1/216‑29 (MQ=255)
                                                                                 ccGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACCGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgacagggaaaattctggggcaaaaacagc                                                                                                                                                    >  2:50751‑M1/1‑172 (MQ=255)
                                                                                  cGGTTTTGTATCAACAACTCACCCGTGACTCGCCAGGGGACAGCCAACAGGCATTGTGTGCAATCACTTTAGCGTTCAGGTACATGCGTAATGTAAAAAATGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGTGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaatgctggtgcaaaaaaagcc                                                                                                                                                   <  2:5241‑M1/251‑81 (MQ=255)
                                                                                                ctactCACCCCGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGAGCGGCCCCTTCACATACATCTTGAGTACTGAGACTGTTTAACCTAGGGTATTATTGTCGTAccatcaatgattaacatgcgcttcggcgcagttttcgtttacaggatgttgcaagggaaacttctggggcaaaaaaagctcgccagttaacgcg                                                                                                                                     >  2:47056‑M1/1‑157 (MQ=255)
                                                                                                     cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaa                                                                                                                                <  1:371860‑M1/251‑100 (MQ=255)
                                                                                                      ccccgGGACTGCCCCGGGGACAGCCAACAGGCATTGGGTGAAACCCCCTTGGCTTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCTGCCCCTTCACATACATCTTTAGTTCTGAGACAGTTTAACCGAGGGAATTATTGTCGTAccatcaacgattaaactgcgctgcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaac                                                                                                                               <  2:358499‑M1/250‑101 (MQ=255)
                                                                                                       cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCCCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacc                                                                                                                              <  1:536714‑M1/251‑102 (MQ=255)
                                                                                                         cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAGTGTAAAAAAGGCCGCGAGCGCGGCCCTTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctc                                                                                                                            <  1:503607‑M1/251‑104 (MQ=255)
                                                                                                           ggattcgccAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcat                                                                                                                          <  1:499708‑M1/247‑106 (MQ=255)
                                                                                                              cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTCGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCACGGCCCCTTCACATACATCTTTAGTACTGAGCATGTTTAACCTCGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgcaagggaacattctggggcaaaaaaagccagccagttacggcgggacaactccaaac                                                                                                                       >  2:348127‑M1/1‑143 (MQ=255)
                                                                                                                    aGGGGACAGCCATCAGTCATTGGGTGCAATCCCCTTAGCGTTCAGGTACAGGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgctgcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatggga                                                                                                                 <  1:98019‑M1/251‑115 (MQ=255)
                                                                                                                          cAGCCAACAGGCATTGGGTGCAATCCCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCTAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaa                                                                                                           <  1:90134‑M1/251‑121 (MQ=255)
                                                                                                                                                                                        ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttattt                                             <  1:87568‑M1/251‑183 (MQ=255)
                                                                                                                                                                                                               acatCTTTAGTTCTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacatttt                      <  1:336156‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                                                   gTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacattttgttttgtgcgtgggaacagc  <  1:465316‑M1/251‑226 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
CCAAAGTTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAG  >  NC_000913/1196614‑1197091

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