Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,999,025	G→C	10.0%	noncoding (33/74 nt)	glyU ←	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,999,025	0	G	C	10.0%	164.7 / 9.9	50	noncoding (33/74 nt)	glyU	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base C (2/3);  total (24/26)
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.98e-01

AACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAAC  >  NC_000913/2998880‑2999173                                                                                                                                                  |                                                                                                                                                     aacTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAg                                                                                                                                                   >  1:215874/1‑149 (MQ=255)    ttGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGctc                                                                                                                                                >  1:30963/1‑149 (MQ=255)      gAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAGGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCtctc                                                                                                                                              >  4:41857/1‑149 (MQ=255)      gAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCtctc                                                                                                                                              <  2:115825/149‑1 (MQ=255)        aaaTAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctata                                                                                                                                            <  3:238574/149‑4 (MQ=255)          aTAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctatacg                                                                                                                                          >  2:229550/1‑144 (MQ=255)            aaGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTa                                                                                                                                        <  4:124949/149‑1 (MQ=255)            aaGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTa                                                                                                                                        >  1:60786/1‑149 (MQ=255)            aaGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTa                                                                                                                                        <  3:256625/149‑1 (MQ=255)             aGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTAc                                                                                                                                       >  3:338806/1‑149 (MQ=255)             aGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTAc                                                                                                                                       <  3:171278/149‑1 (MQ=255)              ggTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTAcc                                                                                                                                      >  3:110705/1‑149 (MQ=255)               gtgtTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctatacga                                                                                                                                         >  1:430668/1‑139 (MQ=255)               gtgtTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctatacga                                                                                                                                         <  2:430668/145‑7 (MQ=255)                                 tACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTc                                                                                                                   <  4:110705/149‑1 (MQ=255)                                 tACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTc                                                                                                                   <  3:15564/149‑1 (MQ=255)                                   cAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCga                                                                                                                 <  2:30963/149‑1 (MQ=255)                                   cAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCga                                                                                                                 <  4:54776/149‑1 (MQ=255)                                    aaaaGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCgag                                                                                                                >  2:170667/1‑149 (MQ=255)                                    aaaaGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCgag                                                                                                                >  4:46258/1‑149 (MQ=255)                                            tgcattgcGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTaa                                                                                                        <  2:60786/149‑1 (MQ=255)                                             tcattgcGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAGGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAg                                                                                                       <  3:41857/148‑1 (MQ=255)                                                 tgcGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCAt                                                                                                   <  4:315416/149‑1 (MQ=255)                                                  gcGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCAtt                                                                                                  >  3:274857/1‑149 (MQ=255)                                                    gTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGtt                                                                                                                                           >  2:257998/1‑106 (MQ=255)                                                    gTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGtt                                                                                                                                           <  1:257998/106‑1 (MQ=255)                                                        aCGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAAc                                                                                            >  4:19459/1‑149 (MQ=255)                                                                        acaaAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATACGTAAGGCATTATAAACCTTACGCTCTCCt                                                                               <  2:464729/146‑1 (MQ=255)                                                                        acaaAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATACGTAAGGCATTATAAACCTTACGCTCTCCt                                                                               >  1:464729/1‑146 (MQ=255)                                                                             aTGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGt                                                                       <  2:215874/149‑1 (MQ=255)                                                                                    aGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACgc                                                                <  4:333211/149‑1 (MQ=255)                                                                                     gTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACgct                                                               >  4:385980/1‑149 (MQ=255)                                                                                       aTTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACgctgc                                                             <  4:331173/149‑1 (MQ=255)                                                                                        ttACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTAGAATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACgctgct                                                            >  4:396271/1‑149 (MQ=255)                                                                                               gAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGctgctgactgct                                                     >  4:128086/1‑149 (MQ=255)                                                                                                   aTAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGctgctgactgctga                                                   >  1:116865/1‑147 (MQ=255)                                                                                                   aTAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGctgctgactgctga                                                   <  2:116865/147‑1 (MQ=255)                                                                                                     aaaTTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTaa                                               <  4:338806/149‑1 (MQ=255)                                                                                                       aTTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAAtt                                             <  3:19459/149‑1 (MQ=255)                                                                                                                         aaTCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCa                                                                                                    >  3:299179/1‑76 (MQ=255)                                                                                                                         aaTCGAACCCTCGTATAGAGCTTGGGAAACTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCa                                                                                                    <  4:299179/76‑1 (MQ=255)                                                                                                                            cGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTtg                        <  3:128086/149‑1 (MQ=255)                                                                                                                                 ccTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGcc                   <  3:46258/149‑1 (MQ=255)                                                                                                                                 ccTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGcc                   >  4:441203/1‑149 (MQ=255)                                                                                                                                   tCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCt                 >  3:368928/1‑149 (MQ=255)                                                                                                                                      acgagagCTTCCCAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGAt                                                  <  4:56517/110‑1 (MQ=255)                                                                                                                                         agagCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGc                                                                                                            <  4:12892/52‑1 (MQ=255)                                                                                                                                         agagCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGc                                                                                                            >  3:12892/1‑52 (MQ=255)                                                                                                                                            gCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTTGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCa        >  2:4663/1‑149 (MQ=255)                                                                                                                                               tGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGt     >  3:464989/1‑149 (MQ=255)                                                                                                                                                 ggAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTaa   <  1:445575/149‑1 (MQ=255)                                                                                                                                                  gAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAAc  >  1:236856/1‑149 (MQ=255)                                                                                                                                                  gAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAAc  >  3:350962/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     AACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAAC  >  NC_000913/2998880‑2999173

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