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

TTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206310‑3206582
                                                                                                                                      |                                                                                                                                          
ttCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaa                                                                                                                                       <  1:321122/140‑1 (MQ=255)
     tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTg                                                                                                                                  >  1:89016/1‑140 (MQ=255)
            gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATAGAGTCGc                                                                                                                           >  2:610248/1‑140 (MQ=255)
             cgcAATATTCATTCTTATTTTCGCGTTGATTCACTAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAACTGGAGTTGcc                                                                                                                          >  1:528260/1‑140 (MQ=255)
             cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGcc                                                                                                                          >  1:997298/1‑140 (MQ=255)
                        ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATGTAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                               <  1:259688/140‑1 (MQ=255)
                        ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                               <  2:89016/140‑1 (MQ=255)
                            tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAAAAGAGGTTGCATGATGAGGGaaaaa                                                                                                           >  2:646647/1‑138 (MQ=255)
                              ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATTGGGTTGCCATGATGAGGGGAAGtaa                                                                                                         >  2:105801/1‑140 (MQ=255)
                              ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAAAGGAGGTGCCATGATGAGGgaaagaaa                                                                                                         >  2:190283/1‑140 (MQ=255)
                               tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                                        <  1:498322/140‑1 (MQ=255)
                                     gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTCGCCATGATGAGCGAAAGTAATAAGCaa                                                                                                  <  1:610248/140‑1 (MQ=255)
                                      acagttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAAc                                                                                                 <  1:881180/136‑1 (MQ=255)
                                         aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                              <  2:281318/140‑1 (MQ=255)
                                                aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                                       <  1:105801/140‑1 (MQ=255)
                                                aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                                       <  1:190283/140‑1 (MQ=255)
                                                   tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATaa                                                                                    <  1:678636/140‑1 (MQ=255)
                                                     ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCAGACAAATAAAAAAAACGGCGATGTCATGAGGGGAGAAAGAAAAAAACAAAAAGCAGAGaaagttt                                                                                  >  2:881180/1‑135 (MQ=255)
                                                             atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAAAGGAGGTGGCCTCATGAGGGAAAGAAATAAGGAAAAGGCAGGGAATAAGATGATagag                                                                          >  2:766625/1‑140 (MQ=255)
                                                              tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                                         <  2:528260/140‑1 (MQ=255)
                                                               atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTTAAAAAAAAAAAAAAATGGGGGTTGCTTGGTGGGGGAAAGGAATAAGAAAAAGGGAGGGAATAAATTGACAGAGAt                                                                        >  2:592783/1‑140 (MQ=255)
                                                               atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGACAAAAAAAAAAAAAGGGGTTGGCATGATGAGCGGAAGTAATAAAAAACAAGCAGGGAAAAAATTGGCAgagaa                                                                        >  2:549941/1‑139 (MQ=255)
                                                               atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTTGGCATGATGAGGGAAAGAAAAAAACAAAAAGCAGTGAAAAAATTGACAgagaa                                                                        >  2:847045/1‑139 (MQ=255)
                                                               atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAAAAAAAAAAAAAAAGGAGGTGGCATGATGAGGGGAAGAAATAAAAAAAAAGCAGGGAAAAAATTGACagaaaa                                                                        >  2:170538/1‑137 (MQ=255)
                                                                gtCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAtt                                                                       <  1:608171/139‑1 (MQ=255)
                                                                     ggTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGGGAAAGAAATAAACAACAAGCAGAGAAAAAATTGAAAGAGATTTTCGc                                                                  >  2:511263/1‑140 (MQ=255)
                                                                            aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAAAGGGGGTGGCCTGATGAGCGAAAATAAAAAACAAAAAGCAGGGAAAAAAATGACAAAGAATGTTGCTAAAttt                                                           >  2:608171/1‑140 (MQ=255)
                                                                                   taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAAATGGAGTTGGCCTGGTGAGGGAAAGAAAAAAACAAAAAGCAGTGAAAAAATTGACAGAGAATGTTGCTAACTTTACCGccc                                                    >  2:73332/1‑139 (MQ=255)
                                                                                   aaaacCCCCATTTTCAGTGGTATTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCa                                                    <  1:592783/139‑1 (MQ=255)
                                                                                         cacaTTTTTAGTGGTAGTTTGTGGCGGGGGATTTTTTCAGAAAAATACAAAAACTAGAGATTGCATGATGAGCGGAAAGAATAAAAAAAAAGCCATGAATAAGATGAAAGAGAATGGTGCTAACTTTACCGCCATGAttt                                              >  2:146289/1‑140 (MQ=255)
                                                                                                  ggTGGTAGTTTGGGGGGGTGATTTTTTCAAGAAAAATACAAAAACTGGAGTTGCCATGAGGAGCAAAAGTAATAACAAACAGGCAGTGAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                                     <  2:997298/139‑1 (MQ=255)
                                                                                                         gTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAACCTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtga                              <  1:830239/140‑1 (MQ=255)
                                                                                                          tttGTGGCGGTGAATTTTTCAAAACAAAAACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                             <  2:778446/140‑1 (MQ=255)
                                                                                                                     gAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa                  <  1:821901/140‑1 (MQ=255)
                                                                                                                               aGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAAAAAACAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa        <  1:549941/140‑1 (MQ=255)
                                                                                                                               aGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa        <  1:146289/140‑1 (MQ=255)
                                                                                                                               aGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaaa        <  1:315836/140‑1 (MQ=255)
                                                                                                                                     aTACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTa  <  1:766625/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
TTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206310‑3206582

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