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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,323	0	G	.	93.5%	119.8 / 1.6	32	noncoding (79/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (15/14);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCG  >  NC_000913/4606178‑4606466                                                                                                                                                  |                                                                                                                                                cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAACAATATCAACGTGGTTTTTGGAGCGAgggg                                                                                                                                              >  2:652769/1‑149 (MQ=255)   ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                                                                                                                                           <  1:244202/150‑1 (MQ=255)           tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAAcc                                                                                                                                 >  1:49076/1‑151 (MQ=255)               gCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGTGCGA‑GGGGGGGACTTGAAccccc                                                                                                                              >  1:47739/1‑150 (MQ=255)                    gTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                                        >  1:687525/1‑151 (MQ=255)                        ggTGGGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTCCCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTc                                                                                                                         <  2:205108/146‑1 (MQ=255)                         gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                                   <  2:687483/151‑1 (MQ=255)                           gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTCACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTACGAGCAATATCAACGAGGTTTTTCGAGCGA‑GGGGGGGACTTGAGCCCCCACGTCCGt                                                                                                                      >  2:502956/1‑146 (MQ=18)                            cGAGGGGGGGGACTTGAACCCCAACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                                <  2:529924/151‑1 (MQ=255)                                 ggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                           <  1:56479/151‑1 (MQ=255)                                    gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac                                                                                                        <  1:469414/151‑1 (MQ=255)                                    gggACTTGAACCCCCACGACCGTAAGAACCCTAACGCCTGAAACTAACTCGACTCCGAAATCCGCCACCTAGGACCAATGCAAGCAAGATCAACGGGGTTTTTGGTGCG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                                         <  2:114037/150‑1 (MQ=21)                                     ggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                                       <  1:528335/151‑1 (MQ=255)                                             aCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                               <  1:456431/151‑1 (MQ=255)                                                  cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                          >  1:276596/1‑151 (MQ=255)                                                         gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                                    >  1:296814/1‑150 (MQ=255)                                                                       acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                     >  1:453120/1‑151 (MQ=255)                                                                        ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGAAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                                      >  2:22161/1‑149 (MQ=255)                                                                           gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                   <  1:737104/149‑1 (MQ=255)                                                                               cTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATAcc                                                              >  1:383119/1‑150 (MQ=255)                                                                                 aGCGCGTCTCCCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCa                                                             <  1:125678/149‑1 (MQ=255)                                                                                        cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGC‑AAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCtt                                                    >  1:125081/1‑151 (MQ=255)                                                                                         tACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCtt                                                    >  1:382026/1‑150 (MQ=255)                                                                                               ttCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaa                                                >  2:726615/1‑148 (MQ=255)                                                                                                 ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                                           >  1:165448/1‑151 (MQ=255)                                                                                                   gCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                         <  1:22161/151‑1 (MQ=255)                                                                                                             cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGa                               <  2:761341/151‑1 (MQ=255)                                                                                                                         aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa                    <  2:683468/150‑1 (MQ=255)                                                                                                                         aaTAACAAAATGGTTTTTGGTCCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa                    <  2:266553/150‑1 (MQ=21)                                                                                                                                aCGTGGTTTTTGGGCCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt              <  2:495849/149‑1 (MQ=255)                                                                                                                                   tGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTc           >  2:564904/1‑149 (MQ=255)                                                                                                                                    ggTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGc        >  2:398362/1‑151 (MQ=255)                                                                                                                                      ttattgggggcaGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTATTGGTTTTCAGc        <  2:379799/139‑1 (MQ=255)                                                                                                                                          tGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCg  >  2:607566/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                CTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCG  >  NC_000913/4606178‑4606466

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