Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,797,388	G→T	15.6%	P186Q (CCG→CAG)	pheT ←	phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,797,388	0	G	T	15.6%	206.3 / 21.5	75	P186Q (CCG→CAG)	pheT	phenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base T (5/7);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAA  >  NC_000913/1797290‑1797487                                                                                                   |                                                                                                    ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGCACCAGCgg                                                                                                    >  2:689684/1‑100 (MQ=255) ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCgg                                                                                                    >  2:29823/1‑100 (MQ=255) ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCgg                                                                                                    >  1:1843079/1‑100 (MQ=255)     aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGc                                                                                                <  1:848797/100‑1 (MQ=255)       gATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                              <  1:60193/100‑1 (MQ=255)         tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGCTGAACCAGCTGCCGCTGGt                                                                                            <  1:45874/100‑1 (MQ=255)          aaCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCCCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt                                                                                           >  1:1974086/1‑100 (MQ=255)           aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt                                                                                           <  1:2129330/99‑1 (MQ=255)           aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                          >  2:1460432/1‑100 (MQ=255)           aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                          >  1:337851/1‑100 (MQ=255)            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCa                                                                                         <  2:642993/100‑1 (MQ=255)            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGtt                                                                                           >  2:1076206/1‑98 (MQ=255)            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                         >  1:375819/1‑100 (MQ=255)             gcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAg                                                                                        <  1:536668/100‑1 (MQ=255)               gggCCGGCTGCCGGGGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAGCGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                          <  2:2300614/96‑1 (MQ=255)               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  1:1828388/100‑4 (MQ=255)               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  1:1776571/100‑4 (MQ=255)               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCa                                                                                      <  1:816041/100‑1 (MQ=255)                 gCAGGCTTCCGGCGCTTCGTCTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCCCTGGTTCaacc                                                                                      <  2:540194/98‑4 (MQ=255)                  cAGGCTTCCGGCGCTTCGACTGTAATCGGCCGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACgg                                                                                   >  1:2061810/1‑100 (MQ=255)                  cAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaaccgg                                                                                    <  2:611820/99‑6 (MQ=39)                  cAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACgg                                                                                   >  2:698813/1‑100 (MQ=255)                   aGGCTTCCGGCGCCTCGCCTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGc                                                                                  <  2:923639/100‑1 (MQ=255)                    ggCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCa                                                                                 >  2:1433539/1‑100 (MQ=255)                     gCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCaa                                                                                >  1:1531573/1‑100 (MQ=255)                     gCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGCTGAACCAGCTGCCGCTGGTTCaaccgga                                                                                   <  2:2185420/97‑7 (MQ=38)                       ttCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACg                                                                              <  1:944772/100‑1 (MQ=255)                       ttCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACg                                                                              <  1:631908/100‑1 (MQ=255)                         ccGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTc                                                                            >  1:859534/1‑100 (MQ=255)                          cGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCa                                                                           >  2:428461/1‑100 (MQ=255)                              gcTTCGCCTGTAATCGGCAGCGTGTCTTCGATGGTCGCACCAACCGGAACGATTTCCTGTTGAACCAGCGGCTGCTGGTTCAGCACGGCAACGTCAcg                                                                         <  2:1514442/98‑1 (MQ=255)                               cTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                                      <  2:359598/100‑1 (MQ=255)                                ttCGACTTTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGCTTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCa                                                                     <  2:1957284/100‑1 (MQ=255)                                  cGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAac                                                                   <  2:216647/100‑1 (MQ=255)                                     cTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGGTCAGCACGGCAACGTCACGCGCAacaca                                                                <  2:1754368/100‑2 (MQ=255)                                           tCGGCAGCCTGTCGACGAAGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCGGGGTGATTACAGCAACGTCACGCGCAACACCa                                                               >  2:798939/1‑95 (MQ=255)                                            ctgCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         <  2:1855977/98‑1 (MQ=255)                                            cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         <  2:1488790/100‑1 (MQ=255)                                            cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         >  1:483983/1‑100 (MQ=255)                                             ggCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                        >  2:1469715/1‑100 (MQ=255)                                             ggCAGCGTGTCGTCGATGGTCGCACCAACCAGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                        >  1:992134/1‑100 (MQ=255)                                               cAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         <  2:1591879/97‑1 (MQ=255)                                                    tgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCa                                                 >  1:1038887/1‑100 (MQ=255)                                                          cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGc                                           <  1:800881/100‑1 (MQ=255)                                                                tCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCATTCGGCACGGtt                                     <  2:1245544/100‑1 (MQ=255)                                                                tCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGtt                                     >  1:920504/1‑100 (MQ=255)                                                                 cGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                    <  1:417836/100‑1 (MQ=255)                                                                 cGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                    >  1:1923109/1‑100 (MQ=255)                                                                  gCACTAACCGGCACGAGTTCCGGTGCAACCTGCGGCAGCTGGTTCAGAACGGCAGCGTCAGGCTCAACAGCGATGATGCCTAAGCAGTCGGCACGGTTTg                                   <  2:218992/100‑1 (MQ=255)                                                                  gCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTg                                   <  2:1367664/100‑1 (MQ=255)                                                                  gCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCGGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTg                                   <  2:564208/100‑1 (MQ=255)                                                                    aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGc                                 >  1:606013/1‑100 (MQ=255)                                                                    aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGc                                 >  1:605954/1‑100 (MQ=255)                                                                    aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGc                                 >  1:1482468/1‑100 (MQ=255)                                                                     ccAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                <  1:1834150/100‑1 (MQ=255)                                                                     ccAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAACGATGCCTAAGCAGTCGGCACGGTTTGGCg                                <  1:1957424/100‑1 (MQ=255)                                                                     ccAACCGGAACAATGTCCGGTTGACCCTGTGGCAGCTAGCTCCGCACGCCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCCGCACGGtttcg                                  <  2:831688/98‑3 (MQ=255)                                                                         ccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCAc                            >  2:970923/1‑100 (MQ=255)                                                                          cGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACg                           >  1:1526453/1‑100 (MQ=255)                                                                           ggAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGCTTGGCGTCACGc                          <  1:1568212/100‑1 (MQ=255)                                                                               cGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt                      >  2:1540817/1‑100 (MQ=255)                                                                               cGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt                      <  1:1767394/100‑1 (MQ=255)                                                                                gATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAtt                     >  1:212256/1‑100 (MQ=255)                                                                                    tCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGa                 >  2:932683/1‑100 (MQ=255)                                                                                    ccgtgctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGa                 >  1:1538626/7‑100 (MQ=38)                                                                                     ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                <  2:1716482/100‑1 (MQ=255)                                                                                     ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                <  1:797117/100‑1 (MQ=255)                                                                                      cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg               <  1:1433547/100‑1 (MQ=255)                                                                                      cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg               <  2:1843065/100‑1 (MQ=255)                                                                                      cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg               <  1:1743437/100‑1 (MQ=255)                                                                                       ggTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg              >  1:2170281/1‑100 (MQ=255)                                                                                        gctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgt             >  2:2377041/3‑100 (MQ=255)                                                                                         ttGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            >  1:618543/1‑100 (MQ=255)                                                                                         ctGAACCAGCTGCCGCTGGTTCAGCACGACAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            >  2:1498314/2‑100 (MQ=255)                                                                                         ctGAACCAGCTGCCGCTGGTTCAGCACGACAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            >  2:2169906/2‑100 (MQ=255)                                                                                           gAACAAGCTGCCGCTGGTTAAGCACGGCAACATCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgtt          >  2:1151954/1‑100 (MQ=255)                                                                                              ccAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtc       >  1:863776/1‑100 (MQ=255)                                                                                                 gCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACGCCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg         <  2:2183109/95‑1 (MQ=255)                                                                                                   ggCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCaa  >  2:1132397/1‑100 (MQ=255)                                                                                                   |                                                                                                    GGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAA  >  NC_000913/1797290‑1797487

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