Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	632,661	0	G	T	27.9%	132.3 / 55.9	69	L273I (CTC→ATC)	ybdH	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base T (5/14);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GAAAACGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCAC  >  NC_000913/632562‑632759                                                                                                    |                                                                                                   gAAAACGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGgag                                                                                                    <  1:1635786/100‑1 (MQ=255)  aaaaCGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGgaga                                                                                                   >  1:246586/1‑100 (MQ=255)  aaaaCGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGgag                                                                                                    >  2:1454401/1‑99 (MQ=255)      cGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACt                                                                                               <  2:596003/100‑1 (MQ=255)      cGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACt                                                                                               <  2:1977597/100‑1 (MQ=255)      cGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACt                                                                                               <  2:1796462/100‑1 (MQ=255)       gCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACtt                                                                                              <  1:2335254/100‑1 (MQ=255)        cTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTtc                                                                                             >  2:871767/1‑100 (MQ=255)          naTACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTtctc                                                                                           >  1:2077633/2‑100 (MQ=255)          gATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTtctc                                                                                           <  2:377962/100‑1 (MQ=255)          gATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTtctc                                                                                           >  1:1872890/1‑100 (MQ=255)                tctcAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTg                                                                                     <  2:208335/97‑1 (MQ=255)                   cAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTgaga                                                                                  >  1:208337/1‑97 (MQ=255)                    aGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGgc                                                                                 <  1:2320041/100‑1 (MQ=255)                    aGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGgc                                                                                 <  1:1247004/100‑1 (MQ=255)                       tAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGgcag                                                                               <  2:1824341/99‑1 (MQ=255)                           tGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACgg                                                                          <  1:1783648/100‑1 (MQ=255)                            ncgcCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGt                                                                         >  2:2161074/2‑100 (MQ=255)                            gcgcCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGt                                                                         >  2:1550064/1‑100 (MQ=255)                            gcgcCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGt                                                                         >  2:1459652/1‑100 (MQ=255)                             cgcCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTc                                                                        <  2:133545/100‑1 (MQ=255)                              gcCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCa                                                                       <  1:871769/100‑1 (MQ=255)                               ccAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAg                                                                      >  1:2199206/1‑100 (MQ=255)                               ccAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAg                                                                      <  2:1289186/100‑1 (MQ=255)                                cAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAgn                                                                     <  1:536566/100‑2 (MQ=255)                                     catcatCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGtt                                                                <  1:439260/100‑1 (MQ=255)                                     catcatCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGtt                                                                <  1:439517/100‑1 (MQ=255)                                        catcCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATg                                                             <  1:1675083/100‑1 (MQ=255)                                        catcCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATg                                                             <  2:2245706/100‑1 (MQ=255)                                        catcCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATg                                                             >  1:328536/1‑100 (MQ=255)                                         atcCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGc                                                            <  2:28141/100‑1 (MQ=255)                                             tGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACgg                                                        <  1:1990331/100‑1 (MQ=255)                                             tGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACgg                                                        >  2:786622/1‑100 (MQ=255)                                             tGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACgg                                                        <  1:786623/100‑1 (MQ=255)                                              gACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGc                                                       >  2:39394/1‑100 (MQ=255)                                              gACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGc                                                       <  2:731757/100‑1 (MQ=255)                                              gACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGc                                                       >  2:384551/1‑100 (MQ=255)                                                cccAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCAt                                                     <  1:384553/100‑1 (MQ=255)                                                 ccAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATg                                                    <  1:548276/100‑1 (MQ=255)                                                 ccAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATg                                                    <  1:1459656/100‑1 (MQ=255)                                                  cAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGa                                                   >  2:34997/1‑100 (MQ=255)                                                   aGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAg                                                  <  2:29415/100‑1 (MQ=255)                                                    gCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGc                                                 <  1:1528645/100‑1 (MQ=255)                                                      aaGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTg                                               >  1:179197/1‑100 (MQ=255)                                                       aGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTgc                                              >  1:342635/1‑100 (MQ=255)                                                         gcgcTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTgccg                                            >  1:2001509/1‑100 (MQ=255)                                                           gcTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTgccgcc                                          <  1:1728052/100‑1 (MQ=255)                                                                 cACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACAcgc                                    >  1:637919/1‑100 (MQ=255)                                                                  aCCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACAcgcn                                   <  1:2161081/100‑2 (MQ=255)                                                                  aCCAGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACAcgcg                                   >  1:2264947/1‑100 (MQ=255)                                                                   ccAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGt                                  >  1:529424/1‑100 (MQ=255)                                                                     aGAATTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTaa                                >  1:371433/1‑100 (MQ=255)                                                                     aGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTaa                                >  1:1021459/1‑100 (MQ=255)                                                                      gAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTaaa                               <  2:1731352/100‑1 (MQ=255)                                                                      gAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATACACGGCATGAGCTGCCGCCACACGAGTa                                 >  2:2094708/1‑98 (MQ=255)                                                                        aTTCCGTAGGCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAAc                             <  2:1346792/100‑1 (MQ=255)                                                                         ttCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGTTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACg                            <  2:1273929/100‑1 (MQ=255)                                                                           ccGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGAt                          >  1:864760/1‑100 (MQ=255)                                                                             gTATGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCg                        <  2:2199198/100‑1 (MQ=255)                                                                                ggCGACTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGccc                     <  2:1414299/100‑1 (MQ=255)                                                                                ggCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGccc                     >  1:45049/1‑100 (MQ=255)                                                                                  cGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAg                   <  2:1143365/100‑1 (MQ=255)                                                                                     cTTTGGTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGacc                <  1:2143067/100‑1 (MQ=255)                                                                                     cTTTGGTGCCGGTGAAACACTTCTCGGTTTGCGTCAGCACGGACAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGacc                >  2:201428/1‑100 (MQ=255)                                                                                       ttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGaccac              <  1:1075911/100‑1 (MQ=255)                                                                                         ggTGCCGTGGATAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGAccacca            <  2:2144343/100‑1 (MQ=255)                                                                                          gTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCaccaa           >  2:2085714/1‑100 (MQ=255)                                                                                            gCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCaccaacc         <  2:342633/100‑1 (MQ=255)                                                                                              cGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCAt       >  1:370225/1‑100 (MQ=255)                                                                                                 ggAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATccc    <  2:2131342/100‑1 (MQ=255)                                                                                                 ggAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATccc    >  2:1374348/1‑100 (MQ=255)                                                                                                 ggAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATccc    <  1:1550068/100‑1 (MQ=255)                                                                                                   agaAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCccac  <  2:179195/100‑1 (MQ=255)                                                                                                    |                                                                                                   GAAAACGCTGATACGCTCCAGTTAATTGCGCCAGCACATCATCCTGACCCAGCAAGGCGCTTTGCACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCAC  >  NC_000913/632562‑632759

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