Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT44.1% 28.7 / 19.4 34intergenic (‑179/‑28)ttdR/ttdAtranscriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (4/15);  new base T (15/0);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 2.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206295‑3206581
                                                                                                                                             |                                                                                                                                                   
cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATa                                                                                                                                           <  2:672046/151‑1 (MQ=255)
           gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                                <  1:85703/151‑1 (MQ=255)
             cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATAAAAAAAATgg                                                                                                                                 >  1:415326/1‑148 (MQ=255)
               ttCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGTAGtt                                                                                                                            <  2:540686/150‑1 (MQ=255)
                  ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                        <  2:155724/151‑1 (MQ=255)
                      ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAaaaaaa                                                                                                                                     >  2:661845/1‑135 (MQ=255)
                          ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGTTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGTCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCAtgatga                                                                                                                  <  2:299010/149‑1 (MQ=255)
                          ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGc                                                                                                                <  2:127040/151‑1 (MQ=255)
                           gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                               <  1:408727/151‑1 (MQ=255)
                                      attcTTCTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGCTTCCACACATCTCCAACCACCTTTTTAGTGGTAGTTTGTGGCGGTGTAATTTTTCCGTCCAAAACAAAAACTCTAGATGCCCCGATTAGTGCAAGtaata                                                                                                      >  1:568488/1‑149 (MQ=255)
                                            aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                              <  1:684959/151‑1 (MQ=255)
                                               ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGGGTTGGCATGATGAGCGgaaagaaaaaacaa                                                                                                 >  2:415711/1‑134 (MQ=255)
                                                 cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAGACAAATAAAAAAACAGGGGGTGCCATGGAGAGGGAAAGtaataaaaaa                                                                                                 >  1:490141/1‑139 (MQ=255)
                                                  gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTg                                                                                        <  2:391475/151‑1 (MQ=255)
                                                       gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAACTGGAGTTGCCATGATGAGCGAAAGGAAAAAACAAAAGGCAGGGaa                                                                                      >  1:700436/1‑148 (MQ=255)
                                                           cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGTATTTTTTTAGACAAATAAAAAAAATGGAGTTGCCATGATGGGCGAAAGGAATAAACAAAAAGGAGGGAATaa                                                                                   >  1:87715/1‑147 (MQ=255)
                                                             cGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAAAGGAGTTGTCATGATGAGCGGAAAGAATAAAAAAAAAGCAGTGAATaa                                                                                   >  1:449086/1‑145 (MQ=255)
                                                               aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACag                                                                           <  2:541531/151‑1 (MQ=255)
                                                                          aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCCCACAAATACAAAAAAATGAGATTGCCTGATGAGCGGGAGTAATAAAGAAAAAGCAGGGAATAAAtttggaaagagt                                                                      >  1:10472/1‑135 (MQ=255)
                                                                            atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCAGAAAAAAAAAAAAACTGGAGTTGGCAAGATGAGCGAAAGAAATAAAAAACAAGCCGAGAATAAAATGACagaagatgttg                                                                  >  2:344970/1‑139 (MQ=255)
                                                                             tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCACACAAAAAAAAAAAATGGAGGTGCCATTATGAGCGAAAAAAATAAAAAAAAAGCAGGGAATAAATTGACAGAGATTTTTGCtta                                                              >  1:21212/1‑148 (MQ=255)
                                                                                       ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCa                                                   <  2:651812/151‑1 (MQ=255)
                                                                                               atatAAAACCACATTTTTAGTGGTAGTTTGGGGCGGTGAATTTTTCCAAACAAATACAAAAACGGGAGTTGCCTTGATGACCAAAATTAATAAGCAACGGCCAGGGAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTc                                            <  2:695176/150‑1 (MQ=255)
                                                                                               ataaAAACCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATATAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                                           <  2:181727/151‑1 (MQ=255)
                                                                                                           aTTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAACTGGAGTTGCCAAGATGAGCGAAAGTAATAAACAAAAGGCAGTGAATAAGTTGACAGAGATTTTCGCTAACTTTACCGCCATGATTTCTACCCGAATGcc                                >  1:423969/1‑150 (MQ=255)
                                                                                                            tttttAGTGGTAGTTTGTGGCGGTGACTTTTTTCAGACAAAAACAAAAAATGGGGTTGGCCCTATGTGCGCAAGGAAAAAGCAAAAAGCGGTGGAAAAATTGACAGAGATTTTCGCTAACTTTATCGCCCTCATTTCTAcccca                                     >  1:294112/1‑142 (MQ=255)
                                                                                                                tAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                           >  1:517445/1‑150 (MQ=255)
                                                                                                                        gTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGGTGCCCTGATGAGCGAAAGTAAAAAGCAAAAAGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTACCCGAATGACTGATGACGtggggg                  >  2:597485/1‑151 (MQ=255)
                                                                                                                        gTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGtggtgg                  <  1:671702/151‑1 (MQ=255)
                                                                                                                          ttGTGGCGGTGAAATTTTTCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAAAAAAAAACAAGGAGTGAATAAGTTGACAGAGATTTTCGCTAAATTTACCGCCATGATTTTTACCCGAAAGCCTGATGACgtgggg                   >  2:634495/1‑146 (MQ=255)
                                                                                                                                      aaTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc      >  1:584506/1‑150 (MQ=255)
                                                                                                                                      aaTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAATTAATAAGCAACAGGCATTAAATAATTTCAAAGAGTTTTTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACa     <  2:691266/151‑1 (MQ=255)
                                                                                                                                        tttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGc   >  2:665870/1‑151 (MQ=255)
                                                                                                                                         ttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt  >  1:248266/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206295‑3206581

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