Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,094	(G)8→9	noncoding (72/87 nt)	leuV ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,086	1	.	G	96.8%	100.2 / ‑3.0	31	noncoding (80/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (12/18);  total (12/19)
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.91e-01

TCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGAC  >  NC_000913/4605947‑4606217                                                                                                                                              |                                                                                                                                     tCACAAACAGTTGCGGCATTAGATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATCATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑Ga                                                                                                                                       <  2:754518/140‑1 (MQ=255)    caAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggg                                                                                                                                  <  1:593188/140‑1 (MQ=255)         aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gggggggg                                                                                                                             >  2:521101/1‑140 (MQ=255)         aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gggggggg                                                                                                                             >  1:30286/1‑140 (MQ=255)         aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gggggggg                                                                                                                             >  1:64095/1‑140 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  2:801252/1‑138 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  2:662133/1‑138 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  2:381492/1‑138 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  1:368523/1‑138 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  1:737628/1‑138 (MQ=255)           ttGCGGCATTACATTTTGCTACCATCGTGCCTGACACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggggga                                                                                                                           >  2:341160/1‑138 (MQ=255)            tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gggggggggac                                                                                                                          <  1:381492/140‑4 (MQ=255)                       aTTTTGCTATAATTGTGACTGTCGCTTATTGGTGACATCATGCTGATGATAATGAATTTGAATCAGGGTTCAAAAATAAAGGAAGAAATACGAAAGGACTTAAAGATAATTGGTCC‑TAG‑GGGGGGGGACTTGCATCCCCa                                                                                                                <  2:548846/140‑1 (MQ=255)                           tGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTc                                                                                                            >  2:182037/1‑140 (MQ=255)                             cTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCg                                                                                                          >  1:73107/1‑140 (MQ=255)                                  aTCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                     <  2:823204/140‑1 (MQ=255)                                       gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                <  1:225619/140‑1 (MQ=255)                                       gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                <  1:521101/140‑1 (MQ=255)                                       gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                <  2:756/140‑1 (MQ=255)                                         cTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                              <  2:73107/140‑1 (MQ=255)                                          tGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                              <  1:555754/140‑1 (MQ=255)                                                  aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                                                     >  1:291177/1‑140 (MQ=255)                                                          catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                             <  2:227564/140‑1 (MQ=255)                                                                gCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                       >  1:531157/1‑140 (MQ=255)                                                                 cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                      <  2:291177/140‑1 (MQ=255)                                                                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                    <  1:801252/140‑1 (MQ=255)                                                                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                    >  2:321255/1‑140 (MQ=255)                                                                        tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                               <  1:57104/140‑1 (MQ=255)                                                                        tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                               <  2:531157/140‑1 (MQ=255)                                                                         aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                              >  1:584720/1‑140 (MQ=255)                                                                               tttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                        <  1:466735/140‑1 (MQ=255)                                                                                  gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                     <  2:368523/140‑1 (MQ=255)                                                                                       aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTc                                                >  2:141669/1‑140 (MQ=255)                                                                                         ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAt                                              >  1:177708/1‑140 (MQ=255)                                                                                                   tCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAACCTTACtttta                                    >  2:407093/1‑139 (MQ=255)                                                                                                    cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtttttt                                   <  1:341160/140‑1 (MQ=255)                                                                                                                       gAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                >  2:439978/1‑140 (MQ=255)                                                                                                                          gTAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg             <  1:141669/140‑1 (MQ=255)                                                                                                                           tAAAGATAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa            <  1:662133/140‑1 (MQ=255)                                                                                                                                 tAATTGGTGCGGG‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg      <  1:182037/140‑1 (MQ=255)                                                                                                                                 tAATTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg      <  1:439978/140‑1 (MQ=255)                                                                                                                                   aTTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg    <  1:321255/140‑1 (MQ=255)                                                                                                                                   aTTGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg    >  1:689547/1‑140 (MQ=255)                                                                                                                                     tGGTGC‑GAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggggga  >  1:281386/1‑138 (MQ=255)                                                                                                                                              |                                                                                                                                     TCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGAC  >  NC_000913/4605947‑4606217

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