Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,362,600	+C	noncoding (27/76 nt)	pheU ←	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,362,596	1	.	C	98.2%	213.0 / ‑4.3	56	noncoding (31/76 nt)	pheU	tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (17/38);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAAT‑CCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGTTTCTGTCAGGT  >  NC_000913/4362459‑4362738                                                                                                                                           |                                                                                                                                                aGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTAc                                                                                                                                    <  2:442850/151‑1 (MQ=255)     cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGAc                                                                                                                                >  2:431533/1‑151 (MQ=255)       gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTg                                                                                                                              <  1:219951/151‑1 (MQ=255)           aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCt                                                                                                                          <  1:407509/151‑1 (MQ=255)            aTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTa                                                                                                                         <  1:22772/151‑1 (MQ=255)              tGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATc                                                                                                                       >  1:262510/1‑151 (MQ=255)                 tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCgg                                                                                                                    >  2:111111/1‑151 (MQ=255)                  ttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCgg                                                                                                                    <  1:264436/150‑1 (MQ=255)                  ttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCgg                                                                                                                    <  1:895193/150‑1 (MQ=255)                   tGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCggg                                                                                                                   <  2:545113/150‑1 (MQ=255)                   tGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCggg                                                                                                                   <  1:61014/150‑1 (MQ=255)                       ggCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGTAGGTCCGTTTATATGAATTTGATGCCCGGACTCGGAATCCAACCAAGTACACGGGGACTTTCAAT‑CCCACTGCTCTACCGACTGAGCTATCCGGGCAAc                                                                                                               <  2:225306/150‑1 (MQ=25)                       ggCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACg                                                                                                              <  2:758980/151‑1 (MQ=255)                            aaCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGgcg                                                                                                         <  1:310091/151‑1 (MQ=255)                                cTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                                                                                     >  1:99782/1‑151 (MQ=255)                                   aaaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                                                                  <  1:159471/151‑1 (MQ=255)                                     aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAcc                                                                                                <  1:22440/151‑1 (MQ=255)                                      tCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCg                                                                                               >  2:341780/1‑151 (MQ=255)                                          ggCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTaa                                                                                            <  2:56000/150‑1 (MQ=255)                                             gTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTg                                                                                        >  1:147369/1‑151 (MQ=255)                                                  gAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAt                                                                                   >  2:401319/1‑151 (MQ=255)                                                  gAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAt                                                                                   <  2:822566/151‑1 (MQ=255)                                                           aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc                                                                          <  2:691099/151‑1 (MQ=255)                                                                aCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAt                                                                     >  1:377504/1‑151 (MQ=255)                                                                aCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAt                                                                     >  1:349585/1‑151 (MQ=255)                                                                       tGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGt                                                               >  1:539623/1‑150 (MQ=255)                                                                        ggAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaa                                                             <  1:71720/151‑1 (MQ=255)                                                                          aGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaaaa                                                           <  1:358355/151‑1 (MQ=255)                                                                          aGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaaaa                                                           <  2:329932/151‑1 (MQ=255)                                                                              ccGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAgag                                                        <  1:694188/150‑3 (MQ=255)                                                                                   tataTGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCa                                                   <  2:389433/150‑1 (MQ=255)                                                                                   tataTGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCa                                                   <  2:374800/150‑1 (MQ=255)                                                                                      atGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTg                                               <  1:154475/151‑1 (MQ=255)                                                                                        gAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGc                                              <  2:101247/150‑1 (MQ=255)                                                                                        gAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCt                                             <  2:262478/151‑1 (MQ=255)                                                                                          aTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTa                                           <  1:507547/151‑1 (MQ=255)                                                                                              ggTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGAt                                       <  1:494878/151‑1 (MQ=255)                                                                                                tGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTg                                     >  1:501453/1‑151 (MQ=255)                                                                                                 gCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGc                                    <  1:754815/151‑1 (MQ=255)                                                                                                  cccGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                                   <  1:530011/151‑1 (MQ=255)                                                                                                      gACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCa                                >  1:919301/1‑150 (MQ=255)                                                                                                          cGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAAt                             <  1:501946/149‑1 (MQ=255)                                                                                                          cGGAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATc                            <  1:373419/150‑1 (MQ=255)                                                                                                           ggAATCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGc                          >  2:225783/1‑151 (MQ=255)                                                                                                             aaTCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGctt                        >  1:210603/1‑151 (MQ=255)                                                                                                             aaTCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGctt                        <  1:713609/151‑1 (MQ=255)                                                                                                              aTCGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGcttc                       <  1:166408/151‑1 (MQ=255)                                                                                                                cGAACCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGcttctt                     >  2:387755/1‑151 (MQ=255)                                                                                                                  aaCCAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCtttt                   >  2:847454/1‑151 (MQ=255)                                                                                                                    ccAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt                 <  1:705029/151‑1 (MQ=255)                                                                                                                     cAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg                <  1:354442/151‑1 (MQ=255)                                                                                                                     cAAGGACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg                >  2:297045/1‑151 (MQ=255)                                                                                                                        ggACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTgtgt             <  1:91550/151‑1 (MQ=255)                                                                                                                         gACACGGGGATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGtt            <  1:894516/151‑1 (MQ=255)                                                                                                                                ggATTTTCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGTTTCTGTCa     <  2:673881/151‑1 (MQ=255)                                                                                                                                    tttCAATCCCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGTTTCTGTCAGGt  <  1:298603/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                AGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAAT‑CCC‑CTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTGTGTTTCTGTCAGGT  >  NC_000913/4362459‑4362738

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