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	95.7%	163.0 / 0.0	47	noncoding (80/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (25/20);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGA  >  NC_000913/4605963‑4606216                                                                                                                             |                                                                                                                                    cATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGa                                                                                                                      <  2:181289/139‑1 (MQ=255)        ttttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCa                                                                                                               <  2:145687/139‑1 (MQ=255)          ttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACg                                                                                                             <  1:401550/139‑1 (MQ=255)                ccATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATG‑AAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTaa                                                                                                      >  1:564460/1‑139 (MQ=255)                     gTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                  >  2:521806/1‑139 (MQ=255)                     gTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                  >  1:778883/1‑139 (MQ=255)                     gTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                  >  2:884223/1‑139 (MQ=255)                        ccTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                               <  2:109018/139‑1 (MQ=255)                        ccTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                               <  1:521806/139‑1 (MQ=255)                          tGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                             <  2:1099040/139‑1 (MQ=255)                            gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                           >  2:351628/1‑139 (MQ=255)                            gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                           >  1:715114/1‑139 (MQ=255)                                      cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                 <  1:846930/139‑1 (MQ=255)                                        ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                                                               >  1:204868/1‑139 (MQ=255)                                          catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                             >  1:26691/1‑139 (MQ=255)                                                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                    >  1:510403/1‑139 (MQ=255)                                                     tgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                 >  1:491992/1‑139 (MQ=255)                                                         aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              >  1:157487/1‑139 (MQ=255)                                                         aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  1:884223/139‑1 (MQ=255)                                                              aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  2:1024693/95‑1 (MQ=255)                                                              aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     >  1:1024693/1‑95 (MQ=255)                                                                ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                       >  1:434458/1‑139 (MQ=255)                                                                ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                       >  2:531333/1‑139 (MQ=255)                                                                ctGAATCAGGGATGAAAACTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGAGGGGGGACTTGAACCCCCACGTCCGTAAGGGCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCAACCTTc                                                       >  1:1032211/2‑139 (MQ=21)                                                                  gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              >  2:135473/1‑130 (MQ=255)                                                                  gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  1:135473/130‑1 (MQ=255)                                                                  gAATCAGGGATG‑AAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                    <  2:564460/139‑1 (MQ=255)                                                                   aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                    <  2:434458/139‑1 (MQ=255)                                                                       aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  1:1018366/125‑1 (MQ=255)                                                                       aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              >  2:1018366/1‑125 (MQ=255)                                                                         ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                               <  2:73746/90‑1 (MQ=255)                                                                         ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                               >  1:73746/1‑90 (MQ=255)                                                                           aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATc                                            <  1:128768/139‑1 (MQ=255)                                                                             gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt                                          >  1:763781/1‑139 (MQ=255)                                                                                            aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATc                           >  2:1083455/1‑139 (MQ=255)                                                                                                       gAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     >  2:3992/1‑54 (MQ=18)                                                                                                       gAAGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  1:3992/54‑1 (MQ=18)                                                                                                        aaGTAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt               <  2:715114/139‑1 (MQ=255)                                                                                                           tAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg            <  2:204868/139‑1 (MQ=255)                                                                                                           tAAAGATAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg            <  1:351628/139‑1 (MQ=255)                                                                                                                aTAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg       >  2:4364/1‑139 (MQ=255)                                                                                                                 tAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTAATTTTTTTTTTATCGCCTCGTTTGGTGGGAggggg      >  1:1004209/1‑139 (MQ=18)                                                                                                                 tAATTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGCCCGTAAGAAAACTAACACCTGAACCTACCGCGTCTACCAATTCCCCCACCTTCGCACGGCCATCTAAATTTTTTTGATATCGCCTGGTTTGTTGCGAggggg      >  2:249740/1‑139 (MQ=17)                                                                                                                  aaTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                <  2:458803/80‑1 (MQ=255)                                                                                                                  aaTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                >  1:458803/1‑80 (MQ=255)                                                                                                                   aTTGGGGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg    <  1:249740/139‑1 (MQ=255)                                                                                                                     tGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGa  >  1:371454/1‑139 (MQ=255)                                                                                                                     tGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGa  >  1:806312/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                    CATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGA  >  NC_000913/4605963‑4606216

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