Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,1780GA42.9% 46.9 / 23.4 42A66T (GCG→ACG) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (6/18);  new base A (18/0);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 5.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  NC_000913/3998030‑3998328
                                                                                                                                                     |                                                                                                                                                        
aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAgcg                                                                                                                                                        <  2:864295/151‑1 (MQ=255)
      ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAgcggcggc                                                                                                                                                   <  1:474322/150‑1 (MQ=255)
        ccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGga                                                                                                                                                <  1:323792/151‑1 (MQ=255)
         cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGGGTAAGGTACCCACAAAAcggcggggggg                                                                                                                                               >  1:164052/1‑149 (MQ=255)
                cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGGAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGTCGGTTTCCAACAAAacagcgggcggg                                                                                                                                               >  1:609422/1‑139 (MQ=255)
                   agcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATc                                                                                                                                    <  2:950267/151‑1 (MQ=255)
                            cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGAGGTTTA‑CAAAAAAAAGGCGGCGGGGGGGCG‑GCCTCGGA‑TCgggg                                                                                                                          >  2:951025/1‑150 (MQ=255)
                              ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCggg                                                                                                                           <  2:894697/148‑1 (MQ=255)
                                 gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACt                                                                                                                      <  2:295901/150‑1 (MQ=255)
                                    ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGAAGGTTTCCCACAAAACGGCGGGGGGGGTGCG‑TTCTCGTA‑TCCGGGAAatgatg                                                                                                                  >  2:146995/1‑151 (MQ=255)
                                       gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCAACAAAACGGCGGCGGGGATGGG‑GCATCGTA‑TCGGGCAACTGATggg                                                                                                                >  2:548709/1‑150 (MQ=255)
                                       gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAACGGCGGCGGGGATGGG‑GCCTCGTA‑TCGGGCAACAGAAggg                                                                                                                >  2:624164/1‑150 (MQ=255)
                                           gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGa                                                                                                           <  1:945603/151‑1 (MQ=255)
                                             gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGc                                                                                                         <  2:20862/151‑1 (MQ=255)
                                                      tAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGAAGGTTTCCCACAAAACGGCGGGGGGGATGGG‑GCATCGTA‑TCGGGGAACAGATGGGGACGAGCCCgggggg                                                                                                 >  2:293253/1‑150 (MQ=255)
                                                       aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGt                                                                                                <  2:949435/150‑1 (MQ=255)
                                                       aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTACCAAAAAAAAAGCGGCGGGGAATCGCGCTTCTAA‑TGGGCCAAC‑GATGGGGACGAGCCAggg                                                                                                    >  2:966092/1‑146 (MQ=255)
                                                          aCGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTTTCAACAAAACGGGGGGGGGGGAGGG‑GCATCGTATTGGGGGAACT‑ATGGGGACGAGGCCggggggggtg                                                                                             >  1:635818/1‑146 (MQ=255)
                                                               cgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTc                                                                                       <  2:777745/151‑1 (MQ=255)
                                                                         cATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc                                                                             <  1:416662/151‑1 (MQ=255)
                                                                               aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg                                                                        <  1:868320/150‑1 (MQ=255)
                                                                                   ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCACCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                                                   <  1:50026/151‑1 (MQ=255)
                                                                                   ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                                                   <  2:86030/151‑1 (MQ=255)
                                                                                    cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCCACAAAACGGCGGCGGGGGAGCG‑TCCTCGTA‑TCGGGCCACTGATGGGGACGGGCCCGGGGGGGGTGGGGGTCGGCACCTTCCCCCGGCCGgcg                                                                  >  2:828444/1‑151 (MQ=255)
                                                                                        ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCAACAAAACGGCGGCGGAGGTGCG‑TCATCGTA‑TCCGGCCACCCATGGGGGCGAGCCAAGGGGGGATGGTGGTCGGGCCCCTCCCCCGGGTGGCGcccc                                                              >  2:469126/1‑151 (MQ=255)
                                                                                         tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTATCCACAAAAAGGCGGCGGGGGGTGG‑TCCTTTTA‑TCGGGGGACCGAAGGGCACGAGCCCGGGGGGGTTGTGGGGCTGGACCCTcccccgggcg                                                                     >  2:739716/1‑135 (MQ=255)
                                                                                             gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGttt                                                          <  1:441477/150‑1 (MQ=255)
                                                                                               gCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTTACGCTTTCCAACCAAAAGGCGGCGGAGCGGCG‑CCCCCGTA‑ACGGGGGACTGAAGGGGGCGAGGCCGGGGGGGGTGTTGGGCGGCACCCTCCCCGGGCCGGGGCAcc                                                              >  1:406234/1‑144 (MQ=255)
                                                                                                 gTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTTCAACAAAACGGCGGCGGGGATGCG‑TCATCGTA‑TCGGGGAACCGAAGGGGACGAGGCAGGGGGGGATGTGGGGCGGCACCCTCCACCGGCCGGCGCACCGgtt                                                          >  2:500296/1‑146 (MQ=255)
                                                                                                    gAAAACTGCTCGCCATACTCGATTATGGC‑GGTGCCGGTTATCAACAAAACGGCGGCGGGGGAGGG‑GTATCGTA‑ATGGGGAACTGATGGGGACGAGGCCGGGGGGGGTGGGGGGGGGGGCCTTCCACCGGCCGgcggccccgttggt                                                       >  1:666113/1‑135 (MQ=255)
                                                                                                          tGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAACGGCGGCGGGGATGCG‑GCATCGTA‑TCGGGCCACTGATGGGGACGAGCCCGGGCGGGATGTGGGGCGGCACCCTCCCCGGGGTGGCGCACCGTTTTCTGCGTGCGCAcc                                            >  2:21435/1‑151 (MQ=255)
                                                                                                            cTCGCCATACTCGATTATGGC‑GGTGACGTTTTTCAAAAAAACCGCCGCGGGGGAGCG‑TCATCGTA‑TCTGGCGACTGAATGGGGCGGGCCAGGGGGGGGTGTGGGGTGGGGCCCTCCCCCGGCCGGGGCCCCCGTTGGTGCGGGCGcccc                                            >  1:81483/1‑149 (MQ=255)
                                                                                                              cGCCATACTCGATTATGGC‑GGTGGCGTTTTACAACAAAACGGGGGCGGGGGAGGG‑TCATCGTA‑TCGGGGAAAAGATGGGGACGAGCCCGGGGGGGATGTGGGGGGGGACCCTCCCCCGGCCGGCGCCCCGgtt                                                          >  1:493495/1‑133 (MQ=255)
                                                                                                              cGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACAata                                         <  2:9001/150‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc                                    <  2:772396/151‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc                                    <  2:275455/151‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGC‑GGTGAAGTTTTTCCAAAAAAAGGCGGGGGGGAAGCG‑TCCTCGTA‑TCGGGCAACCGATGGGGACGAGCCCGGGGGGGATGTGGGGCGGGGCCCTCCACCGGCCGGGGCACCCCTTTCtgcgggggcccc                                            >  1:47312/1‑135 (MQ=255)
                                                                                                                    aCTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc                                  <  1:641625/151‑1 (MQ=255)
                                                                                                                            aTGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTg                           >  2:250347/1‑150 (MQ=255)
                                                                                                                             tGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGCCTCGCGCACCGTTTTCTGCGTGCGCGCCATATGGACGCCAATCTGc                          >  2:24931/1‑150 (MQ=255)
                                                                                                                                c‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGc                       >  2:534165/1‑150 (MQ=255)
                                                                                                                                        gTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc               <  2:847344/151‑1 (MQ=255)
                                                                                                                                             ccAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc              >  2:931969/1‑147 (MQ=255)
                                                                                                                                              caacaaAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc              >  2:662790/1‑146 (MQ=255)
                                                                                                                                               aacaaAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc        >  1:777708/1‑151 (MQ=255)
                                                                                                                                                     gcggcggcgGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa  <  1:406395/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  NC_000913/3998030‑3998328

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