Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,999,024	G→C	10.0%	noncoding (34/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,024	0	G	C	10.0%	165.4 / 9.4	50	noncoding (34/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 = 5.27e-01

AACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAA  >  NC_000913/2998880‑2999172                                                                                                                                                 |                                                                                                                                                     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:171278/149‑1 (MQ=255)             aGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTAc                                                                                                                                      >  3:338806/1‑149 (MQ=255)              ggTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTAcc                                                                                                                                     >  3:110705/1‑149 (MQ=255)               gtgtTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctatacga                                                                                                                                        >  1:430668/1‑139 (MQ=255)               gtgtTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTCCCAAGctctatacga                                                                                                                                        <  2:430668/145‑7 (MQ=255)                                 tACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTc                                                                                                                  <  3:15564/149‑1 (MQ=255)                                 tACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTc                                                                                                                  <  4:110705/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                                                                                                                                          <  1:257998/106‑1 (MQ=255)                                                    gTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGtt                                                                                                                                          >  2:257998/1‑106 (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)                                                                                                                                                 |                                                                                                                                                     AACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAAGGATTGCATTGCGTAAACGCTTTTTATTTACAACAAAATGGGGAAGTATTACGGCGAAGATAAATTGGAGCGGGCGAAGGGAATCGAACCCTCGTATAGAGCTTGGGAAGCTCTCGTTCTACCATTGAACTACGCCCGCTTCGAGATGCGTAAGGCATTATAAACCTTACGCTCTCCTTAGCAAGTGCCACGCTGCTGACTGCTGATTAATTCGCCATCAGCATTTTGGCTTGCTGCCCTGAGGCGGCAAGTAA  >  NC_000913/2998880‑2999172

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