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

AGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑AAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTC  >  NC_000913/1531446‑1531931
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aGGATATTGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCAc                                                                                                                                                                                                                                              <  1:230227/251‑1 (MQ=255)
           gatAAAGAATGTGCACAAGCTATCAGAGCGCACTGGCTGATCGAGCACATTCTGCATGGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTTGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAATACATTAAGGGTGAAGAGGAAAAGAAAGAAAGCTGCTTTAAACACAGGGAAAGATCATCAG‑AATTTCGCTTTTTGTACTa                                                                                                                                                                                                                                   <  2:168178/249‑1 (MQ=255)
                   aTTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTc                                                                                                                                                                                                                           >  1:5268/1‑251 (MQ=255)
                                aTCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGttt                                                                                                                                                                                                              <  1:48813/251‑1 (MQ=255)
                                         cACTGGCTGATCGCGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGAACTAAACAAt                                                                                                                                                                                                                             >  1:88338/1‑227 (MQ=255)
                                                                                                     aGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTAc                                                                                                                                         >  1:216047/1‑251 (MQ=255)
                                                                                                      gCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTCAGAGATTGCAAAGACATTAAGGGCGTAGAGGAAAAGTAAGAAGGCTGCGTTAACCACAGGGAAAGCTCATCAG‑GAGTTCACTTTTTGAACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAGGCTGTTTTTTATTGCTTACACATGATCAACTGCTCCTTAc                                                                                                                                         >  1:269104/1‑250 (MQ=37)
                                                                                                      gCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        <  1:323739/251‑1 (MQ=255)
                                                                                                                    ggAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGAAGAAc                                                                                                                          >  1:301635/1‑251 (MQ=255)
                                                                                                                    ggAAACGCAGACGAAATAAGAGCTGGAAGTAAGGATATGGTACTGAATTTATTAAGAGATTGCAAAGACACTAAGGGTGGAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAATGATCATCAGTATGTT‑GCTTTTTGTACTAAACAATTCGCGTTGTATGTTTAAAAATTGTGATATTCCTTATTACCTAAAGAGGTTTTTTATTGCTTACTCATGTTCAAATACTCCTTACATAATTGAGGAGaag                                                                                                                          <  2:77540/251‑2 (MQ=255)
                                                                                                                                     aataTCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTaaaaa                                                                                                         <  1:194895/251‑1 (MQ=255)
                                                                                                                                                        aTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAAt                                                                                                                     >  2:110353/1‑220 (MQ=255)
                                                                                                                                                        aTGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAAt                                                                                                                     <  1:110353/220‑1 (MQ=255)
                                                                                                                                                                  aTTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        >  1:286028/1‑191 (MQ=255)
                                                                                                                                                                  aTTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACa                                                                                                                                        <  2:286028/191‑1 (MQ=255)
                                                                                                                                                                              aTTGCAAAGACATTAAAGGTGAAGAGGACAAGTAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGCTATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATT‑TCACACGATTACAGACaa                                                                >  2:127291/1‑250 (MQ=18)
                                                                                                                                                                                  acgacgcccttccgaactagagGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCTCTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAGTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTg                                                           <  2:25366/233‑1 (MQ=18)
                                                                                                                                                                                               ggTGAAGAGGAAAAGAAAGAATGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGTTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACAt                                              <  1:156257/251‑1 (MQ=255)
                                                                                                                                                                                                    agagGAAAAGAAAGAAGGCTGCGTTAAACACATGGAAAGATCATCAG‑GAGTTCACTTTTTTTACTAAACAATTTGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGTGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGtt                                         <  2:5695/251‑1 (MQ=255)
                                                                                                                                                                                                    agagGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTgagatcggaagagcaacgt                                         >  1:25366/1‑233 (MQ=255)
                                                                                                                                                                                                           aagaaagaAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGat                                  >  1:48062/1‑251 (MQ=255)
                                                                                                                                                                                                                      cTGCGTTAAACACAGGGAAAGATCATCAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGtt                       >  1:115762/1‑251 (MQ=255)
                                                                                                                                                                                                                                           atcatcAG‑GAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTc  <  1:193008/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AGGATATGGATGCTAAAGAATTTGCACATGCTATCAGAGCGCACTGGCTGATCGAGCACAGTCTTCATTGGGTGTTAGATGTAAAAATGAATGAAGATGCCAGCCGGATAAGAAGAGGAAACGCAGCCTAAATAATATCTGGAATAAAGAAGATGGCACTGAATTTATTAAGAGATTGCAAAGACATTAAGGGTGAAGAGGAAAAGAAAGAAGGCTGCGTTAAACACAGGGAAAGATCATCAG‑AAGTTCACTTTTTGTACTAAACAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTAAAGCTGTTTTTTATTGCTTACACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATC‑CCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTC  >  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: ‑