Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176353,419,1480CT26.2% 95.7 / 15.1 42intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator TtdR/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (4/27);  new base T (11/0);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 3.19e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG  >  NC_017635/3419099‑3419196
                                                  |                                                
aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTc                                                  <  2:2054577/50‑1 (MQ=255)
  aaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCa                                                <  2:786572/50‑1 (MQ=255)
  aaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAAATTTTTcc                                                 >  2:1691426/1‑50 (MQ=38)
   aaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAg                                               <  1:719587/50‑1 (MQ=255)
    atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTccga                                               >  1:411343/1‑48 (MQ=38)
    atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAg                                               >  1:799176/1‑49 (MQ=255)
    atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAAATTTTTCAGa                                              >  2:1642578/1‑50 (MQ=38)
     tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAc                                             <  2:1510728/50‑1 (MQ=255)
      ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAc                                             >  1:330619/1‑49 (MQ=255)
       taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACa                                            >  1:1924491/1‑49 (MQ=38)
        aaaaCCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAGACaa                                           >  2:644533/1‑49 (MQ=38)
         aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACaaa                                          >  2:1918104/1‑49 (MQ=255)
          aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACaaaa                                         >  1:384629/1‑48 (MQ=255)
             cacaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaa                                     <  1:216750/50‑1 (MQ=255)
              acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCTGACAAAAACaa                                     >  1:2052985/1‑49 (MQ=38)
                 tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAAc                                 <  2:1701947/50‑1 (MQ=255)
                  ttttAGTGGTAGTTTGTGGCGGTGGATTTTTTCCGACAaaaaaaaaaa                                  >  2:1087219/1‑48 (MQ=25)
                    ttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTg                               <  2:1831646/49‑1 (MQ=255)
                     tAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTgg                              <  2:1374345/49‑1 (MQ=255)
                      aGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAACAAAAACTgggg                            >  1:1001872/1‑48 (MQ=255)
                      aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                             <  2:1995844/49‑1 (MQ=255)
                         ggTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAAAGGAGTTg                         >  1:932621/1‑50 (MQ=255)
                             gTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCat                     <  1:955076/50‑1 (MQ=255)
                              tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                    <  2:183092/50‑1 (MQ=255)
                              tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                    <  2:30102/50‑1 (MQ=255)
                              tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                    <  2:970388/50‑1 (MQ=255)
                                tgtgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                   <  1:983295/49‑1 (MQ=255)
                                tgtgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                   <  1:264208/49‑1 (MQ=255)
                                  tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgatg                 <  2:1728749/49‑1 (MQ=255)
                                  tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgatg                 <  1:1819028/49‑1 (MQ=255)
                                  tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgat                  >  2:1590504/1‑48 (MQ=255)
                                  tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCAtgatga                <  1:317642/50‑1 (MQ=255)
                                    gCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGc              <  1:2068931/50‑1 (MQ=255)
                                     cGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg             <  2:1663936/50‑1 (MQ=255)
                                     cGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg             <  1:5858/50‑1 (MQ=255)
                                       gTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa           <  1:1996796/50‑1 (MQ=255)
                                       gTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa           <  2:64251/50‑1 (MQ=255)
                                       gTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa            >  1:455938/1‑49 (MQ=255)
                                           aaTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAg         <  2:915574/49‑1 (MQ=255)
                                             tttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa      <  2:58026/50‑1 (MQ=255)
                                              ttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa      >  2:346892/1‑49 (MQ=255)
                                                ttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata    <  2:1420594/49‑1 (MQ=255)
                                                ttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat     >  2:1717937/1‑48 (MQ=255)
                                                 tCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa   <  2:299584/49‑1 (MQ=255)
                                                 tCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa   <  2:1920774/49‑1 (MQ=255)
                                                  ccAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAg  >  1:1721994/1‑49 (MQ=255)
                                                  ccAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAg  <  1:1160851/49‑1 (MQ=255)
                                                  |                                                
AAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG  >  NC_017635/3419099‑3419196

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