Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT32.5% 50.9 / 21.3 40intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/25);  new base T (13/0);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 8.73e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NC_000913/3206297‑3206566
                                                                                                                                           |                                                                                                                                   
cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTcc                                                                                                                                    <  2:511094/140‑1 (MQ=255)
   aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGa                                                                                                                                 >  2:1023783/1‑140 (MQ=255)
             ttCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACaaa                                                                                                                       >  1:1201166/1‑140 (MQ=255)
               cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaaa                                                                                                                     >  1:490584/1‑140 (MQ=255)
                ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAAc                                                                                                                    <  2:490584/140‑1 (MQ=255)
                ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGATTGTGGCGGTGAATTTTTCCAGACAAATACAAAAAc                                                                                                                    <  1:1135898/140‑1 (MQ=255)
                  tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCACACAAAAACAAAAAATg                                                                                                                  >  2:471458/1‑140 (MQ=255)
                  tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTg                                                                                                                  <  1:478333/140‑1 (MQ=255)
                    ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                <  2:173957/140‑1 (MQ=255)
                         gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCACACAAATACAAAAAATGGAggttg                                                                                                           >  1:692250/1‑136 (MQ=255)
                              atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                      <  1:153669/140‑1 (MQ=255)
                              atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                      <  1:471458/140‑1 (MQ=255)
                                   cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGTTGGCATGATGAgg                                                                                                 >  1:1179409/1‑139 (MQ=255)
                                     ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                               <  2:975208/140‑1 (MQ=255)
                                          aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAAAGGAGGTGGCAAGATGAGCGaaaaga                                                                                          >  2:352521/1‑137 (MQ=255)
                                          aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAAAAAAAAAAAAAAAATGGGGGTTCCATGATGAGCGgaaaaa                                                                                          >  2:290786/1‑137 (MQ=255)
                                            tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                        <  1:197094/140‑1 (MQ=255)
                                            tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                        <  1:376502/140‑1 (MQ=255)
                                                      aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                              <  1:54272/140‑1 (MQ=255)
                                                       ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                             <  2:1260558/140‑1 (MQ=255)
                                                           cGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTg                                                                         <  1:130469/140‑1 (MQ=255)
                                                                tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGCAATTTTTCAGACAAAAACAAAAAATGGAGGTTCCCTGGTGAGCGAAAATAATAAACAAAAACCAATGAATaa                                                                    >  2:130469/1‑140 (MQ=255)
                                                                  ttAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGTAATTTTTCCTGCAAAAACAAAAAATTGAGGTGTGATGATGAGCGAAAATAATAAACAAAAAGCAGGGGATaaat                                                                  >  2:732972/1‑138 (MQ=255)
                                                                          atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAAACAAAAAAAAAAAATGGAGGTTCCCTGTTGAGGGGAAATAAAAAACAAAAAGCCGTGAAAAAAtttgaagag                                                          >  2:712829/1‑133 (MQ=255)
                                                                          atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGGCATGATGAGCGAAAAAAAAAAACAAAAGGCAGTGAATAAAATGACagaa                                                          >  2:143172/1‑139 (MQ=255)
                                                                           tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                         <  2:376788/140‑1 (MQ=255)
                                                                            atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAAAGGAGGTGCCCTGATGAGGGAAAGGAAAAAAAAAAAAGCAGGGAAAAAATTGaaaaaaaa                                                        >  2:1089430/1‑133 (MQ=255)
                                                                                   gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTGGCCATGATGGGGGAAAGGAAAAAGCAAAAAGCAGAGAATAAGATGACAGAGATTGTCGCt                                                 >  2:1146766/1‑140 (MQ=255)
                                                                                   gTGGAAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATCCAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCt                                                 <  2:546744/140‑1 (MQ=255)
                                                                                          aaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                          <  1:143172/140‑1 (MQ=255)
                                                                                                taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTTGCCATGATGAGCGAAAGAAAAAAACAACAAGCAGGGAAAAAATTGACAAAGATTGTCGCTAACTTTACCGccc                                    >  2:159819/1‑139 (MQ=255)
                                                                                                    aCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                <  1:352521/140‑1 (MQ=255)
                                                                                                    aCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                <  1:159819/140‑1 (MQ=255)
                                                                                                    aCCACATTTTTAGTGGTAGTTTGGGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                <  1:1245751/140‑1 (MQ=255)
                                                                                                      cacaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAttt                              <  1:599936/140‑1 (MQ=255)
                                                                                                        caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                            <  1:1023783/140‑1 (MQ=255)
                                                                                                         aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTa                           <  2:1179409/140‑1 (MQ=255)
                                                                                                          tttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTAc                          <  1:290786/140‑1 (MQ=255)
                                                                                                               aGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAGAAAAAAGCAACAGGCAGTGAATAAATTGACAAAGATTGTCGCTAACTTTTCCCCCCTGATTTTTTCCCCaa                     >  2:450078/1‑140 (MQ=255)
                                                                                                               aGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGaa                     <  1:538006/140‑1 (MQ=255)
                                                                                                                gTGGTAGTTGGGGGAGGTGATTTTTCCCAGACAAAAACAAAAGTTGGTTTTCCCATAATGAGAGAAGGTAAAAAGAACCAGGCAGTCAGTAAGTTGACAAAGATTTTCTCTACCTTTCCCGCCATGATTTCTACCCGAAt                    <  2:692250/140‑1 (MQ=255)
                                                                                                                   gTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCATGATGAGGGGAAAGAATAAACAAAAAGCAGTGAATAAATTGACAGAGATTGTCGCTAACTTTTCCGCCATGATTTCTACCCCAATGcc                 >  2:1125265/1‑140 (MQ=255)
                                                                                                                     aGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGGGGTGGCATGATGAGGGGAAATAAAAAACAAAAAGCAGGGAAAAAATTGACAGAGAATGTCGCTTACTTTACCCCCCTGATTTTTACCCGAATGCCtg               >  2:697733/1‑140 (MQ=255)
                                                                                                                     aGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGGCCTGATTAGCGAAAAGAATAAACAAAAAGCAGTGAAAAAATTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTACCCCAATGCCtg               >  2:538006/1‑140 (MQ=255)
                                                                                                                                  gAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa  <  2:366512/140‑1 (MQ=255)
                                                                                                                                  gAATTTTTCAAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa  <  2:1201166/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGA  >  NC_000913/3206297‑3206566

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: