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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,208	0	G	.	100.0%	163.7 / NA	43	noncoding (79/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/23);  total (20/23)

ATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606074‑4606341                                                                                                                                        |                                                                                                                                      aTAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                    <  1:742181/139‑1 (MQ=255)   aaTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                                                                                                 >  2:1336978/1‑139 (MQ=255)        gTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTg                                                                                                                           >  2:606980/1‑139 (MQ=255)         tGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACCGTCATCTTACTTTTTTTGATATCTCTTCGTTTGGTGCGAgggggggccttga                                                                                                                           >  1:416321/1‑133 (MQ=18)           cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                    <  2:142054/128‑1 (MQ=255)           cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                    >  1:142054/1‑128 (MQ=255)           cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAAc                                                                                                                        >  1:577221/1‑139 (MQ=255)                 ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                   <  1:606980/139‑1 (MQ=255)                   ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 >  2:1208848/1‑139 (MQ=255)                    gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCACCTTACTTTTTTTGATATCGCCTCGTTGGGTGCG‑GGGGGGGGACTTGAACCCTCACGt                                                                                                                >  1:898831/1‑139 (MQ=21)                     ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTc                                                                                                               <  2:663006/139‑1 (MQ=255)                           gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                         >  1:577869/1‑139 (MQ=255)                            aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                        <  1:436700/139‑2 (MQ=255)                                 ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                   <  1:1336978/139‑1 (MQ=255)                                   aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                 <  1:1208848/139‑1 (MQ=255)                                         gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTg                                                                                           >  2:1324184/1‑139 (MQ=255)                                               acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGa                                                                                                                          >  2:365870/1‑102 (MQ=255)                                               acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGa                                                                                                                          <  1:365870/102‑1 (MQ=255)                                                 acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTCTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                                                                   <  1:1019124/139‑1 (MQ=255)                                                    aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                <  2:577869/139‑1 (MQ=255)                                                        ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAAcc                                                                                                                       <  2:1503263/96‑1 (MQ=255)                                                        ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAAcc                                                                                                                       >  1:1503263/1‑96 (MQ=255)                                                                 aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCCCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                   >  1:522583/1‑139 (MQ=255)                                                                  gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:522583/139‑1 (MQ=255)                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 <  2:1170592/91‑1 (MQ=255)                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 >  1:1170592/1‑91 (MQ=255)                                                                            cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:1530858/139‑1 (MQ=255)                                                                             aaTTACGCAACCTAAACACAGTCATCTTACTTTTATTGATATCGGCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                       <  1:502512/139‑1 (MQ=17)                                                                                   gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 <  1:1324184/139‑1 (MQ=255)                                                                                          tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggactt                                                                                                                             >  1:1444985/1‑53 (MQ=255)                                                                                          tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggactt                                                                                                                             <  2:1444985/57‑5 (MQ=255)                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 >  1:700229/1‑60 (MQ=255)                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 <  2:700229/60‑1 (MQ=255)                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                  >  1:365083/1‑139 (MQ=255)                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  >  2:1230588/1‑106 (MQ=18)                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:1230588/106‑1 (MQ=18)                                                                                                   cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCCACGTgg                                 >  1:524065/1‑139 (MQ=255)                                                                                                    aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                >  1:1613284/1‑139 (MQ=255)                                                                                                          cTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 >  2:821770/1‑52 (MQ=255)                                                                                                          cTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                 <  1:821770/52‑1 (MQ=255)                                                                                                            ttttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                        >  1:626637/1‑139 (MQ=255)                                                                                                            ttttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                        >  1:157276/1‑139 (MQ=255)                                                                                                            ttttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                        <  2:898831/139‑1 (MQ=255)                                                                                                             tttttGATACCGCCACGTTGGGGGCG‑GGGGGGGGACTTGAACCCCCACGCCCGTAAGGCGACTAACACCTGAAGCTAGCGCTTCTACCGATACCGCCACCCTCGCCCAGTGCGAGCAATATCAACGTGGTTATTGGTGc                       <  1:949374/139‑1 (MQ=18)                                                                                                               tttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   >  2:371737/1‑109 (MQ=255)                                                                                                               tttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   <  1:371737/109‑1 (MQ=255)                                                                                                                   atatCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                 <  2:416321/139‑1 (MQ=255)                                                                                                                    tatCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg                <  2:577221/139‑1 (MQ=255)                                                                                                                         ccTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                   >  2:1380351/1‑35 (MQ=17)                                                                                                                         ccTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                   <  1:1380351/35‑1 (MQ=17)                                                                                                                         ccTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGAc           >  2:809351/1‑139 (MQ=255)                                                                                                                                  gTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAcccc  <  2:626637/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      ATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606074‑4606341

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