Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,531,598 A→T 20.9% intergenic (‑1094/‑218) yncO ← / → ydcC protein YncO/H repeat‑associated putative transposase YdcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,531,5980AT20.9% 37.2 / 2.5 19intergenic (‑1094/‑218)yncO/ydcCprotein YncO/H repeat‑associated putative transposase YdcC
Reads supporting (aligned to +/- strand):  ref base A (9/6);  new base T (2/2);  total (11/8)
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.72e-01

TCATTCAGGCAGAAGAAAGAAGATAAATCAGCAGAAGGTGTAAGCATCCGATATTATATTTCATCAAAGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACA  >  NC_000913/1531379‑1531779
                                                                                                                                                                                                                           |                                                                                                                                                                                     
tCATTCAGGCAGAAGAAAGAAGATAAATCAGCAGAAGGTGTAAGCATCCGATATTATATTTCATCAAAGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGAACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAaaagaa                                                                                                                                >  1:67350/1‑275 (MQ=255)
                             agcagAAGGTGTAAGCATCCGATATTATATTTCATCAAAGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCGCTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAaaagaaag                                                                                                                              >  1:152147/1‑248 (MQ=255)
                                    ggTGTATGCATCCGATATTATATTTCATCAAAGGATATGGCTGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAAGGAATGAAGATGCCAGCCGGATAGGAAGAGTAAACGCAACATAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGAGAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAGCACAGGGAAAGgtcaaaaga                                                                                            >  1:119883/1‑270 (MQ=255)
                                        tAAGCATCCGATATTATATTTCATCAAAGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCGGATCTAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGtt                                                                                        >  1:124816/1‑275 (MQ=255)
                                                    atcatattTCATCAAAGGAGATGGTTGGTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTACATTGGGTGTAAGATGTAAAAATGAATGAAGATGCCGGCCGGATAAGAAGAGGAGACGCAGCCTAAATAATATCTGGAATAAagaaga                                                                                                                                                                                       <  2:59592/165‑1 (MQ=255)
                                                                                        tggaaaaggCTATCAGAGCGAACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGata                                        <  2:67350/267‑1 (MQ=255)
                                                                                               tGCTATCAGAGCGCACTGGCTGATCGAGAACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCCGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAACGACAGTAAGGGTGAAGAGCAAAAGAAAGCAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTCTACTAAACAATTCGCCTTTACTGTTTATAAATTGAGATATTCCtta                                 >  1:240050/1‑275 (MQ=255)
                                                                                                          cgcACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTa                                                                                                                  >  1:178000/1‑183 (MQ=255)
                                                                                                          cgcACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGAAAAGATCATCAGAAGTTCa                                                                                      >  1:111873/1‑211 (MQ=255)
                                                                                                              cTGGCTTATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGTAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAaca                                                                                                              <  2:178000/183‑1 (MQ=255)
                                                                                                              cTGGCTGATCGAGCAACGGATTCATTGGGTGTTAGATGTAANAATGAATGAAGATGACAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACGGAATTTATTAAGATATTGCAAAGACATTAAGGGTGAAGAGTAAAAGAAAGAAGGCTGCGTTAAACACAAGGAAAGACCATCCGAAGGTCACttt                                                                                  <  2:111873/211‑1 (MQ=255)
                                                                                                                              aGTATTCATTGGGTGTAAGATGGAAAAATGAATGAAGATGCCAGCCGGGTAAGAAGAGGAAACGCAGCCTAAATAATATCTTGAATAAAGAAGATGGTACTGAAGTTATTAAGAGATTTCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTGATTACCTAAAGCTGTTTTTTATTGCTTaca    <  2:119883/273‑1 (MQ=255)
                                                                                                                                                                                     gagGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTGAAACACAGGCAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCt                           >  1:40592/1‑195 (MQ=255)
                                                                                                                                                                                         aaaCGCAGCCTAAATAATATCTGGAATAAAAAAGATGGAACTGTATTTATTAAGGGATTGCAAAGACCTTAAAGGTGAAGAGGAAAAGAAAGAAGGGTTCGTTAAACACAGGGAAAGAGCATCAGAAGTTCACTTTTTGTACAAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCGTATTACCTAAAg                       <  2:40592/195‑1 (MQ=37)
                                                                                                                                                                                                  cTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTc                                     >  1:19207/1‑172 (MQ=255)
                                                                                                                                                                                                      ataataTCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCtta                                 <  2:19207/172‑1 (MQ=255)
                                                                                                                                                                                                               ggAATAAAGAAGTTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTa      >  1:19428/1‑190 (MQ=255)
                                                                                                                                                                                                               ggAATAAAGAAGTTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTa      >  1:19421/1‑190 (MQ=255)
                                                                                                                                                                                                                   tAAAGAAGTTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTAcaca  <  2:19421/190‑1 (MQ=255)
                                                                                                                                                                                                                   tAAAGAAGTTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTAcaca  <  2:19428/190‑1 (MQ=255)
                                                                                                                                                                                                                      agaagaTGGCATTGAATTTATTAAGAGATTGCAAAGAAATTAAGGGTGAAGAGGAAGAGAAAGAAGGCTGCGTTATACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGACTTTTATGTTTAACAATTGAGCTATTCCTTATTACCTAAAGCTGttt                 >  2:245242/1‑172 (MQ=255)
                                                                                                                                                                                                                          gaTGGCACTGAATTTATTAAGAGAATGCAAAGACATGAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTAGTACCTAAAGCTGTTTTTTa             <  1:245242/172‑1 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                     
TCATTCAGGCAGAAGAAAGAAGATAAATCAGCAGAAGGTGTAAGCATCCGATATTATATTTCATCAAAGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACA  >  NC_000913/1531379‑1531779

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