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	88.2%	152.1 / 13.6	51	noncoding (80/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base G (23/22);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGA  >  NC_000913/4605950‑4606216                                                                                                                                          |                                                                                                                                     caAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAGAAATCAAGGAAGGAACAAGAAAGGACGTAAAGATAATTGGTGCGA‑ggg                                                                                                                                   >  2:323827/1‑140 (MQ=255)      aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATAAAGGAAGAGACATGAATGGAAGTAAAGATAATTGGTGCGA‑ggg‑ggggg                                                                                                                             >  2:271316/1‑140 (MQ=255)        ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑ggg‑gggggga                                                                                                                           >  1:433034/1‑138 (MQ=255)                ttACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAAcc                                                                                                                   >  2:438272/1‑140 (MQ=255)                    aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCa                                                                                                               <  2:433034/140‑1 (MQ=255)                            aCCATCGTGCCTGGCGCTTATTCCTGGCATCATGCTGCTGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTa                                                                                                       <  1:323827/140‑1 (MQ=255)                             ccATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAAGTGGAGTGA‑GGG‑GGGGGTACTTGCACCCCCACGTCCGTaa                                                                                                      >  2:366420/1‑140 (MQ=255)                              cATCGTGCCTGGCACTTTTTGCCGGCATCATACTGCTGAAAATGAATTCGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAg                                                                                                     <  2:231311/140‑1 (MQ=38)                              cATCGTGCCTGGCACTTCTTGCTGGCATCATTCTGCTGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAg                                                                                                     <  2:380703/140‑1 (MQ=39)                              cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAg                                                                                                     <  1:657385/140‑1 (MQ=255)                               aTCGTGCCTGGCACTTATTGCTTGCATCAGGCTGATGTTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAgg                                                                                                    <  1:306284/140‑1 (MQ=39)                               aTCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAgg                                                                                                    <  2:302997/140‑1 (MQ=255)                                tCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGa                                                                                                   >  1:565506/1‑140 (MQ=255)                                     ccTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                             <  2:701140/140‑1 (MQ=255)                                     ccTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                             >  1:701140/1‑140 (MQ=255)                                     ccTGGCACTTATTGCTGGCATAATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGG‑GGGCGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                              >  2:323341/1‑140 (MQ=255)                                         gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACTCCTACGTCCGTAAGGACACTAacac                                                                                          >  1:380703/1‑140 (MQ=255)                                         gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTAAACGCCAACGTCCGTAAGGACACTAacac                                                                                          >  1:231311/1‑140 (MQ=39)                                               aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGACGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                                                    >  2:618724/1‑140 (MQ=255)                                                   cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCCTAGGGACACGAACACCTGAAGCTAg                                                                                >  1:320456/1‑140 (MQ=38)                                                    tGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGCAACGACACTTACACCTGAAGCTAgc                                                                               >  1:318200/1‑140 (MQ=38)                                                     ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                <  1:50307/138‑1 (MQ=255)                                                     ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                >  2:50307/1‑138 (MQ=255)                                                     ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                             >  1:504079/1‑140 (MQ=255)                                                       catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑TCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                            >  1:628007/1‑140 (MQ=255)                                                       catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                            >  1:669565/1‑140 (MQ=255)                                                       catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                            >  1:116069/1‑140 (MQ=255)                                                           atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa                                                                        >  2:91265/1‑140 (MQ=255)                                                            tGCTGATGATAATGAATTTGAATCAGGGAGGAAAAATGAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCAGAAGCTAGCGCGTCTAc                                                                       <  1:366420/140‑1 (MQ=25)                                                                     tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGTCACTAACGCCTGAATCTAGCGTGTCTACCAATTCCGc                                                              >  1:246123/1‑140 (MQ=25)                                                                     tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCAAGTCCGTAAGGACTCTACCACCTGAACCTAGGGCGTCTACCAATTCCGc                                                              >  2:326953/1‑140 (MQ=25)                                                                      aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                             <  2:116069/140‑1 (MQ=255)                                                                        tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                           >  2:97188/1‑140 (MQ=255)                                                                           aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCACCACCtt                                                        >  2:143610/1‑140 (MQ=255)                                                                             ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                      <  1:438272/140‑1 (MQ=255)                                                                              tGAATCAGGGATGAAAACTCAACGAAGAAACTAGAAAGGTAGTAAAGAAAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                     <  2:320456/140‑1 (MQ=18)                                                                               gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                    <  2:669565/140‑1 (MQ=255)                                                                               gAATCAGGGATGAAAAAGCAAGTAAGAAACATGAAAGGAAGTAAAGATTATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                    <  2:246123/140‑1 (MQ=18)                                                                                aaTCAGGGATGAAAAATCAAGGGAGAAACAAGAAAGGTAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                   <  1:326953/140‑1 (MQ=255)                                                                                      ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                   <  1:155735/134‑1 (MQ=255)                                                                                      ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                   >  2:155735/1‑134 (MQ=255)                                                                                        aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                           >  1:688687/1‑140 (MQ=255)                                                                                           aaaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTAc                                        >  1:78720/1‑140 (MQ=255)                                                                                                   aGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTACCGCGTCTACCAATTCCTCCACCTTCACACACTCATCTTACTTTTTTTg                                >  2:331356/1‑140 (MQ=21)                                                                                                       agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGGGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGatat                            >  2:448147/1‑140 (MQ=255)                                                                                                            cAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCt                       <  2:565506/140‑1 (MQ=255)                                                                                                                aaaGGAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                   <  1:97188/140‑1 (MQ=255)                                                                                                                    gAAGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt               >  2:416484/1‑140 (MQ=255)                                                                                                                      aGTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc             <  1:618724/140‑1 (MQ=255)                                                                                                                       gTAAAGATAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg            <  1:448147/140‑1 (MQ=255)                                                                                                                       gTAAAGATAATTGGTACGA‑GGGTGGAGG‑ATTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg            <  1:331356/140‑1 (MQ=21)                                                                                                                        tAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg          <  2:504079/140‑1 (MQ=255)                                                                                                                              tAATTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCGTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg     >  2:512905/1‑140 (MQ=255)                                                                                                                                aTTGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg   <  2:78720/140‑1 (MQ=255)                                                                                                                                 ttGGTGCGAGGGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggggg  <  1:416484/140‑2 (MQ=255)                                                                                                                                          |                                                                                                                                     CAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGG‑GGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGA  >  NC_000913/4605950‑4606216

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