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	100.0%	71.9 / NA	22	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/7);  total (15/7)

GGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCC  >  NC_000913/3982561‑3982797                                                                                                                     |                                                                                                                           ggTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGATGTGGGGGTTCAAGTCCCCCCCC‑CTCGcacca                                                                                                           >  2:301816/1‑133 (MQ=255)                      cccATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTa                                                                                      >  2:6614/1‑133 (MQ=255)                                gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                                                            >  1:383263/1‑133 (MQ=255)                                  aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                                                          >  2:448948/1‑133 (MQ=255)                                     tgtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACAACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                                                       >  1:485798/1‑133 (MQ=255)                                       tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGt                                                                     >  1:496709/1‑133 (MQ=255)                                                  aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAg                                                          <  1:28415/133‑1 (MQ=255)                                                       gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                                                     >  2:493834/1‑133 (MQ=255)                                                               ggTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAAAAGTGAACTAAAAATTCAAAAGGTAATATTTGTGGGAGTAACGCAGCtttg                                             >  1:324966/1‑131 (MQ=18)                                                                     cgcgCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgcgc                                       >  1:43048/1‑133 (MQ=255)                                                                       cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCTCGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa                                     <  1:293568/133‑1 (MQ=255)                                                                       cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACAACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa                                     <  2:485798/133‑1 (MQ=255)                                                                        gcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc                                    <  2:45474/133‑1 (MQ=255)                                                                                cAGGTGATAGTAGCCTTACGGTGGTGGGGGTTAAAGT‑CCCCCCCGCTCGCACCGCGACTTTAAAGAATTGAACTAAAAATTTAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTgg                            <  2:279580/133‑1 (MQ=255)                                                                                     gtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCa                       <  1:361449/133‑1 (MQ=255)                                                                                        aGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTg                    >  2:456889/1‑133 (MQ=255)                                                                                        aGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTg                    >  1:355000/1‑133 (MQ=255)                                                                                              cTTACGGACGTGGGGGTTCCAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTAAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAACGCAGCTTGGTAACGCCACTGGTTTAGGACGTGTGGGTCgg              >  2:247937/1‑133 (MQ=255)                                                                                              cTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAATTAAAAATTCAAAAAGCAGAATTTCGGCGAGAAGCGCAGCTTGGAACCGCAACTGGTTTGGGACCAGTGGGCCgg              >  2:400136/1‑133 (MQ=38)                                                                                              cTTACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCgg              >  1:502781/1‑133 (MQ=255)                                                                                               ttACGGACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGa             <  1:448948/133‑1 (MQ=255)                                                                                                   ggACGTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGtt         <  1:6614/133‑1 (MQ=255)                                                                                                       gTGGGGGTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGaa     <  1:312185/133‑1 (MQ=255)                                                                                                          ggggTTCAAGTCCCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATcc  >  1:68655/1‑133 (MQ=255)                                                                                                                     |                                                                                                                           GGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCC‑CTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCC  >  NC_000913/3982561‑3982797

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: ‑