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,8610AG93.5% 96.3 / ‑6.1 31intergenic (+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 (11/18);  minor base T (2/0);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTA  >  NC_000913/1196621‑1197074
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
tagttaTGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccat                                                                                                                                                                                                              <  1:88176‑M1/250‑5 (MQ=255)
    taTGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACCCCCGCTCAAATATTGGTAAATTGCCGGTGTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCTCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaac                                                                                                                                                                                                          <  2:180146‑M1/250‑9 (MQ=255)
          ttCACTGATATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCACTTCTTGTTAAATTGCCGGATTTGCATCAACTATACACACGGGACTCTCCAGGGGCCAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACCTGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATACTTAGTACTGAGACTGTTTAACCTAGGGAATTATGTTCGAAccatcaacgattaa                                                                                                                                                                                                    >  2:453186‑M1/1‑236 (MQ=255)
           tCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaac                                                                                                                                                                                                  >  2:357430‑M1/1‑235 (MQ=255)
             aCTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactg                                                                                                                                                                                                <  1:357430‑M1/251‑19 (MQ=255)
                    ttaGGACAATTCTGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGACCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcg                                                                                                                                                                                         >  2:118516‑M1/1‑226 (MQ=255)
                           aaTTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagt                                                                                                                                                                                  >  2:397920‑M1/1‑219 (MQ=255)
                            aTTATGTATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagtt                                                                                                                                                                                 <  2:340965‑M1/251‑34 (MQ=255)
                            aTTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGAATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagtt                                                                                                                                                                                 >  1:347450‑M1/1‑218 (MQ=255)
                                     tACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttac                                                                                                                                                                        >  1:186867‑M1/1‑209 (MQ=255)
                                              ttACACACGCTAAAATATTGTTAAATTGCCGGATTTGTATCAAATACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGAACATGCGGAATGTAAAAAAAGCGGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattacactgcgcttcggcgcagttttcgtttacaggatgttg                                                                                                                                                               >  2:148954‑M1/1‑200 (MQ=255)
                                                   cccGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGCCTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagg                                                                                                                                                          >  2:479711‑M1/1‑195 (MQ=255)
                                                         cAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaa                                                                                                                                                     <  1:510273‑M1/250‑62 (MQ=255)
                                                                    aaaTTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaa                                                                                                                                          >  1:459946‑M1/1‑178 (MQ=255)
                                                                         gCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaa                                                                                                                                     <  1:624410‑M1/250‑78 (MQ=255)
                                                                                ttGTATCAACTACCCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgcc                                                                                                                              <  2:252104‑M1/250‑85 (MQ=255)
                                                                                    aTCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagtt                                                                                                                          <  1:460084‑M1/250‑89 (MQ=255)
                                                                                             ccccccgggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACGGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgtggaaggggaaattcctggggcaaaaaaaccccgccagttaaggcggtaa                                                                                                                <  2:527844‑M1/248‑99 (MQ=255)
                                                                                                  cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggggaactc                                                                                                           >  2:322785‑M1/1‑148 (MQ=255)
                                                                                                          gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcggggaaacccatcctatt                                                                                                   >  1:412300‑M1/1‑140 (MQ=255)
                                                                                                                                  gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataat                                                                                  <  2:397796‑M1/244‑129 (MQ=255)
                                                                                                                                  gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaaactcatcctattggggaaaaaagaataaaaaatcgg                                                                           >  1:397796‑M1/1‑116 (MQ=255)
                                                                                                                                       cAATCACCTTAGCGTTCAGGTCCATGGGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTAAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcgggg                                                                      <  2:459946‑M1/251‑141 (MQ=255)
                                                                                                                                            ccctTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatc                                                                 <  1:284569‑M1/250‑146 (MQ=255)
                                                                                                                                                     ttCAGGTGCATGCGGAATGTAAAAAAGGCCGCTATCGCGGACCCTACACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatggggaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctccca                                                         <  2:329700‑M1/250‑154 (MQ=255)
                                                                                                                                                         ggTACAGGCGGAATGTAAAAAAGGCCGGGAGCGCGGCCCCTTCACAGACATCTTTAGTACTGGGACTGTTTAACCTAGGGAAGTATTGTCGTAccataaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtata                                                    <  2:624563‑M1/251‑159 (MQ=255)
                                                                                                                                                                       gTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaat                                      <  1:396218‑M1/251‑173 (MQ=255)
                                                                                                                                                                            aaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTTTCGTAccatcaacgcttaaactgcgcttcggcgccgttttcgttt                                                                                                                                                                          >  2:225248‑M1/1‑74 (MQ=255)
                                                                                                                                                                            aaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgttt                                                                                                                                                                          <  1:225248‑M1/114‑41 (MQ=255)
                                                                                                                                                                                      ggcgcgGCCCCTTCACTTACCTCTGTAGTACGTAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggttttatcatctcccagtatatccatactaacaataaggttatttactc                        <  2:379396‑M1/249‑187 (MQ=255)
                                                                                                                                                                                                              ttAGTACTTAGACTGTTTAACCTAGGGATTTATTGGCGTAccacaaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaaggaaaaattttggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacattttgtt  <  2:347450‑M1/248‑209 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
TAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTA  >  NC_000913/1196621‑1197074

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