Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,306,1360AC52.2% 32.1 / 37.3 23G575G (GGT→GGG) deaDATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (11/0);  new base C (0/12);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 7.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05

GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA  >  NC_000913/3305986‑3306282
                                                                                                                                                        |                                                                                                                                                  
gAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGa                                                                                                                                                    <  1:401220/151‑1 (MQ=255)
   acgacgacgACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGacc                                                                                                                                                  >  2:302982/1‑150 (MQ=255)
                          caCGACGGGGGGCGCGGCCTCCACGACGTTCGCCGCTAAAACGACGACCATCCCCACGGCCCCCTTCACGGCGTTCCCCGCTGAAGTT‑ACGCCCGCCTTCACGACGTTCGCCACCAAAA‑CCCCGCCCACCGCCACGACGCTCACCGCCa                                                                                                                            <  2:109701/149‑1 (MQ=255)
                          caCGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGt                                                                                                                          >  2:474356/1‑151 (MQ=255)
                                           cctccccgactTCCGCCGCAAAACCGCCGCCCACCCCCCCGCCCCCCTCCCCGGCGTCCCCCGCTAAAGTT‑ACGACCCCCTCCCCGCCGTCCCCCCCCGAAC‑CCCCGCCCACCGCCACGACGCTCACCGCCAGTATGCGGCtgtg                                                                                                               <  1:711037/139‑1 (MQ=255)
                                                   agttcccccgcAAAACCGCCTACCATCCCCCCGGCCACCTCCCGGGGGTTCCCCGGAGAATTT‑CCGCCCCCCTCCACGAC‑TTCCCCCCCAAAACCCCCGCCCGCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTaa                                                                                                  <  1:543883/143‑1 (MQ=255)
                                                    cgTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACt                                                                                                >  1:646103/1‑151 (MQ=255)
                                                     gTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTg                                                                                               >  1:357946/1‑151 (MQ=255)
                                                          ccgcTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATg                                                                                           >  2:74806/1‑150 (MQ=255)
                                                          ccgaaaaaaCGCCCACCACCCCCCCGCCCCCCCTCGCGGCTTCCCCCGCTAAGGTT‑CCGCCCCCCTTCCCGCCGTCCCCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCa                                                                                             <  2:897004/143‑1 (MQ=255)
                                                                 aaCGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCgg                                                                                   >  1:364087/1‑151 (MQ=255)
                                                                     acgacCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCtt                                                                                >  1:99355/1‑150 (MQ=255)
                                                                            tCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTgaga                                                                         >  2:769235/1‑150 (MQ=255)
                                                                                  cGGCCCCCCCCACGGCGTCCCCCGTTAAAGTT‑CCGCCCCCCTTAACGACGTCCCCCCCCAAAC‑CCGCCCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCg                                                                    <  2:604324/149‑1 (MQ=255)
                                                                                      cccctccccGGCGTCCCCCGCT‑AAGTTGGCGCCCCCCTCCACGACG‑TCCCCCCCGAAACCCCCGCCCCCCGCCCCGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGt                                                                 <  2:760428/140‑1 (MQ=255)
                                                                                             aCGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTaa                                                         >  2:57093/1‑149 (MQ=255)
                                                                                                      cccgtgaaattt‑ACGCCCCCCTCCCCGCCTTTCGCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcagc                                              <  1:132022/141‑1 (MQ=255)
                                                                                                      aCCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcag                                               >  2:938620/1‑150 (MQ=255)
                                                                                                                 ttccgccccccTCCACGCCGTCCCCCCCCAAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCAcc                                      <  1:583482/140‑1 (MQ=255)
                                                                                                                                   cgTTCCCCCCCGAAACCCACGCCCCCCCCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGc                      <  1:851104/148‑1 (MQ=255)
                                                                                                                                         cccccGAAC‑CCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAg                    <  2:323043/140‑1 (MQ=255)
                                                                                                                                            ccGAAC‑CCACGCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggtg          <  1:290139/150‑1 (MQ=255)
                                                                                                                                              gAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggagg         >  1:945491/1‑149 (MQ=255)
                                                                                                                                                   ccACGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgag   <  1:474438/151‑1 (MQ=255)
                                                                                                                                                     aCGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgaga  >  2:357787/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTT‑ACGACCGCCTTCACGACGTTCGCCACCGAAA‑CCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA  >  NC_000913/3305986‑3306282

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