Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4440TA32.4% 55.3 / 10.7 37intergenic (‑188/‑19)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (5/20);  new base A (12/0);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 3.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGC  >  NC_000913/3206315‑3206580
                                                                                                                                 |                                                                                                                                        
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATg                                                                                                                                >  1:821877/1‑140 (MQ=255)
  ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGGAATTTTTCATACAaaaaaaaaaaaaggg                                                                                                                              >  1:185709/1‑135 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                             <  2:131771/140‑1 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                             <  1:574757/140‑1 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                             <  1:341692/140‑1 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAAACAAAAAAAAAAAATGGAGGTGc                                                                                                                         >  2:574757/1‑140 (MQ=255)
                 cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGTCATGATGAgg                                                                                                               >  2:295439/1‑139 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                             <  2:539091/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                             <  2:644956/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                             <  1:663910/140‑1 (MQ=255)
                       tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGGCATTATGAGGGAAAga                                                                                                         >  1:781801/1‑139 (MQ=255)
                       tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCACAAAAAAAAAAAAAATGGAGTTGGCCTGATGAGGGaaaaa                                                                                                         >  2:856190/1‑138 (MQ=255)
                       tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAAATTTTTCCGACAAAAAAAAAAAATGGAGGTGCCCTGATGAGCGAAAgg                                                                                                         >  1:569859/1‑139 (MQ=255)
                           ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTTCGACAAAAAAAAAAAATGGGGGTGGCATGATTAGAGGAAGtaaaa                                                                                                     >  2:663910/1‑138 (MQ=255)
                              gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGCGAAAGgaaaaaaa                                                                                                  >  1:635023/1‑138 (MQ=255)
                                 ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAAc                                                                                               <  1:856190/140‑1 (MQ=255)
                                    aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                            <  2:177405/140‑1 (MQ=255)
                                       cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAg                                                                                         <  1:873230/140‑1 (MQ=255)
                                                ttAATTACATACTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGGAAGGAATAAAGAAAAAGCAGTGAATaaat                                                                                >  1:639986/1‑138 (MQ=255)
                                                    ttACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAAAAGGAGGTGGCATGATGAGGGAAAATAATAAAAAAAAAGCAGGGAAAAAATTGaa                                                                            >  2:12582/1‑139 (MQ=255)
                                                     tACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACa                                                                           <  1:271023/140‑1 (MQ=255)
                                                        atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGGAAAGAAAAAACAACAAGCAGTGAATAAATTGACagag                                                                        >  2:775348/1‑140 (MQ=255)
                                                         tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                                       <  2:428582/140‑1 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAACAAAAAATGGAGGTGGCATGATGAGGGGAAGGAATAAAAAAAAAGCAGGGAATAAATTGACAgagaa                                                                      >  2:271023/1‑139 (MQ=255)
                                                                  tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAAAAAAAAAAAAAAAGGAGGTGGCATGATTAGGGAAAAAAAAAAACAACAAGCAGTGAAAAAGTTGACAGAGAATGTCGCTa                                                              >  2:812131/1‑140 (MQ=255)
                                                                              taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGAGGTGCCATGATGAGGGAAAAGAATAAACAACAAGCAGTGAAAAAAATGACAGAGATTGTCGCTAACTTTatccccc                                                  >  2:464784/1‑134 (MQ=255)
                                                                                        tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTAc                                        <  1:699514/140‑1 (MQ=255)
                                                                                              gTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                                  <  1:12582/140‑1 (MQ=255)
                                                                                                    ttttGTGGCGGTGAATTTTTCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtga                            <  2:635023/139‑1 (MQ=255)
                                                                                                    gTTTGTGGCGGTGAATTTTTCAAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCtga                            <  2:781801/140‑1 (MQ=255)
                                                                                                      ttGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                          <  1:198980/140‑1 (MQ=255)
                                                                                                          ggCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgt                      <  2:501048/140‑1 (MQ=255)
                                                                                                                gAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa                <  1:775348/140‑1 (MQ=255)
                                                                                                                aaaTTTTCCCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACGGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa                <  2:259392/139‑1 (MQ=255)
                                                                                                                 aaTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGAt               <  1:360535/140‑1 (MQ=255)
                                                                                                                     tttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc           >  2:306441/1‑140 (MQ=255)
                                                                                                                        caaGACAAATAAAAAAAGTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaa        <  1:437628/138‑1 (MQ=255)
                                                                                                                         aaaaaaaaTACAAAAACGGGAGTTGCCATGAAGAGCGAAAGTAATAACCAACAGGAAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa       <  2:569859/135‑1 (MQ=255)
                                                                                                                          aaaCAAATAAAAAAACTGGAGTTGCCATGAGGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACGGAGTTTGTCGTTAACTTTTCCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa      <  2:754159/138‑1 (MQ=255)
                                                                                                                          aGAAAAATACAAAAACGGGAGTTGCAATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa      <  2:475113/140‑1 (MQ=255)
                                                                                                                           gACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc     >  1:802866/1‑140 (MQ=255)
                                                                                                                              aaaTACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGc  >  2:360624/1‑140 (MQ=255)
                                                                                                                              aaaTACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGc  >  1:372300/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                        
TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGC  >  NC_000913/3206315‑3206580

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