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

CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  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)
           gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGcc                                                                                                                                       >  2:616417/1‑150 (MQ=255)
             tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTTGAAAATGcccc                                                                                                                                     >  2:93128/1‑148 (MQ=255)
               aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGATGCcgcg                                                                                                                                    >  1:407850/1‑149 (MQ=255)
                agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGaaa                                                                                                                                <  1:31194/151‑1 (MQ=255)
                agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGaaa                                                                                                                                <  2:334715/151‑1 (MQ=255)
                       gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCttt                                                                                                                         <  1:717413/151‑1 (MQ=255)
                             gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTg                                                                                                                   <  1:48370/151‑1 (MQ=255)
                                   gggtTCTACGATATCACCGCTTCCATCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCgg                                                                                                            <  2:168053/149‑1 (MQ=255)
                                       tCTACGATATCACTGCTTACAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACGTTCCGTTTTGTGTGGCGGGTAGGACTTAGAATCTATGAGATTCCGCGAAAAATTTTTCCTTTGGAGCggg                                                                                                           >  2:623934/1‑147 (MQ=255)
                                       tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATa                                                                                                         <  1:223503/150‑1 (MQ=255)
                                       tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATa                                                                                                         <  1:466538/150‑1 (MQ=255)
                                          aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGATGCCGCGGAAATTTTTTCCTTTTGCGCGGGTaaaa                                                                                                      >  1:207499/1‑150 (MQ=255)
                                             atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTc                                                                                                  <  2:197312/151‑1 (MQ=255)
                                              tatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCCCGGAAATTTTTTTCTTTGGTGCGGGTAAAAttttc                                                                                                 >  2:420649/1‑149 (MQ=255)
                                                  aCCGCTTCCAGCGATGGCCGTGCGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGTGGGTTCAGTCCCTCGAAGTTGGCCACCAACTTTTTTTCTTTTTGGCGGGTACAATTTCCATc                                                                                              >  2:130088/1‑150 (MQ=255)
                                                           aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTTCATCAACTCTcc                                                                                      >  2:76644/1‑149 (MQ=255)
                                                           aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGATGCCGCGGAAATTTTTTTCTTTGGCGCGGGTAAAATTTCCCTCAACTCTcc                                                                                      >  1:313977/1‑149 (MQ=255)
                                                                    cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGGTCTGGGAGATGCCGCGAAAATTTTTTTCTTTGGGGCGGGTAAAATTTTCATCAACTCTCCTGCGCTgg                                                                              >  2:98932/1‑148 (MQ=255)
                                                                           tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTTAGTCTATGGGAGATGCCGCGGAAATTCTTTTCTTTGGCGCGGGACAAAATTTCATCAACACTCCCTCGCTGGCACACCCa                                                                     >  2:514691/1‑149 (MQ=255)
                                                                           tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGGCTCTGGGAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCaa                                                                    >  2:212412/1‑150 (MQ=255)
                                                                             gATAAAAGCTGGGTATCTCGCGTGTCGGAAGTGATCGACATTCCGTTTTGTGGGGCGGGGGGGGCTAGGCCTCTGGGAGATGCCCCGAAAATTTTTTGTTTTGGGGCGGGCAAAATTTTCCTCAAACCTCCTTGGCTTGCAGAccccaga                                                                  >  2:68061/1‑146 (MQ=255)
                                                                                  aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACattaatt                                                            <  1:748791/151‑1 (MQ=255)
                                                                                       gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTc                                                        <  2:747554/150‑1 (MQ=255)
                                                                                       gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCg                                                       <  1:461817/151‑1 (MQ=255)
                                                                                            tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTTCCTCAACTCTCCTGCGCTGGCCGACCCAACATTAATTACTCGCCCgg                                                  >  2:203518/1‑151 (MQ=255)
                                                                                               cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGcc                                                <  1:723992/150‑1 (MQ=255)
                                                                                                gcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGa                                              <  1:231543/151‑1 (MQ=255)
                                                                                                  gTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGACTCTGGGAGAAGCCGCGAAAAATTTTTTCTTTTGGGCGGGAAAAAATTTTCTCCACTCCCCTGCCCTGGGGGAACCCCACTTAATTATTCTCCTCGcccagc                                            >  2:425809/1‑147 (MQ=255)
                                                                                                             gATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGt                                  <  1:612460/150‑1 (MQ=255)
                                                                                                              aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGGCCCTTGAAGATGGCCCGACAATTTTTTCCTTTGGCGGGGATAAAAATTCCCTCAACTCTCCTGCGCTGGCAGACCCCACATTAATTTTCCCCCTGGCCCATCGCTTTg                                      >  2:502857/1‑145 (MQ=255)
                                                                                                               tCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTATCTGGGAGATGCCGCGGAAATTTTTTTCTTTTGCGCGGGTAAAATTTCCCTCAACTCTCCTGCGCTGGCAGAACCCACATTTATTACTCCCCTGGCCGATCGCTTTTGCGTg                                 >  2:70789/1‑149 (MQ=255)
                                                                                                                cGACATTCCGTTTTGTTTTGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTTGCCGATCGCTTTGGCGTGCa                               <  2:627978/150‑1 (MQ=255)
                                                                                                                   cATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTTTCTGGGAGATGCCGCGAAAATTTTTTCCTTTTGGGCGGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTTCTCGCCTGGGCCATCCCTTTGGCGTGCAgtg                            >  2:661856/1‑150 (MQ=255)
                                                                                                                        cGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt                      <  1:455448/151‑1 (MQ=255)
                                                                                                                         gtgttGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgtg                     <  2:115901/148‑1 (MQ=255)
                                                                                                                                       gTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAc       >  1:341807/1‑151 (MQ=255)
                                                                                                                                            aTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt  >  1:754902/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                          
CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT  >  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: