Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,331	A→G	noncoding (71/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,331	0	A	G	94.7%	102.5 / ‑3.2	40	noncoding (71/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (18/18);  total (21/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.30e-01

TTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAAT  >  NC_000913/4606183‑4606471                                                                                                                                                        |                                                                                                                                             tgatatcgcctcgtttggtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCt                                                                                                                                            >  2:289866‑M2/26‑151 (MQ=255)      tcgcctcgtttggtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAAc                                                                                                                                       >  1:607080‑M2/21‑151 (MQ=255)         cctcgtttggtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAAcccc                                                                                                                                    >  1:538170‑M2/18‑151 (MQ=255)                tggtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTcc                                                                                                                             >  1:179456‑M2/11‑151 (MQ=255)                tggtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTGGAACCCCCACGTcc                                                                                                                             >  2:652661‑M2/11‑151 (MQ=255)                  gtgcgaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCCACGTcc                                                                                                                             <  2:150262‑M2/142‑1 (MQ=255)                      gaggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATGCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg                                                                                                                       <  1:157817‑M2/147‑1 (MQ=255)                        ggGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg                                                                                                                       >  1:150302‑M2/3‑150 (MQ=255)                          ggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                                                                                   >  1:457490/1‑151 (MQ=255)                                ggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGCACACTAACAcc                                                                                                            >  1:538627/1‑151 (MQ=255)                                      gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                                                                       <  1:635237/151‑1 (MQ=255)                                           cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCAACGTGGTTCTT‑GCTGCGA‑‑GGCGCGGTCCTTGAACCCCCACGTCCGTAAGGATACTTACTCCTGGAGCGAg                                                                                                   >  2:60726/1‑150 (MQ=21)                                           cccACGTCCATAAGAACACTAACACCTGAAGCTAGCGCTTCTACCAATTCCGCCACCTTCGCCCAGTGCGAACAATATCAACGTGGTTTTT‑GGTGCGAGGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                                    <  2:484964/151‑1 (MQ=25)                                                gTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                                                                              >  1:639868/1‑151 (MQ=255)                                                         aaCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                                    >  1:91516/1‑151 (MQ=255)                                                           cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                                  >  2:397561/1‑151 (MQ=255)                                                                     ggAAGCTGGCGCGTCTACCAATTCCGCAACCTTGGCCCAGTGCGGGCAATATCAACGGGGTTTTT‑GG‑GGGG‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                        <  1:486698/149‑1 (MQ=18)                                                                               gcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCAACTTCCGTAAGACCACTAACACCTGAAGCTAGGGGGTCTACCAATTCAGCAACCTTCGCATACCAt                                                                 >  2:18231/1‑149 (MQ=21)                                                                                cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                                              <  1:501228/151‑1 (MQ=255)                                                                                  tgtACCAATGCTGACAACGTCGGCAAGTGCCAGCAAAATCAACGTGGATTTT‑GGGGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt                                                             <  2:349268/148‑1 (MQ=17)                                                                                      ccAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTa                                                        <  1:218230/150‑1 (MQ=255)                                                                                      ccAAGTCCGCCACCTTCGGCCAGTGCGAAGAAAGAAAACGTGGTTTTTGGGGGGGG‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt                                                             <  1:441627/146‑1 (MQ=17)                                                                                         aTTCCGCCACCTTCGCCCAATGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                                                    <  1:622087/151‑1 (MQ=255)                                                                                            ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGa                                                  >  2:306122/1‑150 (MQ=255)                                                                                               ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                              >  1:256888/1‑151 (MQ=255)                                                                                                    ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAcc                                         >  1:88728/1‑151 (MQ=255)                                                                                                          cAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCAc                                       <  2:668514/148‑1 (MQ=255)                                                                                                                    aaTATAATCGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCt                          <  2:353003/150‑1 (MQ=255)                                                                                                                    aaTAACAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa                         <  2:203753/151‑1 (MQ=255)                                                                                                                       atCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGATCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt                       <  1:226241/150‑1 (MQ=255)                                                                                                                       atCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTg                      >  1:553165/1‑151 (MQ=255)                                                                                                                        tCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTgg                     >  1:52269/1‑151 (MQ=255)                                                                                                                          aaCGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt                   <  1:542769/151‑1 (MQ=255)                                                                                                                          aaCGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt                   <  1:479137/151‑1 (MQ=255)                                                                                                                          aaCGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt                   >  2:148202/1‑151 (MQ=255)                                                                                                                                 taatt‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCTCATTACTTTAAAAGAATTCCTACCGCGGAGGGGCATTCGAGTGGTTTTc                >  2:21972/4‑147 (MQ=14)                                                                                                                                    tt‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtt           <  2:580436/149‑1 (MQ=255)                                                                                                                                       ggTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCg       <  1:555898/151‑1 (MQ=255)                                                                                                                                         tGCGAG‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGt      >  2:24259/1‑151 (MQ=255)                                                                                                                                              agggggggggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAAt  <  1:658691/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                             TTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAAT  >  NC_000913/4606183‑4606471

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