Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,0900AG46.3% 40.4 / 11.0 41I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (6/16);  new base G (19/0);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 7.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGG  >  NC_000913/2095942‑2096239
                                                                                                                                                    |                                                                                                                                                      
gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt                                                                                                                                                      >  1:867297/1‑151 (MQ=255)
gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt                                                                                                                                                      >  1:888777/1‑151 (MQ=255)
gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGAGCGACATTCCGTTTTGTGTGGCGGGTGGGAtt                                                                                                                                                      <  2:1047027/151‑1 (MQ=255)
  cGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt                                                                                                                                                      >  1:310652/1‑149 (MQ=255)
          aaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTgg                                                                                                                                           >  1:95133/1‑151 (MQ=255)
           aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTGA‑GTatcttg                                                                                                                                           >  2:1053710/1‑148 (MQ=255)
                 acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGc                                                                                                                                    <  1:526183/151‑1 (MQ=255)
                     tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GGTTCTGGAAGATGCcgc                                                                                                                                 >  2:806671/1‑150 (MQ=255)
                        agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGaaa                                                                                                                             <  1:901473/151‑1 (MQ=255)
                               gcTGACGAACTGGTGTTCTACGATATCCCCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTGGGAGCAGTCGCGAAAATTCttt                                                                                                                      >  2:790273/1‑151 (MQ=255)
                                    cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCttt                                                                                                                 <  2:172591/151‑1 (MQ=255)
                                                 tACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGTGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATa                                                                                                       <  2:146816/148‑1 (MQ=255)
                                                   cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTGA‑GTCCCTGGGAGATGCCGCGAAAATTTTTTCCTTTGGCGCGGGTaaaaat                                                                                                  >  2:120450/1‑149 (MQ=255)
                                                    gATATCACCGCTTCGAGCGATGGACGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGACGTGATCGACATTCCGTTTTGTGGGGCGGGGGGGAGTAAGGTAGGTGG‑AGAAGACGCGGAATATCTTTTCCTTGGGGCGGAAAAAAt                                                                                                   >  2:179144/1‑149 (MQ=255)
                                                       atCACCACTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCCCGCGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑TAAACCGGCAGATGCCCCGTAAAATCTTTTCTTTGTGGCGGagacaaatt                                                                                                 >  2:152562/1‑140 (MQ=255)
                                                           ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                                                                                          <  2:231333/151‑1 (MQ=255)
                                                              cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTGGGAGGAGCCGCGAAAAATTTTTCCTTTGGGGCGGATAAAATTTTCATCaa                                                                                          >  2:27837/1‑148 (MQ=255)
                                                                 ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTcc                                                                                    <  1:1000333/151‑1 (MQ=255)
                                                                      gATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTTCATCAACTCTcc                                                                                    >  1:200356/1‑146 (MQ=255)
                                                                            cGTGTGGGAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGTGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTACCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGc                                                                           <  2:931210/149‑1 (MQ=255)
                                                                              tgtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAA‑GGCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTGGCGCGGGGAAAAATTTCATCAACCCTCCTGCGCTGGCaaa                                                                        >  2:55671/1‑148 (MQ=255)
                                                                               gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTTGGATTAA‑GTGTCTGGAAGATGCCCCGGAAATTTTTTCCTTTTGCGCGGGTAAAATTTCCCTCAACTCTCCTGCGGTGGccaacc                                                                      >  2:379403/1‑146 (MQ=255)
                                                                               gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GGCTCTGGGAGAAGCCCCGAAAATTCTTTTCTTTGGGGCGGATAAAATTTCCCTCAACTCTCCTGCGCTGGGAGAcc                                                                      >  2:135556/1‑151 (MQ=255)
                                                                                tgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GTCTCTGGGAGAAGCCGCGAAAATTCTTTTCTTTTGCGCGGATAAAATTTCCCTCAACCCTCCTGCGCTGGCCGAccc                                                                     >  2:692838/1‑151 (MQ=255)
                                                                                      aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GGCCCTGGAAAAAGCCGCGGAAATTGTTTCCTTTTGGGCGGATAAAATTTCCATCAACTCTCCCGCGCCTGCAGACCCaa                                                                   >  2:120164/1‑147 (MQ=255)
                                                                                      aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatt                                                               <  1:696540/151‑1 (MQ=255)
                                                                                        aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatta                                                              <  1:338209/150‑1 (MQ=255)
                                                                                               gggTATCTCGCTTGTCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTGCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACt                                                        <  2:798329/149‑1 (MQ=255)
                                                                                                    tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCCTCAACTCTCCTGCGCTGGCAGACCCAACATTAATTAATCGCCTgg                                                 >  1:1054823/1‑151 (MQ=255)
                                                                                                         cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGAt                                            <  2:888175/151‑1 (MQ=255)
                                                                                                              cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCtt                                       <  1:630962/151‑1 (MQ=255)
                                                                                                              cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCtt                                       <  2:937948/151‑1 (MQ=255)
                                                                                                               ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAA‑GACTCTGGAAGATGCCCCCAAAATTTTTTCCTTTGGCGGGGAGAAAAATTTCCTCAACTCTCCTGCGGTGGGAGACCCCAAATTAATTACTTGCCTTGCCGATCGCttt                                      >  2:683110/1‑151 (MQ=255)
                                                                                                                gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAA‑GTCTCTGGAAGATGGCCCGAAAATTCTTTTCTTTGGCGCGGTGAAAATTTCCATCAAATCTCCCTCGCTGGCAGACCCAACATTTATTACTCGCCTTGGCGAACGCctt                                      >  2:935195/1‑150 (MQ=255)
                                                                                                                gAAGTGATCGAAATTCCGTTTTGTGTGGCGGGCGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGGTAAAATTTCCATCAAATCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTCGCCGATCGc                                         >  2:97883/1‑147 (MQ=255)
                                                                                                                           cATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtgt                          <  2:928614/151‑1 (MQ=255)
                                                                                                                            aTTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCCGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAg                             >  2:882687/1‑147 (MQ=255)
                                                                                                                                     tgtgtgGCGGGTGGGATAAA‑ATCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTAATCCCCTGGCCGCTCGCTTTTGCGTGCAGTGTATTGTGGTc                 >  2:84532/1‑150 (MQ=255)
                                                                                                                                           gCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTg          >  2:263897/1‑151 (MQ=255)
                                                                                                                                                 gggATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt    >  2:274806/1‑151 (MQ=255)
                                                                                                                                                   gATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTgg  <  1:183175/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGG  >  NC_000913/2095942‑2096239

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