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

AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAA  >  NC_000913/3206298‑3206574
                                                                                                                                         |                                                                                                                                           
aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCa                                                                                                                                           >  1:523895/1‑140 (MQ=255)
    cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACa                                                                                                                                       <  2:978583/140‑1 (MQ=255)
    cATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACa                                                                                                                                       >  1:764078/1‑140 (MQ=255)
         ccATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                  <  1:566995/140‑1 (MQ=255)
         ccATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                  <  1:603719/140‑1 (MQ=255)
            ttCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaa                                                                                                                               >  2:924370/1‑140 (MQ=255)
                    tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                       <  2:939516/140‑1 (MQ=255)
                            aaTATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                                                                                                               <  2:841438/140‑1 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                       <  1:552679/140‑1 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                       <  1:924370/140‑1 (MQ=255)
                                    ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                       <  2:178691/140‑1 (MQ=255)
                                     tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                                                                                      <  2:712010/140‑1 (MQ=255)
                                          ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAAAAAAAAAAATGGAGGTGGCATGATTAGGGaaaaaaa                                                                                                 >  2:85830/1‑136 (MQ=255)
                                               gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCAATGATGAGCGAAAGTAATAAGc                                                                                            <  1:863900/140‑1 (MQ=255)
                                                 gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCaa                                                                                          <  2:411582/140‑1 (MQ=255)
                                                     aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                      <  2:882018/140‑1 (MQ=255)
                                                     aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                      <  1:648345/140‑1 (MQ=255)
                                                      ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                     <  1:890046/140‑1 (MQ=255)
                                                       tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                    <  1:85830/140‑1 (MQ=255)
                                                       tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                    <  1:571055/140‑1 (MQ=255)
                                                               tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGGTGCCATGGTGGGCGGAAATAAAAAAAAAAAAGCAGGGAaaaa                                                                            >  2:852947/1‑140 (MQ=255)
                                                                         atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAATGGGGGTGGCATGGTGAGGGGAAGGAATAAAAAAAAAGCAGTGAATAAATTGGAagag                                                                  >  2:926612/1‑140 (MQ=255)
                                                                         atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAAAAAAAAAAAAGGAGGTGCCATGATGAGGGGAAAAAATAAACAAAAAGCAAGGAAAAAAATGACagaa                                                                  >  2:317158/1‑139 (MQ=255)
                                                                         atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATGATGAGCGAAAATAAAAAACAAAAAGCAGTGAATAAAGTGACagaa                                                                  >  2:719709/1‑139 (MQ=255)
                                                                           atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGGAAGAAATAAACAAAAAGCAGGGGATAAGTTGACAgagaa                                                                >  2:275672/1‑139 (MQ=255)
                                                                                                   aCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCAAAACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGAAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                        <  2:665011/140‑1 (MQ=255)
                                                                                                            ttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGGCATGATGAGGGAAAAGAATAAGCAACAGGCAGTGAATAAATTGACAGAGATTGTCGCTAACTTTACCCCCCTGAttttttaccc                               >  2:824635/1‑133 (MQ=255)
                                                                                                              aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                             <  1:824635/140‑1 (MQ=255)
                                                                                                               gTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                            <  1:926612/140‑1 (MQ=255)
                                                                                                                      tttGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                     <  1:317158/140‑1 (MQ=255)
                                                                                                                       ttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                    <  2:665576/140‑1 (MQ=255)
                                                                                                                       ttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                    <  1:939060/140‑1 (MQ=255)
                                                                                                                          tgGCGGTGGATTTTTTCAGACAAAAACAAAAAATGGGGTTGTCCTGATGGGCGAAAATAATAAACAAAAAGCAGAGAAAAAGTTGACAGAGATTGTCGCTAACTTTACCCCCCTGATTTTTTCCCCAATGCCTGATGACg                 >  2:147497/1‑140 (MQ=255)
                                                                                                                            gCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg               <  1:39336/140‑1 (MQ=255)
                                                                                                                               gggAATTTTTTCAGACAAAAACAAAAAATGGAGGTGGCATGATGAGCGGAAGTAATAAACAAAAGGCAGTGAATAAGTTGACAGAGAATGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtggtg            >  2:22526/3‑140 (MQ=255)
                                                                                                                               gTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTTTCCTTATTAGCGAAAATAATAAACAAAAAGAAGTGAATAAATTGACAGAGATTGTCGCTAAATTTTCCGCCATGATTTTTTCCCGAATGCCTGATGACgtgggg            >  2:833310/1‑138 (MQ=255)
                                                                                                                                      tttCCAGAAAAATACAAAAACTGGAGTTGCGATGATGAGCGAAAGAAAAAACCAACGGGCAGTAAATAAGTGGACAGAGATTGCCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc     <  2:250314/140‑1 (MQ=255)
                                                                                                                                         ccAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaa  <  2:797418/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
AAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAA  >  NC_000913/3206298‑3206574

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: