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	97.0%	110.4 / ‑3.2	33	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (17/15);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCG  >  NC_000913/3982541‑3982807                                                                                                                                          |                                                                                                                                     ttGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑cccc                                                                                                                                 <  2:366283/140‑1 (MQ=255)  tgtgATTCCAGTTGTCGTGGGTTCGAATCCCATGAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑ccccc                                                                                                                                >  2:878436/1‑140 (MQ=255)    tgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑ccccccc                                                                                                                              >  1:540757/1‑140 (MQ=255)    tgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑cccc                                                                                                                                 >  1:427317/1‑137 (MQ=255)    tgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑cccc                                                                                                                                 <  2:427317/137‑1 (MQ=255)            gTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGcac                                                                                                                       >  2:979243/1‑140 (MQ=255)              tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGcacca                                                                                                                     >  2:332126/1‑140 (MQ=255)              tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGcacca                                                                                                                     <  2:540757/140‑1 (MQ=255)              tGTCGTGGGTTAGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGcacca                                                                                                                     <  1:948192/140‑1 (MQ=255)                      gTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTa                                                                                                             <  2:316818/140‑1 (MQ=255)                                          cccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTc                                                                                         >  1:38951/1‑140 (MQ=255)                                             attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                                      <  1:878436/140‑1 (MQ=255)                                              ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaa                                                                                     <  2:431045/140‑1 (MQ=255)                                                    gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                                                              <  2:223876/140‑1 (MQ=255)                                                     aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                                                              <  1:668704/140‑1 (MQ=255)                                                       gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                                                            >  2:880431/1‑140 (MQ=255)                                                         tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                                                                          >  1:861684/1‑140 (MQ=255)                                                         tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                                                                          >  1:654825/1‑140 (MQ=255)                                                          gtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                                         >  2:645386/1‑140 (MQ=255)                                                           tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                                                                        >  1:753165/1‑140 (MQ=255)                                                               aaTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAg                                                                    >  2:603605/1‑140 (MQ=255)                                                                    gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                                                               <  2:861684/140‑1 (MQ=255)                                                                           gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt                                                         >  1:209630/1‑140 (MQ=255)                                                                             ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGt                                                      <  2:562546/140‑1 (MQ=255)                                                                               aaTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg                                                    <  2:810640/140‑1 (MQ=255)                                                                                    gTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa                                               >  2:195520/1‑140 (MQ=255)                                                                                         ggcgcTACCTTCAGTTGTTAGTTCCCTTACGGACGTGGGGGTCCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGt                                          <  1:734193/139‑1 (MQ=255)                                                                                             cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                                      >  1:924678/1‑140 (MQ=255)                                                                                             cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                                      >  1:887979/1‑140 (MQ=255)                                                                                             cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                                      <  1:603605/140‑1 (MQ=255)                                                                                                gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGAc                                   <  1:254149/140‑1 (MQ=255)                                                                                                   tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAg                                <  1:332126/140‑1 (MQ=255)                                                                                                            aGTGTCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGa                       <  1:195520/140‑1 (MQ=255)                                                                                                                tCCTTACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGtt                   >  2:135473/1‑140 (MQ=255)                                                                                                                    tACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGaa               >  1:929745/1‑140 (MQ=255)                                                                                                                    tACGGACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGaa               <  2:654825/140‑1 (MQ=255)                                                                                                                        gACGTGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCct           >  1:220550/1‑140 (MQ=255)                                                                                                                            tGGGGGTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCctctct       >  1:588215/1‑140 (MQ=255)                                                                                                                                 gTTCAAG‑TC‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCg  <  1:645386/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     TTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑T‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCG  >  NC_000913/3982541‑3982807

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