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

TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAA  >  NC_000913/3206315‑3206575
                                                                                                                                  |                                                                                                                                    
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCTTCCAaaaaaaaaaaacg                                                                                                                            >  1:293829/1‑137 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                        <  1:440568/140‑1 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAATGGAGGTga                                                                                                                    >  1:132701/1‑139 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGGTgg                                                                                                                    >  1:687112/1‑139 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTCCACACAAATAAAAAAAATGGAGGTGc                                                                                                                    >  1:304011/1‑140 (MQ=255)
        cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGcc                                                                                                                   >  2:409401/1‑140 (MQ=255)
            atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGTTGGCatta                                                                                                               >  1:204697/1‑138 (MQ=255)
                tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                           <  1:248944/140‑1 (MQ=255)
                tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                           <  1:409401/140‑1 (MQ=255)
                tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGACAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                           <  1:243642/140‑1 (MQ=255)
                  attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCCAACAAAAAAAAAAAAATGGGGTTGCCTGATGAggg                                                                                                         >  2:248944/1‑138 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                        <  1:307960/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                        <  1:400853/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                        <  1:131010/140‑1 (MQ=255)
                       tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAACAAAAAAAAAAAAATGAGTTGCCCTGATTAGCGaaaaa                                                                                                    >  1:370679/1‑138 (MQ=255)
                       tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATTGAGGTGCCATGATGAGCGaaaaa                                                                                                    >  2:400853/1‑138 (MQ=255)
                        aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATTGAGGTGCCATGATGAGGGaaaaga                                                                                                   >  2:307960/1‑137 (MQ=255)
                         ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGAGTTGCCCTGATGAGCGAAAGtaa                                                                                                  >  1:84603/1‑140 (MQ=255)
                            tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                               <  2:402735/140‑1 (MQ=255)
                                              tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATaa                                                                             <  1:429055/140‑1 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAAGAGGTTGCCAGAAGAGAGGAAGGAAAAAGCAACACACAAGGAATAAAATGAAAGAGAt                                                                 >  2:471885/1‑140 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTTGGCCTGATGAGGGGAAATAATAAGCAAAAAGCAGTGAATAAATTGACAgagaa                                                                 >  2:326186/1‑139 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCAAACAAAAACAAAAAAATGAGGTGCCATTTTGAGGGAAAGAAAAAAAAAAAAAGCAGTGAATAAGTTGACAgagaa                                                                 >  2:728556/1‑139 (MQ=255)
                                                            tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGGCATGATGAGGGAAAAGAAAAAAAAAAAAGCAGTGAAAAAAATTACAgagaatt                                                               >  2:114771/1‑137 (MQ=255)
                                                            tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGGCATGATGAGGGAAAATAAAAAACAAAAAGCAGTGAAAAAAATGACAGAGAttt                                                               >  2:337016/1‑139 (MQ=255)
                                                            tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTg                                                               <  1:101902/140‑1 (MQ=255)
                                                                  tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAATGGAGGTGGCATGATGAGGGAAAAGAAAAAGAAAAAGGCAGGGAAAAAATTGAAAGAGATTTTTGCTa                                                         >  2:218482/1‑140 (MQ=255)
                                                                    gCAAAAAATATAAAACCACATTTTTAGTAGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGGCCAGATGAGCGAAAGGAATAAACAACAAGCAGGGAATAAATTGAAAGAGAATGTTGCTaaa                                                       >  2:709944/1‑139 (MQ=255)
                                                                                      caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATAGAGTTGGCCTGATGGGGGAAAGGAAAAAGCAAAAAGCAGGGAATAAGGTTAAAGAGATTGTTGCTAACTTTACCGCCCTGAttttt                                     >  2:268701/1‑138 (MQ=255)
                                                                                            tAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                               <  1:409884/140‑1 (MQ=255)
                                                                                              gTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAAt                             <  2:84603/140‑1 (MQ=255)
                                                                                                 gTAGTTTGTGGCGGTAAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAAAAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGcc                          <  1:479834/140‑1 (MQ=255)
                                                                                                     tttGTGGGGGTGATTTTTTCCAGAAAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAACCAACAGGCAGTGAATAGGTGGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                      <  2:304011/140‑1 (MQ=255)
                                                                                                     tttGTGGCGGTGATTTTTTCAAAACAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAAGAAATAACCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                      <  2:687112/140‑1 (MQ=255)
                                                                                                        ggggCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc                   <  1:268701/138‑1 (MQ=255)
                                                                                                                     ttttCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaaa       <  1:337016/140‑1 (MQ=255)
                                                                                                                          cAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa  <  2:597395/140‑1 (MQ=255)
                                                                                                                          aaGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa  <  2:204697/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAA  >  NC_000913/3206315‑3206575

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