Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,531,689 A→G intergenic (‑1185/‑127) 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,6890AG91.3% 57.8 / ‑2.7 23intergenic (‑1185/‑127)yncO/ydcCprotein YncO/H repeat‑associated putative transposase YdcC
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (11/10);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTC  >  NC_000913/1531446‑1531931
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aGGATATTGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCAc                                                                                                                                                                                                                                              <  1:230227/251‑1 (MQ=255)
           gatAAAGAATGTGCACAAGCTATCAGAGCGCACTGGCTGATCGAGCACATTCTGCATGGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTTGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAATACATTAAGGGTGAAGAGGAAAAGAAAGAAAGCTGCTTTAAACACAGGGAAAGATCATCAGAATTTCGCTTTTTGTACTa                                                                                                                                                                                                                                   <  2:168178/249‑1 (MQ=255)
                   aTTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTc                                                                                                                                                                                                                           >  1:5268/1‑251 (MQ=255)
                                aTCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGttt                                                                                                                                                                                                              <  1:48813/251‑1 (MQ=255)
                                         cACTGGCTGATCGCGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGAACTAAACAAt                                                                                                                                                                                                                             >  1:88338/1‑227 (MQ=255)
                                                                                                     aGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTAc                                                                                                                                         >  1:216047/1‑251 (MQ=255)
                                                                                                      gCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTCAGAGATTGCAAAGACATTAAGGGCGTAGAGGAAAAGTAAGAAGGCTGCGTTAACCACAGGGAAAGCTCATCAGGAGTTCACTTTTTGAACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAGGCTGTTTTTTATTGCTTACACATGATCAACTGCTCCTTAc                                                                                                                                         >  1:269104/1‑250 (MQ=37)
                                                                                                      gCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        <  1:323739/251‑1 (MQ=255)
                                                                                                                    ggAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGAAGAAc                                                                                                                          >  1:301635/1‑251 (MQ=255)
                                                                                                                    ggAAACGCAGACGAAATAAGAGCTGGAAGTAAGGATATGGTACTGAATTTATTAAGAGATTGCAAAGACACTAAGGGTGGAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAATGATCATCAGTATGTTGCTTTTTGTACTAAACAATTCGCGTTGTATGTTTAAAAATTGTGATATTCCTTATTACCTAAAGAGGTTTTTTATTGCTTACTCATGTTCAAATACTCCTTACATAATTGAGGAGaag                                                                                                                          <  2:77540/251‑2 (MQ=255)
                                                                                                                                     aataTCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTaaaaa                                                                                                         <  1:194895/251‑1 (MQ=255)
                                                                                                                                                        aTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAAt                                                                                                                     >  2:110353/1‑220 (MQ=255)
                                                                                                                                                        aTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAAt                                                                                                                     <  1:110353/220‑1 (MQ=255)
                                                                                                                                                                  aTTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        >  1:286028/1‑191 (MQ=255)
                                                                                                                                                                  aTTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        <  2:286028/191‑1 (MQ=255)
                                                                                                                                                                              aTTGCAAAGACATTAAAGGTGAAGAGGACAAGTAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGCTATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTTCACACGATTACAGACaa                                                                >  2:127291/1‑250 (MQ=18)
                                                                                                                                                                                  acgacgcccttccgaactagagGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCTCTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAGTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTg                                                           <  2:25366/233‑1 (MQ=18)
                                                                                                                                                                                               ggTGAAGAGGAAAAGAAAGAATGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGTTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACAt                                              <  1:156257/251‑1 (MQ=255)
                                                                                                                                                                                                    agagGAAAAGAAAGAAGGCTGCGTTAAACACATGGAAAGATCATCAGGAGTTCACTTTTTTTACTAAACAATTTGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGTGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGtt                                         <  2:5695/251‑1 (MQ=255)
                                                                                                                                                                                                    agagGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTgagatcggaagagcaacgt                                         >  1:25366/1‑233 (MQ=255)
                                                                                                                                                                                                           aagaaagaAGGCTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGat                                  >  1:48062/1‑251 (MQ=255)
                                                                                                                                                                                                                      cTGCGTTAAACACAGGGAAAGATCATCAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGtt                       >  1:115762/1‑251 (MQ=255)
                                                                                                                                                                                                                                           atcatcAGGAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTc  <  1:193008/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAGAAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTC  >  NC_000913/1531446‑1531931

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