Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,034,140	C→A	7.7%	L345L (CTG→CTT)	lysS ←	lysine tRNA synthetase, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,034,140	0	C	A	7.7%	157.5 / 3.4	52	L345L (CTG→CTT)	lysS	lysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base A (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGC  >  NC_000913/3033993‑3034287                                                                                                                                                    |                                                                                                                                                    cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTcc                                                                                                                                                     <  7:245928/148‑1 (MQ=255)          aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGcc                                                                                                                                            <  8:26271/148‑1 (MQ=255)          aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGcc                                                                                                                                            <  1:35365/148‑1 (MQ=255)                     tCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTc                                                                                                                                 >  4:144947/1‑148 (MQ=255)                         cTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGGTTCCggg                                                                                                                             >  6:17814/1‑148 (MQ=255)                            aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGa                                                                                                                          >  2:142417/1‑148 (MQ=255)                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGAt                                                                                                               <  3:262574/148‑1 (MQ=255)                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGAt                                                                                                               <  7:330297/148‑1 (MQ=255)                                               aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCAc                                                                                                       <  1:142417/148‑1 (MQ=255)                                                  tcttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACg                                                                                                      >  6:15378/1‑146 (MQ=255)                                                   cttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGc                                                                                                     <  5:15378/146‑1 (MQ=255)                                                        aaGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACTTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCAAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGt                                                                                                >  2:295828/1‑146 (MQ=255)                                                        aaGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCAAGATCCGCCATGTCGGGTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGt                                                                                                >  2:291377/1‑146 (MQ=255)                                                         aGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAg                                                                                             >  4:10087/1‑148 (MQ=255)                                                          gATCTGGGTTACGATACGGCCCAGACCCCAGCGCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGTAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGAGCGCTTCACGCTTTGTCAg                                                                                             <  4:203428/147‑1 (MQ=255)                                                            tctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGttt                                                                                          >  1:279565/1‑148 (MQ=255)                                                                ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGCCAGttttt                                                                                        <  8:951/146‑1 (MQ=255)                                                                     cGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACgg                                                                                 <  3:144947/148‑1 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGatat                                                                                                                       >  3:295398/1‑94 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGatat                                                                                                                       <  4:295398/94‑1 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGc                                                                                                            >  6:236417/1‑105 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGc                                                                                                            <  5:236417/105‑1 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACGGAAGTCCAGc                                                                 >  8:244023/1‑148 (MQ=255)                                                                                     ccAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGc                                                                 >  6:188171/1‑148 (MQ=255)                                                                                       aGCTCTTCTCAACGTGGATGCCGATAGATTCAGAAATTGCTTTCGCAGAGGCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGt                                                               <  7:69595/148‑1 (MQ=255)                                                                                        gCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGt                                                               >  7:254761/1‑147 (MQ=255)                                                                                         ctctTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCa                                                             <  8:60210/148‑1 (MQ=255)                                                                                             tctcAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTc                                                         >  1:319202/1‑148 (MQ=255)                                                                                              ctcAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCg                                                        >  6:138491/1‑148 (MQ=255)                                                                                                             gATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTcc                                         >  8:147696/1‑148 (MQ=255)                                                                                                              aTAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCg                                        >  4:100778/1‑148 (MQ=255)                                                                                                              aTAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGGCACGGCGCCGTAGGTCACTTCCg                                        <  5:17814/148‑1 (MQ=255)                                                                                                                 gATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCt                                     <  5:188171/148‑1 (MQ=255)                                                                                                                      aGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCg                                <  6:139376/148‑1 (MQ=255)                                                                                                                        cAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCAGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTAcc                                 >  7:318756/1‑145 (MQ=255)                                                                                                                          aTTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCg                                <  3:176790/144‑1 (MQ=255)                                                                                                                          aTTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCg                                >  4:176790/1‑144 (MQ=255)                                                                                                                            tGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAgg                                                 <  8:323051/125‑1 (MQ=255)                                                                                                                            tGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAgg                                                 >  7:323051/1‑125 (MQ=255)                                                                                                                              cTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATc                        <  5:138491/148‑1 (MQ=255)                                                                                                                                   gCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCtgtg                   <  8:254761/148‑1 (MQ=255)                                                                                                                                    cAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGc                  >  7:208775/1‑148 (MQ=255)                                                                                                                                     agagTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCAGCTGGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGc                  >  1:85587/1‑147 (MQ=255)                                                                                                                                       agTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCag               <  7:147696/148‑1 (MQ=255)                                                                                                                                          cGAAGTTGTCCAGATCAGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCagcg             >  1:301074/1‑145 (MQ=255)                                                                                                                                            aaGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGCTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTAc          >  5:266868/1‑148 (MQ=255)                                                                                                                                              gTTGTCAAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGAGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACgg        <  1:295828/148‑1 (MQ=255)                                                                                                                                              gTTGTCAAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGAGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACgg        <  1:291377/148‑1 (MQ=255)                                                                                                                                               ttGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACgg        >  4:213302/1‑147 (MQ=255)                                                                                                                                               ttGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGa       <  3:175118/148‑1 (MQ=255)                                                                                                                                                 gTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGaa      >  2:273083/1‑147 (MQ=255)                                                                                                                                                  tCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGaa      >  4:226718/1‑146 (MQ=255)                                                                                                                                                  tCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAAc     <  8:208775/147‑1 (MQ=255)                                                                                                                                                   ccAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAg   >  7:190933/1‑148 (MQ=255)                                                                                                                                                    cAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGc  >  7:139228/1‑148 (MQ=255)                                                                                                                                                    cAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGc  >  4:63177/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGAAGTTGTCCAGATCCGCCATGTCGGTTTCCGGGCGATATTTCTTGATCGCTTCACGCATGGTCAGTTTTTCGAACGGTTTACCGAAGTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGC  >  NC_000913/3033993‑3034287

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