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,8610AG95.2% 55.5 / ‑5.6 21intergenic (+79/+6)C0293/ymfDsmall regulatory RNA C0293/putative SAM‑dependent methyltransferase YmfD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (7/13);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAGTTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGA  >  NC_000913/1196617‑1197080
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                           
aaGTTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAc                                                                                                                                                                                                                       >  2:340875‑M1/1‑250 (MQ=255)
                          aGGACAATGATTAATTACCCCTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGATTTGTATCAACTCCTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTGCACATACATCCTTAGTACTGAGACTGTTTAGCCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggc                                                                                                                                                                                             <  2:117450‑M1/251‑28 (MQ=255)
                                   aTGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttc                                                                                                                                                                                    <  1:340875‑M1/251‑37 (MQ=255)
                                       aTTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgttt                                                                                                                                                                                >  1:202412‑M1/1‑211 (MQ=255)
                                       aTTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCAAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgttt                                                                                                                                                                                >  1:207804‑M1/1‑211 (MQ=255)
                                                      actcgctCAAATATTGTTAAATTTCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACCGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattcaactacgcttcagcgcagttttcgtttacaggatgttgaa                                                                                                                                                                   >  2:373889‑M1/4‑196 (MQ=255)
                                                                 tcttgtTAAATTGCAGGCTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGACAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTAATGTAGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctg                                                                                                                                                      >  2:146605‑M1/3‑185 (MQ=255)
                                                                            cgcCGGGTTTGTATCAATTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcggttttcgtttacagggtgttgaatgggaaaattctcgggcaaaaaa                                                                                                                                            >  2:32228‑M1/2‑174 (MQ=255)
                                                                                                                 aGGGGACAGCCAACAGCCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatggga                                                                                                      <  1:31274‑M1/251‑115 (MQ=255)
                                                                                                                               aGGCATTGTGTGCAATCACCTTAGCGTTCAGATCCAGTCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACAGACATCTTTAGTACTGAGCCTGTTTAACCTAGGGAATTATTGCCGTAccatcaacgattaaactgctctccggcgcagttttcgtttacaggatgttgaacgggaaaattctggggcaaaaaaagcccgccagttactgcgggaaacctcatcctatgggagaacaatgaataat                                                                                        <  2:140653‑M1/251‑129 (MQ=255)
                                                                                                                                   aTTGGGTGCAATCACCTTAGCGTTCAGTTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaa                                                                                    <  1:231299‑M1/251‑133 (MQ=255)
                                                                                                                                          gCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcgg                                                                              <  1:384301‑M1/250‑139 (MQ=255)
                                                                                                                                            aaTCTCCTTGGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTAACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggt                                                                           <  2:176370‑M1/251‑142 (MQ=255)
                                                                                                                                                ccctTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACGGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatc                                                                       <  2:280732‑M1/250‑146 (MQ=255)
                                                                                                                                                                aCATGCGGAATGTAAAAAAGGCCGCGAGCGCCGCCCCTTCACATTCATCGATAGTACTGAGACTGTTTAACCTAGGTAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcttttacaggattttgaaagggaaaattctggggcaaaaaaatcccgccagttacggcgggaaacctcatcctatgggagaacaatgattaatgaaattgcggggttatcatctcccagtttatc                                                        <  2:413416‑M1/250‑161 (MQ=255)
                                                                                                                                                                                         gAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccataaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactc                              <  2:324096‑M1/251‑187 (MQ=255)
                                                                                                                                                                                                gCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttc                                                                                                                                                                                                <  2:336230‑M1/82‑25 (MQ=255)
                                                                                                                                                                                                gCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttc                                                                                                                                                                                                >  1:336230‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                  cgcttctaaaacatctttagtacaTAGACTGTTTAACCTAGGGAATTATTATCGTAccatcaccgattaaactgctcttcggcgcagttttctttttcaggatgtttaaaggaaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcagcctatgggagaacaatgaataatgaaattgctgggttatcatctcccagtatatccatactaacactaaggttatttactcaaaaaggca                     <  2:231072‑M1/241‑196 (MQ=255)
                                                                                                                                                                                                                  ttAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatggtagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacattttgt         <  1:379053‑M1/247‑208 (MQ=255)
                                                                                                                                                                                                                     gTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacattttgttttgtgc  <  1:120741‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                           
AAGTTAGTTATGTATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGA  >  NC_000913/1196617‑1197080

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