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

CGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCT  >  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)
                         ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑g                                                                                                             <  1:487594/151‑1 (MQ=255)
                            gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGc                                                                                                          <  1:578758/151‑1 (MQ=255)
                            gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGc                                                                                                          <  1:473206/151‑1 (MQ=255)
                              gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCAc                                                                                                        <  1:581883/151‑1 (MQ=255)
                               cgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACg                                                                                                       <  1:929060/151‑1 (MQ=255)
                                  gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGTTTACCCACAAAACGGCGGCGGGGGGGCGTCATCGTATCGGGCGAATGAT‑ggg                                                                                                           >  1:390081/1‑144 (MQ=255)
                                      aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAgg                                                                                                <  1:729691/151‑1 (MQ=255)
                                        tGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGGTTTCCAACAAAACGGGGGCGGGGAAGGGGCACCGTATACGGCAACCGGA‑GGGCACGAGcc                                                                                                   >  2:264679/1‑146 (MQ=255)
                                             aGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACCTTTTCTCAAAAAACGGCGGCGGGGATGCGTCCTCGGCTCGGGGAACAGAT‑GGGCCCGAGcccaggggg                                                                                            >  2:940861/1‑141 (MQ=255)
                                                   gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACGGCGGGGGGGGGGGGTCATCGTATCCGGGGACAGAT‑GGGGACCGGCCAGGGGGGGGTGTggg                                                                                    >  2:118728/1‑150 (MQ=255)
                                                     gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCg                                                                                 <  1:852609/151‑1 (MQ=255)
                                                     gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTTCAACAAAACGGCGGCGGGGGAGCGGCCTCGTATTGGGGGAATGAT‑GGGGACGAGCCAGGGGGGGGTGTggg                                                                                    >  2:397336/1‑148 (MQ=255)
                                                            tGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCAACAAAACGGCGGCGGGGAGGGGTCATCGTATCTGGGAACTGAT‑GGGGACGAGCCAGGGGGGGGTGTGGGGCGGGAccct                                                                          >  2:856894/1‑149 (MQ=255)
                                                             gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTc                                                                         <  2:65878/151‑1 (MQ=255)
                                                             gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACACGGGCGGGGGTGCGGCCTCATCTCTGGCCACCGAT‑GGGCGCGAGCCCGGGGGGGTTGTGGGGCggggccct                                                                          >  1:369993/1‑144 (MQ=255)
                                                                tCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCAAAAAAACGGGGGCGGGGATGCGTCATCGGATTGGGGGACCGGT‑GGGGACGAGCCCGGGGGGGGTGTGGGGGGGGACCCTCCAc                                                                      >  2:866691/1‑151 (MQ=255)
                                                                 cGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACg                                                                     <  2:463912/151‑1 (MQ=255)
                                                                  gTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                                                    <  1:383968/151‑1 (MQ=255)
                                                                       gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGGCGTTTTCCCAAAAAACGGCGGCGGGGAAGCGTCCTCGTATCGGGGAACTGAT‑TGGCCCGAGCCAGGGCGGGATGTGGGTCGGCCCCCTCCCCCGGCTg                                                                >  2:17673/1‑150 (MQ=255)
                                                                       gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTTCAACAAAAAGGCGGCGGGGAGGCGGCCTCGTATCGGGGGACTGCT‑GGGGCCCAGCCACGGGGGGGTGGGGGGGGGGACCCTCCACcgg                                                                   >  1:433349/1‑147 (MQ=255)
                                                                                 tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCCAAAAAACGGCGGGGGGGGTGCGTCATCGTATCGGGGAACTGAT‑GGGCACGAGCCAGGGCGGGGTGTGGGGCGGGGCCCTCCACGGGCCGGCGCCCCGttt                                                     >  2:709121/1‑151 (MQ=255)
                                                                                  gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGttt                                                     <  2:50923/150‑1 (MQ=255)
                                                                                  gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                                                    <  1:821268/151‑1 (MQ=255)
                                                                                      gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                                                <  1:197087/151‑1 (MQ=255)
                                                                                          aaaaCTGCTCGCCATACTCGATTATGGCGGTGACAATTTTCAAAAAAACCGAGGCGGGAATGAGTCATCGGAACGGGGGACTGAT‑GGGGGCGAGCCAGGGGGGGGTGTGGGGGGGGCCCCCCCACGGGCCGGGGCCCCGGTtgct                                                  >  2:82335/1‑145 (MQ=255)
                                                                                           aaaCTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCCACAAAAACGCGGCGGGGGTGCGTCATCCTATCGGGCCACTGAT‑GGGGGCGAGCCCGGGGGGGATGTGGGGGGGCCCCCTCCACCGGCggggggccccctt                                                     >  1:235741/1‑131 (MQ=255)
                                                                                            aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTgcgc                                          <  1:715761/151‑1 (MQ=255)
                                                                                            aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAAACGCGGCGGGGATGCGGCATCCTATCGGGCAACTGAA‑GGGCACGAGGCACGGCGGGGTGGGGGGGGGGGCCCTCCCCCGGCCGGCGCACCGCTTGCttcg                                               >  1:231978/1‑146 (MQ=255)
                                                                                                 cTCGCCATACTCGATTATGGCGGTGAAGTTTACCAACAAAAACGCGGCGGGGGTGTGTCGTCGGTTCGGGCAAATGATGGGGCACG‑GCCAGGGGGGGGTGTGGGGGGGCCCCCCCCCCCGGCGGGGGCCCCGTTTTCTGCGGGCGCGccc                                      >  1:673333/1‑149 (MQ=255)
                                                                                                   cGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCTCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                    <  1:322919/150‑1 (MQ=255)
                                                                                                    gCCATACTCGATTATGGCGGTGACGTTTTCCCACAAAACGGCGGCGGGGATGCGTCATCGTATCGGGGAACAGAT‑GGGGACGAGCCAGGGGGGGGTGTGGGGGGGCACCCTCCCCCGGCTGGCGCCCCCGTTTCTGCGTGGGcccc                                       >  2:830310/1‑146 (MQ=255)
                                                                                                       aTACTCGATTATGGCGGTGACGGTTTCCAACAAAACCGCGGCGGGGGTGCGTCATCCTATCGGGGAAATGAT‑GGGGCCGAGCCCGGGGGGGGTGGGGGGGGGGCCCCTCCCCCGGCCGGGGCACCCGTTTCTGCGTGCGCCCCctattg                                 >  1:149832/1‑147 (MQ=255)
                                                                                                       aTACTCGATTATGGCGGGGACGGTTTTCAACAAAACGGCGGCGGAGATGCGTCATCCTGTCGGGCAAAAGAT‑GGGGACGCACCAGGGGGGGATGTGGGGCGGCACCCTCCCCCGGGGGGCGGCCCGCTTGCTGCGTGCGCAccccatag                                 >  1:740745/1‑143 (MQ=255)
                                                                                                         aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                              <  1:843125/150‑1 (MQ=255)
                                                                                                         aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                              <  1:720665/150‑1 (MQ=255)
                                                                                                            cGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGcc                            <  2:833406/149‑1 (MQ=255)
                                                                                                            cGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGACAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCa                           <  1:263399/150‑1 (MQ=255)
                                                                                                                   ggCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCGCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGcc                    <  1:302735/150‑1 (MQ=255)
                                                                                                                          aCGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAtt            <  1:781022/151‑1 (MQ=255)
                                                                                                                           cGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAttt           <  1:642244/151‑1 (MQ=255)
                                                                                                                                   aacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc   >  1:697434/1‑151 (MQ=255)
                                                                                                                                   aacaaAGCGGCGGCGGAGATGCGTCATCATATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc   <  1:759392/151‑1 (MQ=255)
                                                                                                                                    acaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCt  >  1:632078/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAT‑GGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCT  >  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: ‑