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

TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206315‑3206582
                                                                                                                                 |                                                                                                                                          
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGGATTTTTTCAaaaaaaaaaaaaaaaaa                                                                                                                                  >  2:914099/1‑137 (MQ=255)
tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTTACACAaaaaaaaaaaata                                                                                                                                  >  2:425925/1‑137 (MQ=255)
   tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                               <  2:655962/140‑1 (MQ=255)
     tGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGAGtt                                                                                                                             >  1:671050/1‑140 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGAAAAAAAAAAAAAATGGGGGTGc                                                                                                                           >  2:936524/1‑140 (MQ=255)
       gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGc                                                                                                                           <  1:732970/140‑1 (MQ=255)
        cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGcc                                                                                                                          <  2:792675/140‑1 (MQ=255)
         gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                         <  1:494073/140‑1 (MQ=255)
                 cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGGGGTGTCCTGATGAgg                                                                                                                 >  1:1172808/1‑139 (MQ=255)
                 cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAAGAGGTGGCATGATGAgg                                                                                                                 >  1:1087643/1‑139 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                               <  1:936524/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                               <  1:425925/140‑1 (MQ=255)
                   ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                               <  1:896538/140‑1 (MQ=255)
                     cTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaaa                                                                                                             <  2:671050/140‑1 (MQ=255)
                     cTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAATTTTTCAGAAAAAAAAAAAAAATAGAGGTGGCATGATGGGGGaaa                                                                                                             >  2:987361/1‑140 (MQ=255)
                      ttATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAg                                                                                                            <  1:987361/140‑1 (MQ=255)
                                   gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAA‑AAAAAAAATGGAGTTTCCATGATGAGGGAAAAGAATAAacaacaag                                                                                              >  1:660656/1‑138 (MQ=255)
                                    aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                              <  2:320968/140‑1 (MQ=255)
                                    aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                              <  2:790907/140‑1 (MQ=255)
                                     ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                             <  1:900333/140‑1 (MQ=255)
                                          gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGa                                                                                        <  2:789505/140‑1 (MQ=255)
                                               gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAg                                                                                   <  1:32063/140‑1 (MQ=255)
                                                  gacagcATATTCACGGTGGCAAAAAATATAAAACCACATTTTTGGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTg                                                                                <  1:1001171/135‑1 (MQ=255)
                                                        atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCCGACAAAAAAAAAAAATGGGGGTGGCATGATGAGGGAAAAAAAAAAACAACAAGCAGTGAAAAAAATGACagaa                                                                          >  2:32063/1‑139 (MQ=255)
                                                          atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATGATGAGGGAAAGTAATAAACAAAAAGCAGAGAATAAGTTGACAGAGAt                                                                        >  2:932262/1‑140 (MQ=255)
                                                                  tGGCAAAAAATATAAAACCACATTTTTGGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAATGGAGGTGGCATGATGAGGGAAAAAAAAAAACAAAAAGCAGTGAATAAGATGACAGAGAATGTctgtt                                                                >  2:1001171/1‑136 (MQ=255)
                                                                  tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGGGGTGCCATGATGAGGGAAAGAAAAAAAAAAAAGGCAATGAATAAATTGACAGAGATTTTCGCTa                                                                >  2:630075/1‑140 (MQ=255)
                                                                   ggaaaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTaa                                                               <  1:129292/137‑1 (MQ=255)
                                                                   ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGAAAGTAATAAACAAAAAGCAGAGAATAAGTTGACAGAGATTGTCGCTaa                                                               >  2:400424/1‑140 (MQ=255)
                                                                              taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGCCATGATGAGGGGAAGAAATAAAAAACAAGCAGTGAAAAAGTTGACAGAGAATGTTGCTAACTTTACCGccc                                                    >  2:4663/1‑139 (MQ=255)
                                                                                       aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGGCATGATGAGCGAAAATAATAAGCAAAAAGCAGTGAATAAATTGACAGAGATTTTTGCTAACTTTTCCGCCATGAtttttt                                           >  2:1054278/1‑137 (MQ=255)
                                                                                            tAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                                      <  1:1035839/140‑1 (MQ=255)
                                                                                             aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGaa                                     <  1:780465/140‑1 (MQ=255)
                                                                                             aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                                     <  1:3852/140‑1 (MQ=255)
                                                                                               gggTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATg                                   <  2:1087643/139‑1 (MQ=255)
                                                                                                   aGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGGGTTGCCATGATGGGGGGAAAGAAAAAACAACAAGCAGTGAATAAATTGACAGAGAATGTCGCCAACTTTTCCGCCCTGATTTTTTCCCGAAAGCCtg                               >  2:1035839/1‑140 (MQ=255)
                                                                                                           gCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                       <  1:1054278/140‑1 (MQ=255)
                                                                                                              gTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGCCATGATTAGCGAAAGTAATAAACAACAAGCAGTGAATAAATTGACAGAGATTGTCGCTAACTTTACCGCCCTGATTTCTTCCCGAAAGCCTGATGACgtgggg                    >  2:172711/1‑138 (MQ=255)
                                                                                                                aattttttCCAGAAAAAAACAAAAACTGGAGTCGCCATGAGGAGAGAAAGTAAAAACCAACAGGCAGTGAATAGGTTGACAGATTTTGTCCCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa                  <  2:1172808/137‑1 (MQ=255)
                                                                                                                                aTACAAAAACTGGTGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGCTAAACTAAAACAGCTa  >  2:492042/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                          
TGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTA  >  NC_000913/3206315‑3206582

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: ‑