Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,0900AG42.5% 52.6 / 11.7 40I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (4/19);  new base G (17/0);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.94e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095950‑2096240
                                                                                                                                            |                                                                                                                                                          
cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTctct                                                                                                                                                  >  1:526713/1‑151 (MQ=255)
  aaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCCCCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTgg                                                                                                                                                <  2:62669/151‑1 (MQ=255)
   aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGa                                                                                                                                               <  1:721105/151‑1 (MQ=255)
     cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTGTCTTGaa                                                                                                                                              >  2:430622/1‑150 (MQ=255)
           gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGAT‑Gcc                                                                                                                                       >  2:616417/1‑150 (MQ=255)
             tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTTGAAAAT‑Gcccc                                                                                                                                     >  2:93128/1‑148 (MQ=255)
               aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGAT‑GCcgcg                                                                                                                                    >  1:407850/1‑149 (MQ=255)
                agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑aaa                                                                                                                                <  1:31194/151‑1 (MQ=255)
                agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑aaa                                                                                                                                <  2:334715/151‑1 (MQ=255)
                       gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCttt                                                                                                                         <  1:717413/151‑1 (MQ=255)
                             gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTg                                                                                                                   <  1:48370/151‑1 (MQ=255)
                                   gggtTCTACGATATCACCGCTTCCATCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑gg                                                                                                            <  2:168053/149‑1 (MQ=255)
                                       tCTACGATATCACTGCTTACAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACGTTCCGTTTTGTGTGGCGGGTAGGACTTAGAATCTATGAGAT‑TCCGCG‑AAAAATTTTTCCTTTGGAGC‑ggg                                                                                                           >  2:623934/1‑147 (MQ=255)
                                       tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATa                                                                                                         <  1:223503/150‑1 (MQ=255)
                                       tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATa                                                                                                         <  1:466538/150‑1 (MQ=255)
                                          aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑GAAATTTTTTCCTTTTGCGC‑GGGTaaaa                                                                                                      >  1:207499/1‑150 (MQ=255)
                                             atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTc                                                                                                  <  2:197312/151‑1 (MQ=255)
                                              tatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGAT‑GCCCCG‑GAAATTTTTTTCTTTGGTGC‑GGGTAAAAttttc                                                                                                 >  2:420649/1‑149 (MQ=255)
                                                  aCCGCTTCCAGCGATGGCCGTGCGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGTGGGTTCAGTCCCTCGAAGTTGGCCACC‑AACTTTTTTTCTTTTTG‑GC‑GGGTACAATTTCCATc                                                                                              >  2:130088/1‑150 (MQ=255)
                                                           aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTCCTTTTGCGC‑GGATAAAATTTTCATCAACTCTcc                                                                                      >  2:76644/1‑149 (MQ=255)
                                                           aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑GAAATTTTTTTCTTTGGCGC‑GGGTAAAATTTCCCTCAACTCTcc                                                                                      >  1:313977/1‑149 (MQ=255)
                                                                    cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGGTCTGGGAGAT‑GCCGCG‑AAAATTTTTTTCTTTGGGGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTgg                                                                              >  2:98932/1‑148 (MQ=255)
                                                                           tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTTAGTCTATGGGAGAT‑GCCGCG‑GAAATTCTTTTCTTTGGCGCGGGACAAAATTT‑CATCAACACTCCCTCGCTGGCAC‑ACCCa                                                                     >  2:514691/1‑149 (MQ=255)
                                                                           tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGGCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAG‑ACCCaa                                                                    >  2:212412/1‑150 (MQ=255)
                                                                             gATAAAAGCTGGGTATCTCGCGTGTCGGAAGTGATCGACATTCCGTTTTGTGGGGCGGGGGGGGCTAGGCCTCTGGGAGAT‑GCCCCG‑AAAATTTTTTGTTTTGGGGC‑GGGCAAAATTTTCCTCAAACCTCCTTGGCTTGCAG‑Accccaga                                                                  >  2:68061/1‑146 (MQ=255)
                                                                                  aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACattaatt                                                            <  1:748791/151‑1 (MQ=255)
                                                                                       gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTc                                                        <  2:747554/150‑1 (MQ=255)
                                                                                       gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCg                                                       <  1:461817/151‑1 (MQ=255)
                                                                                            tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTTCTTTGGCGC‑GGATAAAATTTTCCTCAACTCTCCTGCGCTGGCCG‑ACCCAACATTAATTACTCGCCCgg                                                  >  2:203518/1‑151 (MQ=255)
                                                                                               cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGcc                                                <  1:723992/150‑1 (MQ=255)
                                                                                                gcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGa                                              <  1:231543/151‑1 (MQ=255)
                                                                                                  gTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGACTCTGGGAGAA‑GCCGCGAAAAATTTTTTCTTTTGGGCG‑GGAAAAAATTTTC‑TCCACTCCCCTGCCCTGGGGGAACCCCAC‑TTAATTATTCTCCTCGcccagc                                            >  2:425809/1‑147 (MQ=255)
                                                                                                             gATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGt                                  <  1:612460/150‑1 (MQ=255)
                                                                                                              aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGGCCCTTGAAGAT‑GGCCCG‑ACAATTTTTTCCTTTGGCGG‑GGATAAAAATTCCCTCAACTCTCCTGCGCTGGCAG‑ACCCCACATTAATTTTCCCCCTGGCCCATCGCTTTg                                      >  2:502857/1‑145 (MQ=255)
                                                                                                               tCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTATCTGGGAGAT‑GCCGCG‑GAAATTTTTTTCTTTTGCGC‑GGGTAAAATTTCCCTCAACTCTCCTGCGCTGGCAG‑AACCCACATTTATTACTCCCCTGGCCGATCGCTTTTGCGTg                                 >  2:70789/1‑149 (MQ=255)
                                                                                                                cGACATTCCGTTTTGTTTTGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTTGCCGATCGCTTTGGCGTGCa                               <  2:627978/150‑1 (MQ=255)
                                                                                                                   cATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTTTCTGGGAGAT‑GCCGCG‑AAAATTTTTTCCTTTTGGGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTTCTCGCCTGGGCCATCCCTTTGGCGTGCAgtg                            >  2:661856/1‑150 (MQ=255)
                                                                                                                        cGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt                      <  1:455448/151‑1 (MQ=255)
                                                                                                                         gtgttGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgtg                     <  2:115901/148‑1 (MQ=255)
                                                                                                                                       gTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAc       >  1:341807/1‑151 (MQ=255)
                                                                                                                                            aTTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt  >  1:754902/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                          
CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  NC_000913/2095950‑2096240

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