Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,6940CT46.5% 34.5 / 26.4 43F191F (TTC→TTT) yafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (6/17);  new base T (20/0);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 2.42e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCA  >  NC_000913/231554‑231839
                                                                                                                                            |                                                                                                                                                  
ccTGCCACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTg                                                                                                                                           <  2:498979/149‑1 (MQ=255)
ccTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGa                                                                                                                                         <  2:400171/151‑1 (MQ=255)
  tGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTCAA‑TTCCCgg                                                                                                                                          >  2:870893/1‑148 (MQ=255)
   gACCCCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGCTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGc                                                                                                                                       <  2:124336/150‑1 (MQ=255)
   gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAg                                                                                                                                        <  2:1106508/149‑1 (MQ=255)
      acCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTTTA‑TGCCTggg                                                                                                                                         >  2:834379/1‑144 (MQ=255)
          gCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAA‑TTCCCTGAGCCg                                                                                                                                     >  1:23805/1‑145 (MQ=255)
               ttGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAA‑TGCCTGGGGGCGTAGAAAGATg                                                                                                                           >  1:1048670/1‑150 (MQ=255)
                tGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGa                                                                                                                          <  2:133404/150‑1 (MQ=255)
                   tggtggTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAAcg                                                                                                                       <  1:786212/150‑1 (MQ=255)
                      tggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAA‑TGCCTGGAGCCGGAGAAGGATGAACGCGtt                                                                                                                   >  1:699279/1‑151 (MQ=255)
                        gtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATTACCTGG‑GCCGGAGAAGGGTGAAAGCGGtat                                                                                                                 >  2:454218/1‑151 (MQ=255)
                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGATTTAA‑TTCCTTGGGGCGGAGAAGGGTGAAAGCGGTTTATCg                                                                                                             >  1:811448/1‑149 (MQ=255)
                                 cacGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAA‑TGCCTGGGGCCGTAGAAGGGTGAAAGCGTTATATTGCctt                                                                                                         >  1:522413/1‑150 (MQ=255)
                                     ccGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCTA‑ATCCTGGGGGCCTAGAAGGGTGGAAGGCTTATATTGCTTTtcg                                                                                                      >  2:974571/1‑149 (MQ=255)
                                     ccGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGcgcg                                                                                                    <  2:699159/151‑1 (MQ=255)
                                       gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCggg                                                                                                  <  1:1022009/151‑1 (MQ=255)
                                           aCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TGCCTGGAGCCGCAGCAGGAGGAACGCGGTATATTGCTTTGCGgggg                                                                                                  >  2:90286/1‑147 (MQ=255)
                                               cAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTc                                                                                           <  2:826217/150‑1 (MQ=255)
                                                   cTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TGCCTGGGGCCGGAGAAGGGTGAAAGCGGTTTATCGGTTTGCGCGGGgaa                                                                                               >  2:151491/1‑142 (MQ=255)
                                                     gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TGCCTGGAGCCGTTGAAGGATGAAAGCGGTTTATCGCTTTTCGCGGGGAAAGTCGCTgcg                                                                                     >  1:853806/1‑150 (MQ=255)
                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAA‑TGCCTGGGGCCGGAGAAGGGTGAAAGCGGTTTATTGCTTTTTGCGGGgaa                                                                                               >  2:1087617/1‑138 (MQ=255)
                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCa                                                                                  <  2:267567/151‑1 (MQ=255)
                                                            gATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTCA‑AGCCTTGGGACCTAGAAGGGTGAACGCGGTATTTCGCATTTCCCGGGGCAAGTCGCTGCGCCAgggg                                                                              >  2:244825/1‑148 (MQ=255)
                                                               gTCTTTAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtt                                                                          <  2:324637/151‑1 (MQ=255)
                                                                 cTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTCAA‑TACCCGGAGGCCGAGAAGGATGAACGCGGTTTATTGCTTTTCGCGGGGAATGTCGCTGCGCCAGGTGGGgttt                                                                        >  2:544242/1‑151 (MQ=255)
                                                                   ataGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGttt                                                                         <  1:255991/148‑1 (MQ=255)
                                                                        aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAA‑TGCCTGGAGCCGGAGAAGGGTGAACGCGGTATATTGCTTTGCGCGGGGAAAGTCGCTGCGCCAGGTGCGGTTTTCCgatg                                                                 >  2:357496/1‑151 (MQ=255)
                                                                          gCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCgatgat                                                               <  1:229460/151‑1 (MQ=255)
                                                                             gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TGCCCGGGGCCGGAGAAAGGGGAACACGGTATATCTCTTTTTGCGGGGAATGTCGAGGGGCCAgggggggttt                                                                        >  2:637518/1‑131 (MQ=255)
                                                                              ttACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TTCCTGGAGGCGGAGAAGGGAGAACGCGTTATATCGCTTTTCGCGGGGAAAGTCGCTGCGCCAGGGGCGTTTTTCCCATGATCa                                                             >  1:1037575/1‑149 (MQ=255)
                                                                               tACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTCAA‑TTCCTGGGGGCGCAGGAAGATGAAAACGTTATATACCTTTGTGCGGGGAAAGTTGCCGCGGCAGGAGCGGTTTTCCCATGATCAgc                                                           >  2:810237/1‑150 (MQ=255)
                                                                                aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAgcgc                                                         <  2:779554/151‑1 (MQ=255)
                                                                                    ccTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTTAA‑AGCCTGGGGGCGGAGAAGGATGGACGCGGTATTTTGCTTTTCGCGGGGAATGTTGCCGCGCCAGGGGCGGTTTTCCgat                                                                  >  2:327337/1‑138 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAA‑TGCCTGGGGCCGGAGAAGGGTGAAAGCGGTTTATTGCTTTGCGCGGGGAAAGTTGCTGCGCCCGGGGGGGTTTTTCGATGATCAGCGCCGGCggg                                                  >  1:839261/1‑151 (MQ=255)
                                                                                        ttGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAA‑TTCCTGGGGCCGGAGAAAGATGAAAGCGGTATATCGCTTTTCGCGGGGAAAGTCGCTGCGCCCGGGGGGGTTTTACGATGATCAGCGCCGGCggg                                                  >  1:1093028/1‑150 (MQ=255)
                                                                                                 aGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTCAA‑AGCCTGGCGCCCTAGAAGGATGAAAGCGTTATATTGCCTTTCGCGGGGAAAGTCGCTGCGCCAGGGGGGGTTTTCCCATTATTCGCGCCGTCCGGGGGTTgg                                           >  2:94650/1‑148 (MQ=255)
                                                                                                       cTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAA‑TGCCTGGAGGCGTAGAAGGATGAAAGCGGTATATCGCTTTTCGCGGGGAATGTCGCTGCGCCAGGGGCGTTTTTCCGATGATCAGCGCCGTCGGGCGGTTGGTCGCCCGc                                   >  1:65303/1‑150 (MQ=255)
                                                                                                              caGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGAtttt                           <  2:495522/151‑1 (MQ=255)
                                                                                                                    cccGGTCATTATGGCGGGAGATTTTTA‑TTCCTGGAGCCGCAGGAGGGTGGAAGCGTTATATCGCCTTGCGCGGGGAATGTTGCGGCGCCAGGGGCGGTTTTCCGAAGGTCAACGCCGCCGGGGGGTTGGTGGCCCCCCCGGTttttttttt                     >  2:250725/1‑150 (MQ=255)
                                                                                                                     ccGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCt                    <  1:333492/151‑1 (MQ=255)
                                                                                                                              aTGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTTCGCGGGAAATGTCGCTGCGCCAGGGGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGGTTTTGTTTTCTACCGTgg             >  1:734649/1‑149 (MQ=255)
                                                                                                                                       gaTTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCa  <  2:901629/151‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAA‑TGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCA  >  NC_000913/231554‑231839

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