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→GGGdeaDATP‑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

GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA  >  NC_000913/3305986‑3306282
                                                                                                                                                        |                                                                                                                                                  
gAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGa                                                                                                                                                    <  1:401220/151‑1 (MQ=255)
   acgacgacgACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGacc                                                                                                                                                  >  2:302982/1‑150 (MQ=255)
                          caCGACGGGGGGCGCGGCCTCCACGACGTTCGCCGCTAAAACGACGACCATCCCCACGGCCCCCTTCACGGCGTTCCCCGCTGAAGTTACGCCCGCCTTCACGACGTTCGCCACCAAAACCCCGCCCACCGCCACGACGCTCACCGCCa                                                                                                                            <  2:109701/149‑1 (MQ=255)
                          caCGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGt                                                                                                                          >  2:474356/1‑151 (MQ=255)
                                           cctccccgactTCCGCCGCAAAACCGCCGCCCACCCCCCCGCCCCCCTCCCCGGCGTCCCCCGCTAAAGTTACGACCCCCTCCCCGCCGTCCCCCCCCGAACCCCCGCCCACCGCCACGACGCTCACCGCCAGTATGCGGCtgtg                                                                                                               <  1:711037/139‑1 (MQ=255)
                                                   agttcccccgcAAAACCGCCTACCATCCCCCCGGCCACCTCCCGGGGGTTCCCCGGAGAATTTCCGCCCCCCTCCACGACTTCCCCCCCAAAACCCCCGCCCGCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTaa                                                                                                  <  1:543883/143‑1 (MQ=255)
                                                    cgTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACt                                                                                                >  1:646103/1‑151 (MQ=255)
                                                     gTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTg                                                                                               >  1:357946/1‑151 (MQ=255)
                                                          ccgcTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATg                                                                                           >  2:74806/1‑150 (MQ=255)
                                                          ccgaaaaaaCGCCCACCACCCCCCCGCCCCCCCTCGCGGCTTCCCCCGCTAAGGTTCCGCCCCCCTTCCCGCCGTCCCCCCCCAAACCCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCa                                                                                             <  2:897004/143‑1 (MQ=255)
                                                                 aaCGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCgg                                                                                   >  1:364087/1‑151 (MQ=255)
                                                                     acgacCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCtt                                                                                >  1:99355/1‑150 (MQ=255)
                                                                            tCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTgaga                                                                         >  2:769235/1‑150 (MQ=255)
                                                                                  cGGCCCCCCCCACGGCGTCCCCCGTTAAAGTTCCGCCCCCCTTAACGACGTCCCCCCCCAAACCCGCCCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCg                                                                    <  2:604324/149‑1 (MQ=255)
                                                                                      cccctccccGGCGTCCCCCGCTAAGTTGGCGCCCCCCTCCACGACGTCCCCCCCGAAACCCCCGCCCCCCGCCCCGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGt                                                                 <  2:760428/140‑1 (MQ=255)
                                                                                             aCGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTaa                                                         >  2:57093/1‑149 (MQ=255)
                                                                                                      cccgtgaaatttACGCCCCCCTCCCCGCCTTTCGCCCCCAAACCCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcagc                                              <  1:132022/141‑1 (MQ=255)
                                                                                                      aCCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTgcag                                               >  2:938620/1‑150 (MQ=255)
                                                                                                                 ttccgccccccTCCACGCCGTCCCCCCCCAAACCCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCAcc                                      <  1:583482/140‑1 (MQ=255)
                                                                                                                                   cgTTCCCCCCCGAAACCCACGCCCCCCCCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGc                      <  1:851104/148‑1 (MQ=255)
                                                                                                                                         cccccGAACCCCCGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAg                    <  2:323043/140‑1 (MQ=255)
                                                                                                                                            ccGAACCCACGCCCCCCGCCACGACCCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggtg          <  1:290139/150‑1 (MQ=255)
                                                                                                                                              gAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGAtggagg         >  1:945491/1‑149 (MQ=255)
                                                                                                                                                   ccACGCCCCCCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgag   <  1:474438/151‑1 (MQ=255)
                                                                                                                                                     aCGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTgaga  >  2:357787/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAAACGACGACGACCGGTAGAATCATCACGACGCGGAGCGCGGCCTTCACGACGTTCGCCGCTAAAACGACGACCATCACCACGGCCACCTTCACGGCGTTCACCGCTGAAGTTACGACCGCCTTCACGACGTTCGCCACCGAAACCACGACCACCGCCACGACGCTCACCGCCAGTATGCGGCTGTGCATCGCCCAGTAACTGCATGTTCATCGGCTTGTTGAGAATGCGAGTGCGCGTAAAGTGTTGCAGCACTTCACCCGGCATACCTTTCGGCAGTTCGATGGTGGAGTGAGA  >  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: