Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT32.3% 40.5 / 18.3 31intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/21);  new base T (10/0);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 2.25e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAA  >  NC_000913/3206298‑3206569
                                                                                                                                          |                                                                                                                                      
aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCa                                                                                                                                      <  1:1457138/140‑1 (MQ=255)
              cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAcaaaaacaaaaa                                                                                                                        >  1:124998/1‑140 (MQ=255)
                  gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTgg                                                                                                                    <  1:1113507/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACTTATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  1:357606/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  2:999090/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  2:539604/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  2:124998/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  1:81076/140‑1 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                  <  1:1304654/140‑1 (MQ=255)
                             atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACTTATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTAAATTTTTTCTGAAAAAAAAAAAAAATGAGTTgtccaga                                                                                                         >  2:357606/1‑134 (MQ=255)
                             atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAAGAGTTGTCatga                                                                                                         >  2:1304654/1‑140 (MQ=255)
                                 tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                     <  1:1112492/140‑1 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                  <  2:1455889/140‑1 (MQ=255)
                                        tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAACTGGAGTTGCCATGATGAGCGAAAga                                                                                              >  1:1263438/1‑139 (MQ=255)
                                                 gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAGGAAAAAGaaa                                                                                     >  1:562702/1‑140 (MQ=255)
                                                      ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                <  2:1333804/140‑1 (MQ=255)
                                                                 ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGGTGAGCGGAAGGAAAAAAAAACAAGCAGGGAATaaat                                                                     >  2:990329/1‑138 (MQ=255)
                                                                          tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                            <  1:1509029/140‑1 (MQ=255)
                                                                           atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAt                                                           <  1:990329/140‑1 (MQ=255)
                                                                           atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATAGAGGTGGCATGATGAGGGAAAGTAAAAAACAAAAAGCAGGGAATAAGATGACAaagaa                                                           >  2:839525/1‑139 (MQ=255)
                                                                               aCGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTc                                                       <  1:1311211/140‑1 (MQ=255)
                                                                                  gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTACACAAAAACAAAAACTGGAGTTGCCATGATGGGCGGAAAAAAAAAACAAAAAGCAGTGAATAAAATTAAAAAGATTTTCGcc                                                    >  2:1412477/1‑139 (MQ=255)
                                                                                         aaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                             <  1:817451/140‑1 (MQ=255)
                                                                                           gatatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTAcc                                           <  1:1279779/139‑1 (MQ=255)
                                                                                                       caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGGCCTGATGGGGGGAAATAATAAACAAAAAGCAGGGAATAAATTGACAGAGATTGTCGCTAAATTTTCCGCCCTGATTTcc                               >  2:1279779/1‑139 (MQ=255)
                                                                                                            ttAGTGGTAGTTTGTGGGGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAAAAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCg                          <  1:1412477/140‑1 (MQ=255)
                                                                                                              ggTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                        <  1:576475/139‑1 (MQ=255)
                                                                                                               gTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCAATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                       <  2:1263438/140‑1 (MQ=255)
                                                                                                                      tttGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGCGAAAGAAAAAAAAAAAAAGCAGGGAAAAAATTGACAGAGAATGTTGCTAACTTTACCCCCCTGATTTTTACCCCAATGCCtgat                >  2:1441637/1‑140 (MQ=255)
                                                                                                                        tgtgGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGGCATGATGAGCGGAAATAATAAAAAAAAAGCAGGGAAAAAATTGACAGAGATTGTCGCTAACTTTACCCCCCCGATTTTTTCCCCAAAGCCtgatga              >  2:440000/1‑140 (MQ=255)
                                                                                                                               gTGAATTTTTTCAGACAAAAACAAAAAAAGGAGGTGGCATGATGAGCGGAAGGAAAAAGCAACAAGCAGGGAAAAAATTGACAGAGATTGTCGCTAAATTTACCGCCATGATTTTTACCCCAATGCCTGATGACgggggg       >  2:1478101/1‑138 (MQ=255)
                                                                                                                               gTGAAATTTTTCACACAAATACAAAAAAAGGGGGTGGCATCAAGAGGGAAAGGAAAAAACAACAAGCAGTGAATAAGTTGACAGAGAATGTCGCTAAATTTAACCCCATGATTTTTACCCCAATACCCGATGAAggggtg       >  2:805066/1‑140 (MQ=255)
                                                                                                                                  gAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa     <  1:1478101/140‑1 (MQ=255)
                                                                                                                                  aaaTTTTTCCAGAAAAATACAAAAACTGGATTTCCCATGATGAGCAAAAGTAAAAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa     <  2:562702/139‑1 (MQ=255)
                                                                                                                                     tttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCCTGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa  <  1:440000/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAA  >  NC_000913/3206298‑3206569

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: