Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT34.4% 35.3 / 18.9 32intergenic (‑179/‑28)ttdR/ttdADNA‑binding transcriptional activator Dan/L(+)‑tartrate dehydratase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/19);  new base T (11/0);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 6.05e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAA  >  NC_000913/3206308‑3206574
                                                                                                                                 |                                                                                                                                            
cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACa                                                                                                                                  <  2:1223400/140‑1 (MQ=255)
   tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                   <  2:186378/136‑1 (MQ=255)
   tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                   >  1:186378/1‑136 (MQ=255)
    cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAAAAAAAACaaaaa                                                                                                                              >  2:166440/1‑140 (MQ=255)
    cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCCGAcaaaaacaaaaa                                                                                                                              >  2:1175533/1‑140 (MQ=255)
       tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTGATTTTTTCATAcaaaaacaaaaaaaa                                                                                                                           >  2:458361/1‑137 (MQ=255)
         ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTAGAcaaaaacaaaaaattga                                                                                                                         >  2:876945/1‑135 (MQ=255)
          tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAg                                                                                                                        <  2:567911/140‑1 (MQ=255)
             ggCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCCAAAAAAAAAAAAAAATGGAGTTg                                                                                                                     >  1:751120/1‑140 (MQ=255)
              gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGc                                                                                                                    <  2:1006738/140‑1 (MQ=255)
               cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTCTACAAAAAAAAAAAATTGAGTTGcc                                                                                                                   >  2:458699/1‑140 (MQ=255)
                   atatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAAAGGAGTTGCCatga                                                                                                               >  1:393846/1‑140 (MQ=255)
                       tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAg                                                                                                           <  1:1175533/140‑1 (MQ=255)
                          ttcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                        <  1:876945/140‑1 (MQ=255)
                                ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaa                                                                                                  <  1:988878/140‑1 (MQ=255)
                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                                 <  1:471548/140‑1 (MQ=255)
                                     gcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGGTGGCATGATGAGCGGAAGgaataaaa                                                                                             >  2:912306/1‑138 (MQ=255)
                                           aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                       <  1:239746/140‑1 (MQ=255)
                                             tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                     <  1:458361/140‑1 (MQ=255)
                                                   aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAAt                                                                              <  1:1126699/140‑1 (MQ=255)
                                                     gacagaTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATaa                                                                            <  1:329101/135‑1 (MQ=255)
                                                        tAATTAACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCCAACAAAAAAAAAAAATGGAGTTTGCATGATGAGCGAAAATAAAAAACAAAAAGCAGGGAaaaaat                                                                          >  2:351791/1‑138 (MQ=255)
                                                                atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGTATTTTTTCATACAAATAAAAAAAAAAGAGTTTTCCTGATGAGCGGAAGGAATAAACCAAAGGGAGGGAAAAAGATTTCagag                                                                  >  2:463288/1‑140 (MQ=255)
                                                                 tatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgaga                                                                 <  1:166440/140‑1 (MQ=255)
                                                                  atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGGGTTGCCATGATGAGCGAAAAGAAAAAAAAAAAAGCAGGGAATAAATTGACAgagaa                                                                >  2:295853/1‑139 (MQ=255)
                                                                  atTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAACAAAAAAAGGAGTTGGCATTATGAGGGAAAGAAATAAAAAAAAAGCAGTGAATAAATTGACAgagaa                                                                >  2:1315982/1‑139 (MQ=255)
                                                                      aCGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACAAAAAAAAAAAATGGAGGTGGCCTGATGAGGGAAAAGAAAAAAAAAAAAGCAGGGAAAAAAATGACAGAGACTGTc                                                            >  2:329101/1‑140 (MQ=255)
                                                                        ggTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATAAAAAAAAAGGAGGTGCCCTGATGAGGGAAAGTAATAAACAACAAGCAGGGAATAAGTTGACAGAGATTGTTGc                                                          >  2:1133432/1‑140 (MQ=255)
                                                                          tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAATGGAGTTGGCATGATGGGCGGAAAAAAAAAACAAAAGGCAGTGAATAAATTGACAGAGATTGTCGCTa                                                        >  2:1198350/1‑140 (MQ=255)
                                                                            gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAAATGGAGGTGGCATGATGGGCGAAAAAAAAAAACAAAAAGCAATGAATAAATTGAAAAAGAGTGTTGCTaaa                                                      >  2:1009296/1‑139 (MQ=255)
                                                                                     ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGcc                                             <  1:1315982/140‑1 (MQ=255)
                                                                                     acaaaaCCACATTTTTAGTGGGAGTTTGCGGCGTAGATTTTTTCCAGCCAAAAACAAAACTTGGAGTTGCCGAGAAGACCGAAAGAAATAAGCAAGAGGCAGTAAATAGGTTGACAGAAATTGTCGAAAATTTTCCCGcc                                             <  2:374525/138‑1 (MQ=255)
                                                                                      taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGGGGTGGCCTGATGAGGGGAAAGAAAAAAAAACAAGCAGTGAAAAAAGTTGCAGAGATTTTTGGTAAATTTTCccccc                                            >  2:1308116/1‑139 (MQ=255)
                                                                                      aaaaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCa                                            <  1:1308116/139‑1 (MQ=255)
                                                                                               aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAACAGGAGTTTGCCTGATGGGCGAAAGGAAAAAACAAAAAGCAGTGAATAAATTGACAGAGATTGTTGCTTACTTTAACCCCATGAtttttt                                   >  2:913562/1‑137 (MQ=255)
                                                                                                    tAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGa                              <  2:1205927/140‑1 (MQ=255)
                                                                                                                 gggCGGTAAATTTTCCCAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACg                 <  1:913562/139‑1 (MQ=255)
                                                                                                                         aaaTTTTTCCAGACAAATAAAAAAACGGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa          <  1:175444/139‑1 (MQ=255)
                                                                                                                            tttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAAAAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa       <  1:351791/140‑1 (MQ=255)
                                                                                                                                 ccAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaa  <  1:883032/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                            
CATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAA  >  NC_000913/3206308‑3206574

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: