Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,494,388 Δ1 bp coding (684/924 nt) ydcI ← DNA‑binding transcriptional dual regulator YdcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,494,3880C.95.7% 96.6 / ‑2.6 23coding (684/924 nt)ydcIDNA‑binding transcriptional dual regulator YdcI
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (11/11);  total (11/12)
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.96e-01

TTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACGATTCAAGAAACAGCAGTTCGTAATTAAGCCCGGTCATCAGTTCAGGATCTGACATCCGAC  >  NC_000913/1494165‑1494605
                                                                                                                                                                                                                               |                                                                                                                                                                                                                         
ttCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTAGCATCTACGCGGGTCAGTATTCCAATCGGTTCGACTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGc                                                                                                                                                                                                >  1:124701/1‑250 (MQ=255)
  cGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCAc                                                                                                                                                                                              <  2:97786/250‑1 (MQ=255)
   gAATAGCGTTAATCATAAACTGGCAACCAGAAGAGAACGTCGCATATACGAGAAACCGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCCCCAGCGTGACATGACGCAGGACGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGATCTGAACt                                                                                                                                                                                             >  1:148994/1‑250 (MQ=255)
            tAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAAGCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCt                                                                                                                                                                                    <  1:15189/250‑1 (MQ=255)
               tcatcaACTGGCAACCATAAGATAACGTCGCATCTACGCGGGGCAGTATTCCAATCGCTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCATCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGCGACAAACCACACGTAATCGTATTCAACCTTAAGTTGACGAGATAGCTAAGCAGACAGCTTTTCGCTACA‑TCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGaa                                                                                                                                                                                 <  1:26864/250‑1 (MQ=255)
                catcaACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATg                                                                                                                                                                               <  1:157482/251‑1 (MQ=255)
                   caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCtg                                                                                                                                                                            >  1:169641/1‑251 (MQ=255)
                                                  aCGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACa                                                                                                                                                <  1:76222/248‑1 (MQ=255)
                                                                    aTCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGc                                                                                                                              >  1:127201/1‑248 (MQ=255)
                                                                    aTCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGc                                                                                                                              <  2:127201/248‑1 (MQ=255)
                                                                                    cATGTCCCGAAGCAGGAAGGCCCACCCACGTGGCATGACACAGGTCGTCTTTTACAGAAACAGAAGGGACAAACCCCACGACATCGGATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAACGTTTCGATACAACC‑GCAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGACTGCTGGCGTGGGGCCGTGCCTTCTGGTGATAAAACGACCGGCCATTCCCGCACCCGGCGTAGCGTTACGt                                                                                                           >  2:111043/1‑251 (MQ=255)
                                                                                           cGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGTCAAACCACACGTAAGCGTATTCAACCGTAAGTTGACTAGATAGCGAAGCAGACAGCGTTGCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGTTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTg                                                                                                    <  1:104282/251‑1 (MQ=255)
                                                                                                    caggGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGACAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGACGTAGCgg                                                                                           >  1:103666/1‑251 (MQ=255)
                                                                                                       ggCCACCAGCGTGGCATGACGCAGGTCGTATTTTACAGCGCCAGAAGGGACAAACCAAACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTg                                                                                        >  1:100282/1‑251 (MQ=255)
                                                                                                                cGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCg                                                                               <  1:204068/251‑1 (MQ=255)
                                                                                                                    gCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGCCAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCTGACa                                                                           <  2:176257/251‑1 (MQ=255)
                                                                                                                     cATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACaa                                                                          <  1:138841/251‑1 (MQ=255)
                                                                                                                              ggTCGTCTTTTACAGCGCCAGAAGGGACAAATCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTc                                                                  >  2:195334/1‑250 (MQ=255)
                                                                                                                                 cgtcTTTTACAGCGCCAGAAGGGACAAATCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACg                                                              <  1:195334/251‑1 (MQ=255)
                                                                                                                                     ttttACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACGATTc                                                          >  2:40500/1‑251 (MQ=255)
                                                                                                                                           aGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACGATTCAAGaaa                                                    <  2:56950/251‑1 (MQ=255)
                                                                                                                                                                                    tAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGACTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTTGGGCGGACAACCAGCTTCAACGATTCACGAAACAGCAGTTCGTAATTAAGCCCGCTCATCAGTTCAGGATCTg           >  2:64761/1‑251 (MQ=255)
                                                                                                                                                                                             gagaTAGCGAAGCAGACAGCGTTTCGATACAACC‑GAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACGATTCAAGAAACAGCAGTTCGTAATTAAGCCCGGTCATCAGTTCAGGATCTGACATCCGAc  >  1:206043/1‑251 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                         
TTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGAATGCTGGCGTGGCGCAGTGCCTTCTGGTGATACAACGACCGGCCATTCCAGCACCCGGCTTAGCGTTACGTTCTCCTGAAGTAGCGGGTGATTAGGGCGGACAACCAGCTTCAACGATTCAAGAAACAGCAGTTCGTAATTAAGCCCGGTCATCAGTTCAGGATCTGACATCCGAC  >  NC_000913/1494165‑1494605

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