Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,982,684	(C)8→9	noncoding (79/87 nt)	leuT →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,676	1	.	C	95.1%	141.3 / ‑0.0	41	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (20/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGA  >  NC_000913/3982539‑3982808                                                                                                                                            |                                                                                                                                      gATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑cc                                                                                                                                    <  1:100486/140‑1 (MQ=255) gATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑cc                                                                                                                                    <  2:1000664/140‑1 (MQ=255)      tgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑ccccccc                                                                                                                               >  2:141355/1‑140 (MQ=255)              gTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGcac                                                                                                                        <  1:836990/140‑1 (MQ=255)                tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGcacca                                                                                                                      <  1:1053211/140‑1 (MQ=255)                tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGcacca                                                                                                                      <  2:665603/140‑1 (MQ=255)                tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGcacca                                                                                                                      <  2:63716/140‑1 (MQ=255)                tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGcacca                                                                                                                      <  1:318771/140‑1 (MQ=255)                  tCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑CCCCCCCCTCGCACCACGa                                                                                                                   <  2:644288/140‑1 (MQ=255)                    gTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGAc                                                                                                                  <  2:134748/140‑1 (MQ=255)                                       gCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaa                                                                                                >  1:1020587/1‑140 (MQ=255)                                             ccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                                         >  2:329574/1‑140 (MQ=255)                                                      gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                                         <  1:295030/132‑1 (MQ=255)                                                      gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                                         >  2:295030/1‑132 (MQ=255)                                                           tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGTCGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                                    <  2:448600/131‑1 (MQ=255)                                                           tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                                                                    >  1:448600/1‑131 (MQ=255)                                                             tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                                         <  2:573124/124‑1 (MQ=255)                                                             tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                                                         >  1:573124/1‑124 (MQ=255)                                                                    gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                                                                  <  1:1029325/140‑1 (MQ=255)                                                                       aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcgc                                                               >  1:110616/1‑140 (MQ=255)                                                                          gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc                                                            >  1:550140/1‑140 (MQ=255)                                                                           tGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt                                                           >  1:461906/1‑140 (MQ=255)                                                                             gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTg                                                         >  2:573397/1‑140 (MQ=255)                                                                               ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGt                                                       <  1:843213/140‑1 (MQ=255)                                                                                     ggTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgcgc                                                  <  2:1020587/140‑1 (MQ=255)                                                                                        aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACt                                              >  2:244272/1‑140 (MQ=255)                                                                                        aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACt                                              <  2:470122/140‑1 (MQ=255)                                                                                             cgcTAGCTTCAGGTGTTAGTGTCCATACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGttt                                         <  1:703148/140‑1 (MQ=255)                                                                                               cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                                       <  2:550140/140‑1 (MQ=255)                                                                                                     tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAg                                 <  2:461906/140‑1 (MQ=255)                                                                                                            ttAGTGTCCTTACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTc                           >  1:137818/1‑140 (MQ=255)                                                                                                              aGTGTCCTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGa                        >  1:149869/1‑140 (MQ=255)                                                                                                                    cTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCg                  >  2:853118/1‑140 (MQ=255)                                                                                                                    cTTACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCg                  >  2:410196/1‑140 (MQ=255)                                                                                                                      tACGGACGTGGGGGTTCAA‑GTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGa                 >  1:444671/1‑140 (MQ=255)                                                                                                                      tACGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGaa                >  2:230344/1‑140 (MQ=255)                                                                                                                       aCGGACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAAt               <  1:244272/140‑1 (MQ=255)                                                                                                                          gACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCg                  >  1:874705/1‑134 (MQ=255)                                                                                                                          gACGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCg                  <  2:874705/134‑1 (MQ=255)                                                                                                                            cGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTc                   >  1:204430/1‑131 (MQ=255)                                                                                                                            cGTGGGGGTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTc                   <  2:204430/131‑1 (MQ=255)                                                                                                                                ggggTTCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCg      >  2:302143/1‑140 (MQ=255)                                                                                                                                    ttCAA‑GTC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGa  >  2:297990/1‑140 (MQ=255)                                                                                                                                    ttCAA‑GTC‑CCCCCCCCCCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGa  >  1:729008/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                      GATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAA‑GT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGA  >  NC_000913/3982539‑3982808

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