Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,1780GA40.4% 61.8 / 26.1 47A66T (GCG→ACG) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (2/26);  new base A (19/0);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.01e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCT  >  NC_000913/3998041‑3998323
                                                                                                                                         |                                                                                                                                                  
cgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGAt                                                                                                                                       <  2:820510/151‑1 (MQ=255)
   cACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCg                                                                                                                                    <  1:678082/151‑1 (MQ=255)
       cagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGAGTTTTACAAAAAAAAGGCGGCGGGGGTGCGTcc                                                                                                                                 >  2:101042/1‑149 (MQ=255)
        agcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc                                                                                                                               <  1:266317/151‑1 (MQ=255)
                         ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATg                                                                                                             <  1:487594/151‑1 (MQ=255)
                            gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGc                                                                                                          <  1:578758/151‑1 (MQ=255)
                            gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGc                                                                                                          <  1:473206/151‑1 (MQ=255)
                              gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                                                                                                        <  1:581883/151‑1 (MQ=255)
                               cgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACg                                                                                                       <  1:929060/151‑1 (MQ=255)
                                  gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGTTTACCCACAAAACGGCGGCGGGGGGGCGTCATCGTATCGGGCGAATGATggg                                                                                                           >  1:390081/1‑144 (MQ=255)
                                      aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAgg                                                                                                <  1:729691/151‑1 (MQ=255)
                                        tGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGGTTTCCAACAAAACGGGGGCGGGGAAGGGGCACCGTATACGGCAACCGGAGGGCACGAGcc                                                                                                   >  2:264679/1‑146 (MQ=255)
                                             aGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACCTTTTCTCAAAAAACGGCGGCGGGGATGCGTCCTCGGCTCGGGGAACAGATGGGCCCGAGcccaggggg                                                                                            >  2:940861/1‑141 (MQ=255)
                                                   gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACGGCGGGGGGGGGGGGTCATCGTATCCGGGGACAGATGGGGACCGGCCAGGGGGGGGTGTggg                                                                                    >  2:118728/1‑150 (MQ=255)
                                                     gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCg                                                                                 <  1:852609/151‑1 (MQ=255)
                                                     gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTTCAACAAAACGGCGGCGGGGGAGCGGCCTCGTATTGGGGGAATGATGGGGACGAGCCAGGGGGGGGTGTggg                                                                                    >  2:397336/1‑148 (MQ=255)
                                                            tGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCAACAAAACGGCGGCGGGGAGGGGTCATCGTATCTGGGAACTGATGGGGACGAGCCAGGGGGGGGTGTGGGGCGGGAccct                                                                          >  2:856894/1‑149 (MQ=255)
                                                             gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTc                                                                         <  2:65878/151‑1 (MQ=255)
                                                             gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACACGGGCGGGGGTGCGGCCTCATCTCTGGCCACCGATGGGCGCGAGCCCGGGGGGGTTGTGGGGCggggccct                                                                          >  1:369993/1‑144 (MQ=255)
                                                                tCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCAAAAAAACGGGGGCGGGGATGCGTCATCGGATTGGGGGACCGGTGGGGACGAGCCCGGGGGGGGTGTGGGGGGGGACCCTCCAc                                                                      >  2:866691/1‑151 (MQ=255)
                                                                 cGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACg                                                                     <  2:463912/151‑1 (MQ=255)
                                                                  gTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                                                    <  1:383968/151‑1 (MQ=255)
                                                                       gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGTTTTCCCAAAAAACGGCGGCGGGGAAGCGTCCTCGTATCGGGGAACTGATTGGCCCGAGCCAGGGCGGGATGTGGGTCGGCCCCCTCCCCCGGCTg                                                                >  2:17673/1‑150 (MQ=255)
                                                                       gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCAACAAAAAGGCGGCGGGGAGGCGGCCTCGTATCGGGGGACTGCTGGGGCCCAGCCACGGGGGGGTGGGGGGGGGGACCCTCCACcgg                                                                   >  1:433349/1‑147 (MQ=255)
                                                                                 tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCCAAAAAACGGCGGGGGGGGTGCGTCATCGTATCGGGGAACTGATGGGCACGAGCCAGGGCGGGGTGTGGGGCGGGGCCCTCCACGGGCCGGCGCCCCGttt                                                     >  2:709121/1‑151 (MQ=255)
                                                                                  gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGttt                                                     <  2:50923/150‑1 (MQ=255)
                                                                                  gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                                                    <  1:821268/151‑1 (MQ=255)
                                                                                      gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                                                <  1:197087/151‑1 (MQ=255)
                                                                                          aaaaCTGCTCGCCATACTCGATTATGGCGGTGACAATTTTCAAAAAAACCGAGGCGGGAATGAGTCATCGGAACGGGGGACTGATGGGGGCGAGCCAGGGGGGGGTGTGGGGGGGGCCCCCCCACGGGCCGGGGCCCCGGTtgct                                                  >  2:82335/1‑145 (MQ=255)
                                                                                           aaaCTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCCACAAAAACGCGGCGGGGGTGCGTCATCCTATCGGGCCACTGATGGGGGCGAGCCCGGGGGGGATGTGGGGGGGCCCCCTCCACCGGCggggggccccctt                                                     >  1:235741/1‑131 (MQ=255)
                                                                                            aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTgcgc                                          <  1:715761/151‑1 (MQ=255)
                                                                                            aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAAACGCGGCGGGGATGCGGCATCCTATCGGGCAACTGAAGGGCACGAGGCACGGCGGGGTGGGGGGGGGGGCCCTCCCCCGGCCGGCGCACCGCTTGCttcg                                               >  1:231978/1‑146 (MQ=255)
                                                                                                 cTCGCCATACTCGATTATGGCGGTGAAGTTTACCAACAAAAACGCGGCGGGGGTGTGTCGTCGGTTCGGGCAAATGATGGGGCACGGCCAGGGGGGGGTGTGGGGGGGCCCCCCCCCCCGGCGGGGGCCCCGTTTTCTGCGGGCGCGccc                                      >  1:673333/1‑149 (MQ=255)
                                                                                                   cGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCTCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                    <  1:322919/150‑1 (MQ=255)
                                                                                                    gCCATACTCGATTATGGCGGTGACGTTTTCCCACAAAACGGCGGCGGGGATGCGTCATCGTATCGGGGAACAGATGGGGACGAGCCAGGGGGGGGTGTGGGGGGGCACCCTCCCCCGGCTGGCGCCCCCGTTTCTGCGTGGGcccc                                       >  2:830310/1‑146 (MQ=255)
                                                                                                       aTACTCGATTATGGCGGTGACGGTTTCCAACAAAACCGCGGCGGGGGTGCGTCATCCTATCGGGGAAATGATGGGGCCGAGCCCGGGGGGGGTGGGGGGGGGGCCCCTCCCCCGGCCGGGGCACCCGTTTCTGCGTGCGCCCCctattg                                 >  1:149832/1‑147 (MQ=255)
                                                                                                       aTACTCGATTATGGCGGGGACGGTTTTCAACAAAACGGCGGCGGAGATGCGTCATCCTGTCGGGCAAAAGATGGGGACGCACCAGGGGGGGATGTGGGGCGGCACCCTCCCCCGGGGGGCGGCCCGCTTGCTGCGTGCGCAccccatag                                 >  1:740745/1‑143 (MQ=255)
                                                                                                         aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                              <  1:843125/150‑1 (MQ=255)
                                                                                                         aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                              <  1:720665/150‑1 (MQ=255)
                                                                                                            cGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGcc                            <  2:833406/149‑1 (MQ=255)
                                                                                                            cGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGACAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCa                           <  1:263399/150‑1 (MQ=255)
                                                                                                                   ggCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCGCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGcc                    <  1:302735/150‑1 (MQ=255)
                                                                                                                          aCGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAtt            <  1:781022/151‑1 (MQ=255)
                                                                                                                           cGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAttt           <  1:642244/151‑1 (MQ=255)
                                                                                                                                   aacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc   >  1:697434/1‑151 (MQ=255)
                                                                                                                                   aacaaAGCGGCGGCGGAGATGCGTCATCATATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc   <  1:759392/151‑1 (MQ=255)
                                                                                                                                    acaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCt  >  1:632078/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCT  >  NC_000913/3998041‑3998323

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: