Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,055,185	Δ1 bp	coding (125/1455 nt)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,055,184	0	C	.	100.0%	69.7 / NA	17	coding (124/1455 nt)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/8);  total (9/8)

ATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACC  >  NC_000913/2055104‑2055258                                                                                 |                                                                           aTAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTAc                                                                            <  1:359626/81‑1 (MQ=255) aTAAACCCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTAc                                                                            >  2:359655/1‑81 (MQ=255)    aaCTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTActg                                                                          <  2:359627/80‑3 (MQ=255)        cGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑Ctgatga                                                                     >  2:359663/1‑80 (MQ=255)         gACGCGCTGTATGAGCAATCTGTAGTCGCATTACGAAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑Ctgatga                                                                     <  2:359628/79‑1 (MQ=255)          aCGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGaaaa                                                                  >  2:359666/1‑81 (MQ=255)                      aGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAAc                                                       <  1:359629/80‑1 (MQ=255)                      aGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAAc                                                       >  1:359654/1‑80 (MQ=255)                       gCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAAc                                                       >  2:359667/1‑79 (MQ=255)                               tAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCttt                                              <  1:359631/80‑1 (MQ=255)                               tAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCt                                                <  1:359630/78‑1 (MQ=255)                                 gTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTg                                           <  1:359632/81‑1 (MQ=255)                                      aTTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTAt                                      >  1:359664/1‑81 (MQ=255)                                        tACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATc                                     >  2:359662/1‑80 (MQ=255)                                                   aTTGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCg                         >  2:359669/1‑81 (MQ=255)                                                    ttGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAATGTCTTTTTGTCTATCCATCACTGACCg                         >  2:359659/1‑80 (MQ=255)                                                    ttGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGt                        <  2:359633/81‑1 (MQ=255)                                                     tGGCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGt                        >  1:359668/1‑80 (MQ=255)                                                       gCAACTATATTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTaa                      <  1:359634/80‑1 (MQ=255)                                                            tataTTACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTgg                 >  2:359665/1‑80 (MQ=255)                                                                ttACAAGTGGCGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACg             <  1:359635/80‑1 (MQ=255)                                                                          cGAATTA‑CTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAAcc  <  2:359636/81‑1 (MQ=255)                                                                                 |                                                                           ATAAACTCGACGCGCTGTATGAGCAATCTGTAGTCGCATTACGCAACGCCATTGGCAACTATATTACAAGTGGCGAATTACCTGATGAAAACGCCCGCAAACAAGGTCTTTTTGTCTATCCATCACTGACCGTAACCTGGGACGGTAGCACAACC  >  NC_000913/2055104‑2055258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑