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

CTCACGACTATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGC  >  KT2440/5334553‑5334808
                                                                                                                                        |                                                                                                                         
cTCACGACTATAGGCGTCCCGTGTGGAACCAAAGTGCCGCATAACCGTATTAAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGAAGATGTGCAGCTCCGGAAGATGCGTACACGCGTTGATTGCAGACGt‑t                                                                                                                           >  6:217074/1‑136 (MQ=255)
                tCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTACAGCGCCTGaa                                                                                                           >  8:192264/1‑136 (MQ=255)
                  ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTCAGCGCCTGAACg                                                                                                         >  5:99247/1‑136 (MQ=255)
                    gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGCCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTAGAGCGCCGGAAccgc                                                                                                       >  6:149930/1‑133 (MQ=255)
                    gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCGCGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAg                                                                                                       >  4:140562/1‑136 (MQ=255)
                    gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGAACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTGTTTGAGCGCCGGAACggg                                                                                                       >  5:143066/1‑134 (MQ=255)
                     tgtgGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTAAAGCGCCAGGACGatg                                                                                                      >  1:231397/1‑134 (MQ=255)
                     tgtgGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GGTTTTGAGCGCCGGCACGCgg                                                                                                      >  5:153849/1‑136 (MQ=255)
                                 gTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                          <  3:216880/136‑1 (MQ=255)
                                 gTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCAGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                          <  8:99585/136‑1 (MQ=255)
                                      gCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTg                                                                                     <  7:116144/136‑1 (MQ=255)
                                       cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTg                                                                                     >  2:52042/1‑135 (MQ=255)
                                       cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTAAAGGGCCGCGACCAGGGATTTGCTAGATTCCTGa                                                                                    >  5:165565/1‑136 (MQ=255)
                                       cATACCCGGATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGa                                                                                    <  6:242845/136‑1 (MQ=255)
                                         tACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAg                                                                                  <  5:18776/136‑1 (MQ=255)
                                          aCCCGTAGTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTAGAATGCAGTCGT‑GTTTTACAGCCCGGGAACCCGCAATTTCCCACATGCCAGaac                                                                                  >  2:3811/1‑134 (MQ=255)
                                             cGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGGTTTTAC‑GCGGCCGGACGCTGAATTTACTAAATGGCTGCAGTAAc                                                                              >  2:123570/1‑136 (MQ=255)
                                               tATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTCAGCGCCGGCACGGGGAATTTGCTACATGGCTGCAGGTAccc                                                                            >  5:33875/1‑136 (MQ=255)
                                                 tttAGGCCCGGACGGACCTGTGCCCGAACATCTCGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                           <  3:140562/135‑1 (MQ=255)
                                                 tttAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                           <  5:39546/135‑1 (MQ=255)
                                                  ttAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCtt                                                                          <  1:122454/135‑1 (MQ=255)
                                                       cccGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTaccacc                                                                    >  8:103011/1‑136 (MQ=255)
                                                         cggacggaCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTTCAGCGCCGGCACGCGGAATTTCCTACATGCCTGAAGTAACCCTTAccccccg                                                                  >  8:110976/1‑135 (MQ=255)
                                                          ggacggacCTGTGCCCGAACATGACGCGCGCTGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCtt                                                                 <  3:161427/136‑1 (MQ=255)
                                                          ggacggacCTGGGCCCGACAATCCCGCGCGATGATGGGCCGCTCCGAAAGAGACCTACACGCGTTCAATGAAGTCGT‑TTTTTACAGCGCCTTAACTAGTAATTTCCTAAATCCATGAATTAACCGTTACCACCCta                                                                 <  8:54048/136‑2 (MQ=255)
                                                           gacggacCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTa                                                                <  2:55845/136‑1 (MQ=255)
                                                            acggccCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTAAAGCGCCTGAAGGAGGAATTTGCTACAAGCCTTAAGTTACCCTTAACACCCGTaa                                                               >  3:205000/1‑136 (MQ=255)
                                                              ggacCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTACTACATGCCTGAAGTAACCCTTACCACCCTTAACg                                                             >  6:114091/1‑136 (MQ=255)
                                                                   tgtgCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTCAGCGCCTGAACGAGGCATTTCCTACATGCCGGAAGTAACCCTTACCACCCTTAACGGccg                                                         >  3:114800/1‑134 (MQ=255)
                                                                             cATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTCAGCGCCTGAACGCGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCATAACCGCCGCACAGCCCAAg                                              >  8:115791/1‑136 (MQ=255)
                                                                                     gcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGaaaaa                                     <  1:213507/136‑1 (MQ=255)
                                                                                              tgCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATg                            <  6:63319/136‑1 (MQ=255)
                                                                                                ccGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCt                          <  7:103011/136‑1 (MQ=255)
                                                                                                      cGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTCAGCGCCGGAACGAGGCATTTACTACAAGCCGTAAGTAACCCATACCACCCATAACCGCCACA‑AGGCCAAGGGGGAAAAAATCACCAAGCGGAccag                    >  8:116753/1‑134 (MQ=255)
                                                                                                      cGGAAGATGACTACACGCGTTGAATGCAGTCGT‑GTTTTTAAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCGTAACGGCCGCA‑AGGCCAAGGGGGAAAAAATCACCAAGCTGACCTg                    >  1:35335/1‑136 (MQ=255)
                                                                                                           gATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTacg               <  7:106076/136‑1 (MQ=255)
                                                                                                           gATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTacg               <  3:26335/136‑1 (MQ=255)
                                                                                                                cTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGAcgcg          <  2:35335/136‑1 (MQ=255)
                                                                                                                   cacGCGTTGAATGCAGTCGT‑GTTTTAAAGCGCCTGAACGAGGAATTTACTAAATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGccc       >  8:184982/1‑136 (MQ=255)
                                                                                                                      gcgTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTc    <  7:192264/136‑1 (MQ=255)
                                                                                                                        gTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGc  <  3:133659/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                         
CTCACGACTATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGT‑TTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCA‑AGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGC  >  KT2440/5334553‑5334808

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