Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	228,966	C→G	9.8%	noncoding (39/77 nt)	aspU →	tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	228,966	0	C	G	9.8%	133.1 / 6.4	41	noncoding (39/77 nt)	aspU	tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base G (2/2);  total (21/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 = 9.79e-01

GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTC  >  NC_000913/228860‑229114                                                                                                           |                                                                                                                                                     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                                                                                         <  3:59131/141‑1 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTaa                                                                                         >  4:59131/1‑141 (MQ=255)                                gCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGttt                                                                             <  3:683543/149‑1 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                            >  4:161703/1‑149 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                            >  2:162972/1‑149 (MQ=255)                                 ccGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                            >  2:672556/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                                                                            >  3:333977/1‑123 (MQ=255)                                                           aGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGtttt                                                                            <  4:333977/123‑1 (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)                                                                                                           cGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTc  <  1:621182/149‑1 (MQ=255)                                                                                                           cGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTc  <  4:4250/149‑1 (MQ=255)                                                                                                           |                                                                                                                                                     GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTC  >  NC_000913/228860‑229114

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: ‑