Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,334,6880TG32.6% 60.1 / 11.8 46intergenic (+231/‑32)folK/panB2‑amino‑4‑hydroxy‑6‑hydroxymethyldihydropteridine pyrophosphokinase/3‑methyl‑2‑oxobutanoatehydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (11/20);  new base G (14/1);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 3.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGG  >  KT2440/5334560‑5334813
                                                                                                                                |                                                                                                                             
cTATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGtttttt                                                                                                                           >  1:262219/1‑133 (MQ=255)
           ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGaa                                                                                                               >  3:86423/1‑134 (MQ=255)
           ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCTGCAc                                                                                                              >  5:117872/1‑135 (MQ=255)
             gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGCACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTGTTTTAAGCGCCGGAACGa                                                                                                            >  6:98627/1‑135 (MQ=255)
             gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTAAGCGCCTGTACGa                                                                                                            >  4:53092/1‑135 (MQ=255)
             gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTAAAGAGCCTGaaagc                                                                                                            >  1:205531/1‑132 (MQ=255)
             gtgtGGCACCAAAGTGCCGCATACCAGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATTTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTAGATGGCCTGAAAGa                                                                                                            >  2:50791/1‑135 (MQ=255)
              tgtgGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGGTTTTAAGCGCCGGAAccc                                                                                                            >  2:252538/1‑132 (MQ=255)
                    aCCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTc                                                                                                    <  7:129655/136‑1 (MQ=255)
                      aaaaGTGACGCATAACATTATTTAGGCCCGGACGGTCCTGTGCGCGAAAATAAAGAGCGATGAGTTGCCGCTCCGGAATATGACTACACGCGTTGAATTCATTCGTTTTTTACAGCGCCTGAACGAGGAATTTc                                                                                                    <  7:120596/133‑1 (MQ=255)
                         aGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACa                                                                                               <  7:116127/136‑1 (MQ=255)
                           tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                              <  2:300455/135‑1 (MQ=255)
                           tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                              <  4:86423/135‑1 (MQ=255)
                           tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                              <  3:53092/135‑1 (MQ=255)
                           tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                              <  5:257604/135‑1 (MQ=255)
                                cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACATCGCCTGAACGCGGTATTTACTAAATTCCTgc                                                                                        >  8:11509/1‑135 (MQ=255)
                                cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAAAGAGGAATTTCCTACATGCCt                                                                                          >  4:147684/1‑134 (MQ=255)
                                        tATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAAccc                                                                                <  8:143601/136‑1 (MQ=255)
                                          tgtcggCCCGGCCTGACCGGGGCCCGAACATCACGCGCGCTGCTGTGCCGCTCCGGAAGAGGACGACACGCGTTGAATGCAGTCGTGTTGTACCGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                               <  6:193936/131‑1 (MQ=255)
                                           ttAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                               <  4:292023/134‑1 (MQ=255)
                                           ttAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                               <  6:209943/134‑1 (MQ=255)
                                            tAGGCCCGGACGGAGCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCtt                                                                              <  4:3594/134‑1 (MQ=255)
                                             aGGCCCGGCCGGACCTGTGACCGAACATCACCCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGGTGAATGCAGTCGTTTTTTTCACAGCCGGAACAGGTAACTTCGTACCTGCCGGCAGTAACCCttt                                                                             >  1:256418/1‑133 (MQ=255)
                                             aGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTa                                                                             <  1:252538/134‑1 (MQ=255)
                                                  cggacggaCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCAccc                                                                       >  6:77916/1‑135 (MQ=255)
                                                       ggacCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTGCGCGCCGGTACGAGGAATTTCCGACAATCCTGTAGTAACCCTTACCAcccaaaaag                                                                 >  8:82357/1‑130 (MQ=255)
                                                         acCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGc                                                               <  8:292932/136‑1 (MQ=255)
                                                                 ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTAAGCGCCAGAAACACGAATTTACTACCTGCCTGAAGTAACACTTACCACCCATAAACGCCGCAAgg                                                        >  4:54514/1‑135 (MQ=255)
                                                                           cgcgcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTAAAGCGCCTGCACGAGGAATTTCCTACATGCCTGAAGTAACCCTTTCCACCCGTAACCGCCGCAAGGCCAAGGGGGa                                              >  8:137623/1‑135 (MQ=255)
                                                                            gcgcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGGCGGAACGCGGGATTTTCTAACTGCCGGAAGTTAACCGTTCCCCCCCTAAAGGGCGCCAGCCCAAgggggc                                              >  6:227163/1‑131 (MQ=255)
                                                                             cgcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTAAGCCGGCGGAACGAGGGATTTCGTAAATGCCTGAAGTACACCTTACCACCCTTTAAGGCCCCAAGGCCAAggggtaa                                             >  6:77927/1‑130 (MQ=255)
                                                                             cgcgATGATGCGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTAAGCGCCGGAACGAGGAATTTCCTACATGCCAGAAGTAACCCATAACACCCATAACGGCCGCAAGGCCAAGGGGGaa                                             >  3:230403/1‑134 (MQ=255)
                                                                              gcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTCAGCGCCCGGACGAGGCATTTCCTACATGCCTGAAGTAACCCATTCAACCCCTAACGCCCGCAAGGGCCAggggggaa                                            >  2:315991/1‑129 (MQ=255)
                                                                              gcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTAAGCGCCGGAACGAGGAATTTGCTACATGCCGGAAGTAAACCATACCACCCCTAACGGCCGCAACGCCAAGGGGGaaa                                            >  5:87920/1‑134 (MQ=255)
                                                                                gatgatGTTCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCATGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGaaaaaa                                         <  1:50791/135‑1 (MQ=255)
                                                                                gatgatGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGaaaaaa                                         <  5:26302/135‑1 (MQ=255)
                                                                                    atgTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCAGAACCAGGAATTTCCGACATGCCTGAAGTAACCCATAACACCCGAAACCGCCGCAAGGCCCAGGGTGAAAAAATCAc                                     >  4:210829/1‑135 (MQ=255)
                                                                                       tgCCGCTCCGGAAGATGACTACACGCGGTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATg                                 <  8:317545/136‑1 (MQ=255)
                                                                                        gCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATg                                 <  2:262219/135‑1 (MQ=255)
                                                                                        gCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTAGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATg                                 <  3:161888/135‑1 (MQ=255)
                                                                                        gCCGCTCCGGAAGATGACTACACGAGTTGAATGCAGTCGTTTTTTAGAGAGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATg                                 <  2:256418/135‑1 (MQ=255)
                                                                                               cGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCCGCGCCTGGACGCGGAATTTCCTACATGGCTGAAGTTACCCTTACCACCCTTAACGGCCGCAAGGCCCAGGGTGATAAAATGACCATGCTGAccg                          >  6:109958/1‑134 (MQ=255)
                                                                                               cGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCCCCGGCACGAGGACTTTCCTACATGCCGGAAGGAACCCATACCACCCCTAACGGCCGCCACGCCAAGGGGGAAAAAAACACCAAGCAGAccc                          >  2:258092/1‑134 (MQ=255)
                                                                                                          tACACGCGTTGAATGCAGTCGTTTTTTTAATCCGCGGCAACAGGAAGTTCGTAAATGCCATTATTTACCCTTACAACCCAGAAAGGCCACAACCCCAAGGGGGAAAAAAATAACATGCAGAACAGCAACGAcccg               >  4:197995/1‑135 (MQ=255)
                                                                                                            cacGCGTTGAATGCAGTCGTTTTTTAAAGCGCCTGAACGAGGAATTTCCTACATGCCAGAAGTAACCCTTACCAACCATAAAGGCCTCAAGGCCAAGGGTGCAAAAATCACCATGCTGACCGGCGAcgacccgacc            >  8:264555/1‑136 (MQ=255)
                                                                                                                     aaTGCAGTCGTTTTTTTCACCGCCTGAACCAGGTATTTACTACATGCCTGAACTAACCCGTACCACCCGTAACCGCCGCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCAGCGAACACCCGACCTTCGCCaa    >  4:310773/1‑135 (MQ=255)
                                                                                                                       tGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAgg  <  3:197995/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                             
CTATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGG  >  KT2440/5334560‑5334813

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