Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,096,090	0	A	G	43.6%	71.9 / 15.3	56	I83V (ATT→GTT)	hisF	imidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (5/26);  new base G (24/0);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 7.81e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTA  >  NC_000913/2095946‑2096228                                                                                                                                                 |                                                                                                                                           cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                                                                        >  2:625165/1‑149 (MQ=255)  tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                                                                        >  1:46358/1‑148 (MQ=255)   ggCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtct                                                                                                                                    <  2:517496/151‑1 (MQ=255)   ggCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGt                                                                                                                                      <  2:1016532/149‑1 (MQ=255)     cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTctct                                                                                                                                  <  2:985701/151‑1 (MQ=255)      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGtct                                                                                                                                    >  2:53046/1‑148 (MQ=255)        aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTCCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGTTTAGGCTCCgg                                                                                                                                >  2:312471/1‑150 (MQ=255)        aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGAGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGttgaagct                                                                                                                                    >  2:372501/1‑140 (MQ=255)         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGaa                                                                                                                              <  1:597714/151‑1 (MQ=255)            ctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTAAGTATCTgg                                                                                                                                >  1:393157/1‑146 (MQ=255)             tacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTg                                                                                                                                 >  2:1011325/1‑144 (MQ=255)             tacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGaa                                                                                                                              >  2:241835/1‑147 (MQ=255)              acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTCTCTGGAAGaggg                                                                                                                         >  1:619209/1‑148 (MQ=255)                gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGAGGGGCTAAGAATCTGGAAGagaaac                                                                                                                       >  1:176831/1‑146 (MQ=255)                  tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGATGCcgc                                                                                                                      >  2:955387/1‑150 (MQ=255)                  tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGAAGCcgcg                                                                                                                     >  1:154558/1‑151 (MQ=255)                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGa                                                                                                                    <  1:60438/151‑1 (MQ=255)                    aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGaa                                                                                                                   <  2:365384/151‑1 (MQ=255)                         ggCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTc                                                                                                              <  2:316166/151‑1 (MQ=255)                               gACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTc                                                                                                          <  1:530596/149‑1 (MQ=255)                               gACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTCTCTGGGAGATGCCGCGAAAAATtttttc                                                                                                         >  2:920308/1‑148 (MQ=255)                                       ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcgc                                                                                                 <  2:932457/150‑1 (MQ=255)                                       ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcg                                                                                                  <  1:619709/149‑1 (MQ=255)                                          gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGAt                                                                                             <  1:616844/151‑1 (MQ=255)                                            tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATa                                                                                            <  1:600540/150‑1 (MQ=255)                                             cTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTGTCTGGAAGATGACGCGAAAAATCTTTCCTTTGGCGCgg                                                                                               >  2:256080/1‑146 (MQ=255)                                              tACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAAGCTCTGGAAGAGGTCGCGTAAAATCTTTTCTTTTGCGCGGGTaaa                                                                                          >  1:651616/1‑150 (MQ=255)                                               aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATaa                                                                                           <  2:368067/148‑1 (MQ=255)                                                cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAAt                                                                                        <  2:630028/150‑1 (MQ=255)                                                 gATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAAATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAAttt                                                                                      >  1:70562/1‑151 (MQ=255)                                                    atCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTcc                                                                                    <  1:917200/150‑1 (MQ=255)                                                     tCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGGTAAAAttt                                                                                      >  2:102389/1‑147 (MQ=255)                                                        ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                                                                               <  2:210975/151‑1 (MQ=255)                                                        ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                                                                               <  1:670968/151‑1 (MQ=255)                                                          gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTGGCGCGGGTAAAATTTCCCTCAAc                                                                              >  1:51769/1‑150 (MQ=255)                                                              ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACtctc                                                                          <  1:955111/150‑1 (MQ=255)                                                                 gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTg                                                                       <  1:853231/150‑1 (MQ=255)                                                                 gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATAAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgc                                                                      <  1:1011391/151‑1 (MQ=255)                                                                 gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTTTTTCCTTTTGCGCGGGTAAAATTTTCATCAACTCTcc                                                                         >  2:515639/1‑148 (MQ=255)                                                                         cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAg                                                              <  2:145423/151‑1 (MQ=255)                                                                                    tAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGCTTAATACTATGGCAAAAGCCGCGGAAATTATTTTCTTTCGAGGGGAAAAAATATCCATCAACACTTCAGCCCGGGCacagcccacctt                                                    >  2:175427/1‑139 (MQ=255)                                                                                       aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGGCTCTGGGAGAAGGCGCGGAAATTTTTTTCTTTTGGGCGGGTAAAAATTTCATCAACTCTCCCTCGCCGGCAGACCCAACattaa                                                  >  1:207902/1‑149 (MQ=255)                                                                                       aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCTGGAAGATGCCCCGGAAATTTTTTTCTTTGGCGCGGGTAAAATTTTCATCAACTCCCCCGCGCTGGCAGACCCAACaataa                                                  >  1:159092/1‑149 (MQ=255)                                                                                           tGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGGAAATTCTTTCCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCCGGCAGACCCAACattaatta                                               >  1:622497/1‑148 (MQ=255)                                                                                           tGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGGAAAATTTTTTCTTTGGCGCGGGTAAAATTTTCATCAAACCTCCTGCGCTGGGAGACCCCACattaatt                                                >  1:333577/1‑147 (MQ=255)                                                                                                aTCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCt                                        <  1:547743/150‑1 (MQ=255)                                                                                                  ctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAAATCTCTGGGAAAAGCCGCGCGAAATCTTTTCTTTGGGGGGGGAAAAAATTCCCTGAACTCTCCTGCGCTGGCCGACACAACATTTATTACtttcc                                         >  1:35204/1‑143 (MQ=255)                                                                                                      cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCTTGGAGAAGCCGCGAAAAATTTTTTCTTTTGCGGGGGTAAAAATTTCATCAAAACTCCTCCGCTGGGCAACCCAAAATTAATTACTCCCCCgg                                      >  1:90553/1‑146 (MQ=255)                                                                                                            ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGGAGATGCCCCGAAAAATCTTTTCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTTCGCTGGCAGACCCAACATTAATTACTCGCCTTGCCGATCGCttt                           >  2:227936/1‑151 (MQ=255)                                                                                                             gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGGAAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCttt                           >  1:845965/1‑150 (MQ=255)                                                                                                             gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTg                          <  2:688666/151‑1 (MQ=255)                                                                                                             gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTg                          <  1:696354/151‑1 (MQ=255)                                                                                                                   aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTTTCTTGAAGATGCCGCGAAAATTCTTTTCTTTTGCGCGGATAAAATTTCCATCAACACTCCTGCGCCGGCAGACCCAACCATAATTACTCTCCCGGGCggacgcgtttgggtg                     >  2:103567/1‑136 (MQ=255)                                                                                                                     cGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAg                  <  1:717787/151‑1 (MQ=255)                                                                                                                       aCATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTGGAAGATGCCCCGGAAATTTTTTTCTTTGGGGCGGGTAAAAATTCCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTACTCTCCTGGCCGATCGCTTTTGCGTGGAgtg                >  1:178324/1‑151 (MQ=255)                                                                                                                                     gtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTa  <  1:596718/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                           CTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTA  >  NC_000913/2095946‑2096228

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