Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,1780GA48.0% 3.0 / 24.1 25A66T (GCG→ACG) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (3/10);  new base A (12/0);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCC  >  NC_000913/3998029‑3998316
                                                                                                                                                     |                                                                                                                                             
aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgc                                                                                                                                              <  1:246082/151‑1 (MQ=255)
aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCCCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgc                                                                                                                                              <  2:697429/151‑1 (MQ=255)
      tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGGAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGGTGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGTTTTACTAAAAAacggcgg                                                                                                                                         >  1:134311/1‑150 (MQ=255)
          cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgag                                                                                                                                    <  1:215868/151‑1 (MQ=255)
           gcgcCACGCAGCAACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgaga                                                                                                                                   <  2:28512/151‑1 (MQ=255)
                     gcaACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCg                                                                                                                         <  1:113822/151‑1 (MQ=255)
                      caACCTTCTGGTGCTGGCGGTCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGAGGTTTCCAACAAAAaggcggcgggg                                                                                                                                    >  2:727435/1‑137 (MQ=255)
                                       cggtcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGCGGTTTCCAAAAAAACCGCGGCGGAGGTGGGGCATCGTATCGGGGCAACT‑‑ATggg                                                                                                     >  1:156534/1‑151 (MQ=255)
                                                         aGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACGGCGGCGGGGGTGCGTCATCGTATC‑TGGGAACTG‑ATGGGGCCGAGCCAgggggg                                                                                      >  2:136870/1‑148 (MQ=255)
                                                          gACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTAtg                                                                                  <  2:786974/151‑1 (MQ=255)
                                                             gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGGTTTCCAAAAAAACAGCGGCGGGGGGGCGTCATCGTCTG‑GGGGAAATG‑ATGGGGACGAGCCAGGGGGGGAtg                                                                                  >  2:203233/1‑148 (MQ=255)
                                                                    cTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCCACAAAACGGCGGGGGGGATGGGTCATCGTATT‑GGGGGAACG‑ATGGGGGCGAGCCCGGGGGGGATGTGGGGC‑gg                                                                         >  2:223216/1‑149 (MQ=255)
                                                                      ggTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGCGACGTTTACCAACATAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCAcc                                                                     <  1:84255/151‑1 (MQ=255)
                                                                              gTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTTCAACAAAACGGCGGCGGGGGGGGGTCCTCGGATT‑GGGGAAACG‑ATGGGGACGAGGCCGGGCGGGGTGGGGGGC‑GGCACCCTCCACgg                                                             >  2:507561/1‑151 (MQ=255)
                                                                               tATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACgg                                                             <  2:786402/150‑1 (MQ=255)
                                                                                     cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTGCCACAAAACGGCGGGGGAGATGTGTTCTAGGGTC‑GGGGTAAAG‑ATGGGGACGAGCCAGGGGGGGATGTGGGGC‑GGCGCCTTCCCcgggcggggg                                                      >  1:668205/1‑149 (MQ=255)
                                                                                              gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTtg                                             <  2:808675/151‑1 (MQ=255)
                                                                                                gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCAACAAAACGGCGGCGGGGAAGGGTCATCGTATC‑CGGGAACTG‑ATGGGGGCGAGCCACGGCGGGATGTGGGTC‑GGGACCTTCCCCCGGCGGGCGCAc                                                   >  2:856453/1‑143 (MQ=255)
                                                                                                  gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                                         <  2:409919/151‑1 (MQ=255)
                                                                                                   tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAGGCTTATCAACAAAACCGCCGCGGAGATGGGTCATCGGAAA‑GGGCGACTGAAGGGGGACG‑GCCAGGGGGGTGGGAGGGGG‑GGGCACTTCCACCGGGCGGCGCACCCttttgct                                          >  1:828143/1‑145 (MQ=255)
                                                                                                    ggAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGGTTTCCAACAAAACGGCGGGGGGGATGAGTCATCGTCTC‑GGGCGAATG‑ATGGGGACGAGCCAGGGCGGGATGTGGGGC‑GGGACCCTTCCCGGGCGGGCGCAcc                                                  >  2:793205/1‑140 (MQ=255)
                                                                                                            gCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAcc                                <  1:127964/150‑1 (MQ=255)
                                                                                                             cTCGCCATACTCGATTATGGCGGTGACGTTTTCCAACAAAACCGCGGCGGGGATGCGTCATCGTATC‑GGGGGAATG‑ATGGGGACGAGCCAGGGGGGGATGTGGGTC‑GGCACCCTCCACGGGGGGGCGCACCCTTTTCTGCGTGCGCAcc                                >  2:569291/1‑149 (MQ=255)
                                                                                                              tCGCCATACTCGATTATGGCGGTGAACTTTTTCAACAAAAAGGCGGGGGGGGTGCGTTCTTGTATT‑GGGGGACTG‑ATGGGGACCGGGCCGGGGGGGATGTGGGGCGGGCCCCTCCCCCGGGC‑GGGGCCCCGttttgttcggggggccc                                >  1:220914/1‑134 (MQ=255)
                                                                                                                              tGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGcc             >  1:594394/1‑151 (MQ=255)
                                                                                                                                   gTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAgg        <  1:690366/151‑1 (MQ=255)
                                                                                                                                     gACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAt      >  1:710619/1‑151 (MQ=255)
                                                                                                                                          ttACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGAC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACTCCAATCTGCCGCAGGATTTcc  >  2:375288/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATC‑GGGCAACTG‑ATGGGCACGAGCCAGGGCGGTATGTGGGTC‑GGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCC  >  NC_000913/3998029‑3998316

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