Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,334,6930AT31.2% 62.6 / 20.1 48intergenic (+236/‑27)folK/panB2‑amino‑4‑hydroxy‑6‑hydroxymethyldihydropteridine pyrophosphokinase/3‑methyl‑2‑oxobutanoatehydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (10/23);  new base T (15/0);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 3.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCGCCAGCCA  >  KT2440/5334571‑5334823
                                                                                                                          |                                                                                                                                  
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                                                                                                        <  4:86423/135‑1 (MQ=255)
                tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                                        <  5:257604/135‑1 (MQ=255)
                tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                                        <  3:53092/135‑1 (MQ=255)
                tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                                        <  2:300455/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                                                                                         <  6:209943/134‑1 (MQ=255)
                                ttAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                                         <  4:292023/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)
                                                                                                                     gttttACAGCGCCTGAACGAGGAATTTCCTACCTGCCTGAAGTAACCCATACCACCCTTACCGGCCTCAAGGCCAAGGGTGAAAAAATCCCCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCgccagcca  <  8:112304/135‑1 (MQ=255)
                                                                                                                       tttACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCgccagcca  <  3:310773/134‑1 (MQ=255)
                                                                                                                          |                                                                                                                                  
CCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCGCCAGCCA  >  KT2440/5334571‑5334823

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