Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	491,743	0	T	G	24.7%	140.4 / 62.9	81	I111S (ATC→AGC)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base G (9/11);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

CGTTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCG  >  NC_000913/491648‑491842                                                                                                |                                                                                                    cGTTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCAc                                                                                                 >  2:1236874/1‑100 (MQ=255)  gTTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACg                                                                                                <  2:1214262/100‑1 (MQ=255)     ggCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGtt                                                                                              <  1:1009969/99‑1 (MQ=255)     ggCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTg                                                                                             >  2:1561520/1‑100 (MQ=255)      gCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGa                                                                                            >  2:1164959/1‑100 (MQ=255)        tttGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGTAACTTACGTCCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGa                                                                                            <  1:1699031/98‑1 (MQ=255)        tttGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATg                                                                                          <  2:127511/100‑1 (MQ=255)        tttGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATg                                                                                          <  1:1603390/100‑1 (MQ=255)         ttGTACCGGTCCCTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGc                                                                                         >  2:1279968/1‑100 (MQ=255)           gTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCa                                                                                       >  1:836044/1‑100 (MQ=255)            tACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCAt                                                                                      >  2:857776/1‑100 (MQ=255)            tACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCAt                                                                                      >  1:1935372/1‑100 (MQ=255)             aCCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATc                                                                                     <  2:1084393/100‑1 (MQ=255)             aCCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATc                                                                                     <  2:78585/100‑1 (MQ=255)             aCCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATc                                                                                     >  2:62661/1‑100 (MQ=255)             aCCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATc                                                                                     <  2:391362/100‑1 (MQ=255)             aCCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATc                                                                                     <  1:1494131/100‑1 (MQ=255)                ggTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCaaa                                                                                  >  2:962601/1‑100 (MQ=255)                ggTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCaaa                                                                                  >  1:681512/1‑100 (MQ=255)                ggTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCaaa                                                                                  >  2:1796253/1‑100 (MQ=255)                 gTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTCACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAAc                                                                                 >  1:1279215/1‑100 (MQ=255)                   ccGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCg                                                                               >  2:1894433/1‑100 (MQ=255)                    cGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCgg                                                                              >  1:1355024/1‑100 (MQ=255)                     gTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCggg                                                                             >  2:455262/1‑100 (MQ=255)                     gTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCggg                                                                             <  1:1788095/100‑1 (MQ=255)                       aaaCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGc                                                                            <  2:872973/99‑1 (MQ=255)                       aaaCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCg                                                                           >  2:311878/1‑100 (MQ=255)                         aCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGAc                                                                         <  1:782675/100‑1 (MQ=255)                         aCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGAc                                                                         <  2:279986/100‑1 (MQ=255)                         aCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGAc                                                                         >  2:1653865/1‑100 (MQ=255)                           cGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACaa                                                                       >  1:1905846/1‑100 (MQ=255)                               cAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGtt                                                                   >  1:1696704/1‑100 (MQ=255)                               cAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGtt                                                                   >  1:588837/1‑100 (MQ=255)                               cAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGt                                                                    >  1:1051956/1‑99 (MQ=255)                                   cTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCtgg                                                               <  1:1614723/100‑1 (MQ=255)                                     gCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCtggtg                                                             >  1:1775917/1‑100 (MQ=255)                                     gCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCtggtg                                                             >  2:1067638/1‑100 (MQ=255)                                        gcgTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCtggtggtg                                                          <  1:1796242/100‑1 (MQ=255)                                         cgTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCtggtggtgg                                                         <  2:823074/100‑1 (MQ=255)                                          gTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGga                                                        >  1:418711/1‑100 (MQ=255)                                           tGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGgac                                                       >  1:1719372/1‑100 (MQ=255)                                           tGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGgac                                                       >  1:1719276/1‑100 (MQ=255)                                            gAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGgacg                                                      <  1:112946/100‑1 (MQ=255)                                              aaCCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGgacgac                                                    >  1:1477572/1‑100 (MQ=255)                                                ccATCAGTGAACATTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGgacgacn                                                   <  1:93/99‑2 (MQ=255)                                                ccATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACct                                                  >  2:889138/1‑100 (MQ=255)                                                ccATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACct                                                  >  2:888950/1‑100 (MQ=255)                                                     aGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTgg                                             >  2:1835717/1‑100 (MQ=255)                                                     aGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTgg                                             >  1:5261/1‑100 (MQ=255)                                                     aGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTgg                                             <  2:1117017/100‑1 (MQ=255)                                                      gTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGc                                            <  1:1835705/100‑1 (MQ=255)                                                      gTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGc                                            >  2:1497644/1‑100 (MQ=255)                                                       tGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCa                                           <  1:973969/100‑1 (MQ=255)                                                        gAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCaa                                          <  1:1828643/100‑1 (MQ=255)                                                          aaCTTACGACCTGGAATACGGCACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAAcc                                        <  1:199740/100‑1 (MQ=255)                                                           aCTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACcg                                       <  1:780308/100‑1 (MQ=255)                                                               aCGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACcggcgg                                   <  1:1614481/100‑1 (MQ=255)                                                                   ccTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCAc                                <  1:1221922/99‑1 (MQ=255)                                                                   ccTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACt                               >  1:594872/1‑100 (MQ=255)                                                                      ggAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATc                            >  2:1113691/1‑100 (MQ=255)                                                                       gAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCg                           <  1:1561509/100‑1 (MQ=255)                                                                       gAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGCTGTACGACCTGCTGGCAACCGGCGGCACTATCg                           <  1:1279961/100‑1 (MQ=255)                                                                           aCGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGc                       <  2:1399631/100‑1 (MQ=255)                                                                            cGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGGAGCg                      >  1:758952/1‑100 (MQ=255)                                                                            cGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCg                      >  2:81944/1‑100 (MQ=255)                                                                               cACCGATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGAcc                   <  2:1905859/100‑1 (MQ=255)                                                                                 ccGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGt                 <  2:809006/100‑1 (MQ=255)                                                                                   gATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTa               >  1:1285951/1‑100 (MQ=255)                                                                                   gATCAGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTa               <  1:1653855/100‑1 (MQ=255)                                                                                      cAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAAc            <  1:864763/100‑1 (MQ=255)                                                                                       aGCTGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACt           >  2:450201/1‑100 (MQ=255)                                                                                        gCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTg          <  1:937747/100‑1 (MQ=255)                                                                                         cTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGa         <  2:129536/100‑1 (MQ=255)                                                                                          tGGAGAGCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGAt        >  2:1733054/1‑100 (MQ=255)                                                                                           ggAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGaag       >  1:129535/1‑98 (MQ=255)                                                                                           ggAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATc       <  2:681513/100‑1 (MQ=255)                                                                                           ggAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATc       <  1:65411/100‑1 (MQ=255)                                                                                           ggAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATc       <  1:215366/100‑1 (MQ=255)                                                                                             agaTCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCcg     <  2:803027/100‑1 (MQ=255)                                                                                              gaTCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCcgt    >  1:1869845/1‑100 (MQ=255)                                                                                              gaTCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCcg     <  2:803039/99‑1 (MQ=255)                                                                                                tCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCcgtcg  >  1:763446/1‑100 (MQ=255)                                                                                                tCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCcgtcg  >  1:763340/1‑100 (MQ=255)                                                                                                |                                                                                                    CGTTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCG  >  NC_000913/491648‑491842

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