Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,844,805	T→C	L281P (CTC→CCC)	adeD →	cryptic adenine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,844,805	0	T	C	100.0%	68.3 / NA	20	L281P (CTC→CCC)	adeD	cryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/11);  total (9/11)

ATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCTCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTC  >  NC_000913/3844676‑3844917                                                                                                                                  |                                                                                                                 aTCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGa                                                                                                         <  2:188752/139‑1 (MQ=255)       tGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCg                                                                                                   <  1:13538/139‑1 (MQ=255)        ggAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGt                                                                                                  >  1:275890/1‑139 (MQ=255)           agaagGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                               <  1:64493/139‑1 (MQ=255)           agaagGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                               <  2:275620/139‑1 (MQ=255)                gACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                               >  2:62444/1‑134 (MQ=255)                gACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                               <  1:62444/134‑1 (MQ=255)                   gCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGga                                                                                       <  2:207977/139‑1 (MQ=255)                    ccGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGgag                                                                                      >  1:462640/1‑139 (MQ=255)                              cAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATg                                                                            <  2:450195/139‑1 (MQ=255)                                                                              ctcAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAAtg                            <  2:462640/139‑1 (MQ=255)                                                                                                ttGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGc                                                                                 <  1:199967/68‑1 (MQ=255)                                                                                                ttGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGc                                                                                 >  2:199967/1‑68 (MQ=255)                                                                                                ttGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCat          >  2:76838/1‑139 (MQ=255)                                                                                                     cAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCg                                        >  2:289611/1‑104 (MQ=255)                                                                                                     cAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCg                                        <  1:289611/104‑1 (MQ=255)                                                                                                        cGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTc  >  2:109151/1‑139 (MQ=255)                                                                                                        cGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTc  >  1:203032/1‑139 (MQ=255)                                                                                                                      cGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGa                                                               <  1:452080/64‑1 (MQ=255)                                                                                                                      cGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGa                                                               >  2:452080/1‑64 (MQ=255)                                                                                                                                  |                                                                                                                 ATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCTCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTC  >  NC_000913/3844676‑3844917

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