Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	228,963	T→A	10.0%	noncoding (36/77 nt)	aspU →	tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	228,963	0	T	A	10.0%	122.3 / 5.5	40	noncoding (36/77 nt)	aspU	tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base A (2/2);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACA  >  NC_000913/228860‑229111                                                                                                        |                                                                                                                                                     gggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACtta                                                                                                         <  1:313004/149‑1 (MQ=32)  ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTTCGTTCCGCCACttat                                                                                                        >  1:268176/1‑149 (MQ=17)  ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACttat                                                                                                        <  2:289411/149‑1 (MQ=32)     aCTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTaa                                                                                                     <  2:50480/149‑1 (MQ=35)               aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGa                                                                                           >  2:567791/1‑149 (MQ=255)               aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGa                                                                                           <  4:271184/149‑1 (MQ=255)                     ttAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAAc                                                                                     >  2:313772/1‑149 (MQ=255)                       aaGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGc                                                                                   >  4:35681/1‑149 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTaa                                                                                      >  4:59131/1‑141 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTaa                                                                                      <  3:59131/141‑1 (MQ=255)                                gCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGttt                                                                          <  3:683543/149‑1 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                         >  2:672556/1‑149 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                         >  2:162972/1‑149 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                         >  4:161703/1‑149 (MQ=255)                                   ggTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttttt                                                                       <  1:567791/149‑1 (MQ=255)                                   ggTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGGTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTGATTAAGAAGCCTCGAGTTAACGCTCGAGGtttttt                                                                       >  1:626372/1‑149 (MQ=255)                                      cATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTc                                                                    <  1:78452/149‑1 (MQ=255)                                           aaCCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTg                                                               >  3:154655/1‑149 (MQ=255)                                                         aaaGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTg                                                 >  2:621182/1‑149 (MQ=255)                                                          aaGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGc                                                >  1:87053/1‑149 (MQ=255)                                                           agcaTTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTaa                                                                                                     >  3:598381/4‑95 (MQ=255)                                                           aGAATTCGGTGGAGCGGTCGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACttattat                                                                                                     <  4:598381/95‑2 (MQ=255)                                                           aGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                         <  4:333977/123‑1 (MQ=255)                                                           aGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                         >  3:333977/1‑123 (MQ=255)                                                               ttCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGaa                                           <  4:154655/149‑1 (MQ=255)                                                               ttCGGTAGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGaa                                           >  2:324007/1‑149 (MQ=255)                                                                 cGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATc                                         <  1:162972/149‑1 (MQ=255)                                                                            tAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCct                              <  1:573202/149‑1 (MQ=255)                                                                                 cAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCa                         >  4:405752/1‑149 (MQ=255)                                                                                    tCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCAttt                      >  4:310862/1‑149 (MQ=255)                                                                                     cGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCAtttt                     <  3:405752/149‑1 (MQ=255)                                                                                      ggTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTa                    <  4:364225/149‑1 (MQ=255)                                                                                       gTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTAc                   >  3:4250/1‑149 (MQ=255)                                                                                              tACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCttt            <  1:672556/149‑1 (MQ=255)                                                                                              cccctgcttgACAGGCAGGGGGTCGCGGGTTCGAGTTCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGa                                            <  2:572562/115‑1 (MQ=38)                                                                                              cccctgcgtgACAGGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGa                                            >  1:572562/3‑117 (MQ=255)                                                                                               ccctgcgtgACAGGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCAtttt                     >  3:447387/2‑139 (MQ=255)                                                                                               ccctgcgtgACAGGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCAtttt                     <  4:447387/138‑1 (MQ=255)                                                                                                 ctgcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTc         <  4:461753/149‑1 (MQ=255)                                                                                                      tgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCaa    <  1:313772/149‑1 (MQ=255)                                                                                                        tCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACa  >  3:525998/1‑149 (MQ=255)                                                                                                        |                                                                                                                                                     GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACA  >  NC_000913/228860‑229111

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