Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,797,391	A→C	16.0%	L185R (CTG→CGG)	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,391	0	A	C	16.0%	186.2 / 28.0	76	L185R (CTG→CGG)	pheT	phenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base C (5/7);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTT  >  NC_000913/1797294‑1797490                                                                                                  |                                                                                                    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:1776571/100‑4 (MQ=255)           gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                         <  1:1828388/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:631908/100‑1 (MQ=255)                   ttCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACg                                                                                 <  1:944772/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                                                            >  1:483983/1‑100 (MQ=255)                                        cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                            <  2:1488790/100‑1 (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:2169906/2‑100 (MQ=255)                                                                                     ctGAACCAGCTGCCGCTGGTTCAGCACGACAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg               >  2:1498314/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)                                                                                                  aGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGtt  <  1:451412/100‑1 (MQ=255)                                                                                                  |                                                                                                    AAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTT  >  NC_000913/1797294‑1797490

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