Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,334,6880TG50.0% 9.5 / 13.6 26intergenic (+231/‑32)folK/panB2‑amino‑4‑hydroxy‑6‑hydroxymethyldihydropteridine pyrophosphokinase/3‑methyl‑2‑oxobutanoatehydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/10);  new base G (13/0);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCG  >  KT2440/5334566‑5334800
                                                                                                                           |                                                                                                                 
gCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGcc                                                                                                       >  4:72819/1‑135 (MQ=255)
gCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGcc                                                                                                       >  3:62024/1‑135 (MQ=255)
   tCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGGTTTTCAGCGcccgg                                                                                                    >  8:50165/1‑132 (MQ=255)
     ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTCAGCGCCTGAAc                                                                                                  >  8:37150/1‑135 (MQ=255)
     ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGCCTGaa                                                                                                   >  3:121693/1‑134 (MQ=255)
      cGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTGAGCGCCCGAGCg                                                                                                 >  1:95510/1‑135 (MQ=255)
        tgtgGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAg                                                                                               <  4:62024/135‑1 (MQ=255)
           ggCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTAAAGCGCCGGAACGAggta                                                                                            >  7:66649/1‑133 (MQ=255)
               ccAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTc                                                                                        <  2:95124/135‑1 (MQ=255)
                 aaaGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCt                                                                                      <  4:121693/135‑1 (MQ=255)
                   aGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTAAACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCAGAACGAGGAATTtactac                                                                                    >  4:66509/1‑135 (MQ=255)
                     tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                  <  1:159702/135‑1 (MQ=255)
                            tACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGaa                                                                           <  5:66397/135‑1 (MQ=255)
                               ccGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGt                                                                         <  8:180287/134‑1 (MQ=255)
                                            cggacggaCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTACGCGCCGGTACCAGGAATTTGCTACCTGCGCGAAGTAAACCATaccaac                                                            >  5:20693/1‑132 (MQ=255)
                                            cggacggaCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCGGCACGATGTAGTTGCTACCTGCCGGAAGTAACCCGTACCAccc                                                           >  5:50123/1‑135 (MQ=255)
                                                    ccTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGc                                                  <  8:27955/135‑1 (MQ=255)
                                                           ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCTGAACGAGGAATTTGCTACATGCCTGAAGTAAACCTTACCACCCTTAACCGCCGCAAgg                                           >  7:193760/1‑135 (MQ=255)
                                                            cGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGc                                          <  6:156463/135‑1 (MQ=255)
                                                                cATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGCCTGGACGAGGAATTTCCGAACTGCCTGAAGTAACCCTTACCACCCTTAACGGCCAAAAGGCCaac                                      >  6:100219/1‑134 (MQ=255)
                                                                  tCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCCGGACGCGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCATAACGGCCGCAAGGCCAAggg                                    >  2:73401/1‑135 (MQ=255)
                                                                          gatgatGTGCCGCTCCGGAAGATGACTACACGCGTTGACTGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGaaaaaa                            <  6:163573/135‑1 (MQ=255)
                                                                                  gCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATAACCATg                    >  6:105667/1‑135 (MQ=255)
                                                                                    cGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCt                  <  7:91671/135‑1 (MQ=255)
                                                                                         cGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGCAGTAACCCTTTCCACCCGTAACGGCCGCAAGGCCAACGGTGAAAAAATCACCATGCTGAcca             >  7:24535/1‑134 (MQ=255)
                                                                                                    tACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCTGCAAGAGGAATTTACTACCAGCCTGAAGTAACCCATACCACCCATTAACGCCGCAAGGCCCAGGGGGAAAAAAAAACCAAGCTGACCGGCGacgccccg  >  2:199906/1‑130 (MQ=255)
                                                                                                                           |                                                                                                                 
GCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCG  >  KT2440/5334566‑5334800

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: