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

AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  NC_000913/3998030‑3998328
                                                                                                                                                     |                                                                                                                                                        
aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcg                                                                                                                                                        <  2:864295/151‑1 (MQ=255)
      ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                                                                   <  1:474322/150‑1 (MQ=255)
        ccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGga                                                                                                                                                <  1:323792/151‑1 (MQ=255)
         cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGTAAGGTACCCACAAAAcggcggggggg                                                                                                                                               >  1:164052/1‑149 (MQ=255)
                cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGGAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGTCGGTTTCCAACAAAacagcgggcggg                                                                                                                                               >  1:609422/1‑139 (MQ=255)
                   agcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc                                                                                                                                    <  2:950267/151‑1 (MQ=255)
                            cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGAGGTTTACAAAAAAAAGGCGGCGGGGGGGCGGCCTCGGATCgggg                                                                                                                          >  2:951025/1‑150 (MQ=255)
                              ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCggg                                                                                                                           <  2:894697/148‑1 (MQ=255)
                                 gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACt                                                                                                                      <  2:295901/150‑1 (MQ=255)
                                    ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGGTTTCCCACAAAACGGCGGGGGGGGTGCGTTCTCGTATCCGGGAAatgatg                                                                                                                  >  2:146995/1‑151 (MQ=255)
                                       gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTCCAACAAAACGGCGGCGGGGATGGGGCATCGTATCGGGCAACTGATggg                                                                                                                >  2:548709/1‑150 (MQ=255)
                                       gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAACGGCGGCGGGGATGGGGCCTCGTATCGGGCAACAGAAggg                                                                                                                >  2:624164/1‑150 (MQ=255)
                                           gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGa                                                                                                           <  1:945603/151‑1 (MQ=255)
                                             gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGc                                                                                                         <  2:20862/151‑1 (MQ=255)
                                                      tAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGGTTTCCCACAAAACGGCGGGGGGGATGGGGCATCGTATCGGGGAACAGATGGGGACGAGCCCgggggg                                                                                                 >  2:293253/1‑150 (MQ=255)
                                                       aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGt                                                                                                <  2:949435/150‑1 (MQ=255)
                                                       aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCAAAAAAAAAGCGGCGGGGAATCGCGCTTCTAATGGGCCAACGATGGGGACGAGCCAggg                                                                                                    >  2:966092/1‑146 (MQ=255)
                                                          aCGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTTCAACAAAACGGGGGGGGGGGAGGGGCATCGTATTGGGGGAACTATGGGGACGAGGCCggggggggtg                                                                                             >  1:635818/1‑146 (MQ=255)
                                                               cgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTc                                                                                       <  2:777745/151‑1 (MQ=255)
                                                                         cATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc                                                                             <  1:416662/151‑1 (MQ=255)
                                                                               aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg                                                                        <  1:868320/150‑1 (MQ=255)
                                                                                   ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCACCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                                                   <  1:50026/151‑1 (MQ=255)
                                                                                   ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                                                   <  2:86030/151‑1 (MQ=255)
                                                                                    cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTCCCACAAAACGGCGGCGGGGGAGCGTCCTCGTATCGGGCCACTGATGGGGACGGGCCCGGGGGGGGTGGGGGTCGGCACCTTCCCCCGGCCGgcg                                                                  >  2:828444/1‑151 (MQ=255)
                                                                                        ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTCCAACAAAACGGCGGCGGAGGTGCGTCATCGTATCCGGCCACCCATGGGGGCGAGCCAAGGGGGGATGGTGGTCGGGCCCCTCCCCCGGGTGGCGcccc                                                              >  2:469126/1‑151 (MQ=255)
                                                                                         tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTATCCACAAAAAGGCGGCGGGGGGTGGTCCTTTTATCGGGGGACCGAAGGGCACGAGCCCGGGGGGGTTGTGGGGCTGGACCCTcccccgggcg                                                                     >  2:739716/1‑135 (MQ=255)
                                                                                             gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGttt                                                          <  1:441477/150‑1 (MQ=255)
                                                                                               gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTTACGCTTTCCAACCAAAAGGCGGCGGAGCGGCGCCCCCGTAACGGGGGACTGAAGGGGGCGAGGCCGGGGGGGGTGTTGGGCGGCACCCTCCCCGGGCCGGGGCAcc                                                              >  1:406234/1‑144 (MQ=255)
                                                                                                 gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCAACAAAACGGCGGCGGGGATGCGTCATCGTATCGGGGAACCGAAGGGGACGAGGCAGGGGGGGATGTGGGGCGGCACCCTCCACCGGCCGGCGCACCGgtt                                                          >  2:500296/1‑146 (MQ=255)
                                                                                                    gAAAACTGCTCGCCATACTCGATTATGGCGGTGCCGGTTATCAACAAAACGGCGGCGGGGGAGGGGTATCGTAATGGGGAACTGATGGGGACGAGGCCGGGGGGGGTGGGGGGGGGGGCCTTCCACCGGCCGgcggccccgttggt                                                       >  1:666113/1‑135 (MQ=255)
                                                                                                          tGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAACGGCGGCGGGGATGCGGCATCGTATCGGGCCACTGATGGGGACGAGCCCGGGCGGGATGTGGGGCGGCACCCTCCCCGGGGTGGCGCACCGTTTTCTGCGTGCGCAcc                                            >  2:21435/1‑151 (MQ=255)
                                                                                                            cTCGCCATACTCGATTATGGCGGTGACGTTTTTCAAAAAAACCGCCGCGGGGGAGCGTCATCGTATCTGGCGACTGAATGGGGCGGGCCAGGGGGGGGTGTGGGGTGGGGCCCTCCCCCGGCCGGGGCCCCCGTTGGTGCGGGCGcccc                                            >  1:81483/1‑149 (MQ=255)
                                                                                                              cGCCATACTCGATTATGGCGGTGGCGTTTTACAACAAAACGGGGGCGGGGGAGGGTCATCGTATCGGGGAAAAGATGGGGACGAGCCCGGGGGGGATGTGGGGGGGGACCCTCCCCCGGCCGGCGCCCCGgtt                                                          >  1:493495/1‑133 (MQ=255)
                                                                                                              cGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACAata                                         <  2:9001/150‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc                                    <  2:772396/151‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc                                    <  2:275455/151‑1 (MQ=255)
                                                                                                                  aTACTCGATTATGGCGGTGAAGTTTTTCCAAAAAAAGGCGGGGGGGAAGCGTCCTCGTATCGGGCAACCGATGGGGACGAGCCCGGGGGGGATGTGGGGCGGGGCCCTCCACCGGCCGGGGCACCCCTTTCtgcgggggcccc                                            >  1:47312/1‑135 (MQ=255)
                                                                                                                    aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc                                  <  1:641625/151‑1 (MQ=255)
                                                                                                                            aTGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTg                           >  2:250347/1‑150 (MQ=255)
                                                                                                                             tGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGCCTCGCGCACCGTTTTCTGCGTGCGCGCCATATGGACGCCAATCTGc                          >  2:24931/1‑150 (MQ=255)
                                                                                                                                cGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGc                       >  2:534165/1‑150 (MQ=255)
                                                                                                                                        gTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc               <  2:847344/151‑1 (MQ=255)
                                                                                                                                             ccAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc              >  2:931969/1‑147 (MQ=255)
                                                                                                                                              caacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc              >  2:662790/1‑146 (MQ=255)
                                                                                                                                               aacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc        >  1:777708/1‑151 (MQ=255)
                                                                                                                                                     gcggcggcgGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa  <  1:406395/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA  >  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: