Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,844,146	T→C	G61G (GGT→GGC)	adeD →	cryptic adenine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,844,146	0	T	C	100.0%	156.4 / NA	44	G61G (GGT→GGC)	adeD	cryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/24);  total (20/24)

TAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCC  >  NC_000913/3843999‑3844289                                                                                                                                                    |                                                                                                                                                tAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGtt                                                                                                                                              >  1:207998/1‑151 (MQ=255)                    gAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATg                                                                                                                           <  2:80465/151‑1 (MQ=255)                          ttAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                         >  2:186198/1‑131 (MQ=255)                          ttAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                         <  1:186195/131‑1 (MQ=255)                           tAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCgg                                                                                                                     <  1:162012/150‑1 (MQ=255)                                         ggCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGAtt                                                                                                       <  2:664680/150‑1 (MQ=255)                                            gACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATg                                                                                                   >  1:279097/1‑151 (MQ=255)                                               gCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTc                                                                                                >  2:460241/1‑151 (MQ=255)                                                 tGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcgc                                                                                              >  1:308881/1‑151 (MQ=255)                                                  gTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                               <  2:271944/149‑1 (MQ=255)                                                    tGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCgcgg                                                                                            <  1:561436/150‑1 (MQ=255)                                                         attatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAc                                                                                      >  1:598300/1‑151 (MQ=255)                                                           tatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcgg                                                                                    >  2:591398/1‑151 (MQ=255)                                                           tatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcg                                                                                     <  1:116282/150‑1 (MQ=255)                                                            atattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggc                                                                                   <  2:291182/151‑1 (MQ=255)                                                             tattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcg                                                                                  >  2:343604/1‑151 (MQ=255)                                                                 atcgATAATGTCTCTATTCTCGACCTGCTCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcgg                                                                                 >  2:239634/4‑148 (MQ=255)                                                                          gTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGttt                                                                      <  2:486634/150‑1 (MQ=255)                                                                              ctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAtt                                                                   >  2:388485/1‑149 (MQ=255)                                                                                     cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTTTTTATGCTc                                                            >  1:100732/1‑149 (MQ=255)                                                                                     cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGATACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGCCCAGGGTTTATTgctcct                                                             >  1:537069/1‑143 (MQ=255)                                                                                         cTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTg                                                       >  2:511459/1‑150 (MQ=255)                                                                                            aTCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTTATGCTCACCCGc                                                      >  2:498320/1‑148 (MQ=255)                                                                                              cAATAGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATAtt                                                 >  1:398463/1‑151 (MQ=255)                                                                                                      gagaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGCTTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAg                                         <  1:44125/151‑1 (MQ=255)                                                                                                         aaaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCat                                      <  2:49594/151‑1 (MQ=255)                                                                                                           aTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGTTCACCTGCATATTGAATCCAGCat                                      >  1:46774/1‑149 (MQ=255)                                                                                                             ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatgat                                   <  2:468713/150‑1 (MQ=255)                                                                                                             ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCCTATTGAATCCAGCatgat                                   >  1:273443/1‑150 (MQ=255)                                                                                                               ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCAtgatga                                 <  2:600372/150‑1 (MQ=255)                                                                                                               ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCAtgatga                                 <  1:611889/150‑1 (MQ=255)                                                                                                               ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGAc                                <  2:510373/151‑1 (MQ=255)                                                                                                                    ccAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCgg                           <  2:49311/151‑1 (MQ=255)                                                                                                                         tgtgATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACt                      >  2:27534/1‑151 (MQ=255)                                                                                                                          gtgATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACt                      <  1:384620/150‑1 (MQ=255)                                                                                                                                aaaGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaaa               <  1:590305/151‑1 (MQ=255)                                                                                                                                 aaGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAc              <  1:416376/151‑1 (MQ=255)                                                                                                                                  aGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAcc             <  1:178920/151‑1 (MQ=255)                                                                                                                                   ggACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCg            <  1:639144/151‑1 (MQ=255)                                                                                                                                          aCATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAccc      >  2:543869/1‑150 (MQ=255)                                                                                                                                          aCATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCt     <  1:107311/151‑1 (MQ=255)                                                                                                                                           cATTGCCGGCGTTGGCGCAGAATCCACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCGCCTGCATATTGAATCCAGCATGATGACGCCGGTCCCTTTTGAAACCGCTACCCt     <  1:78554/150‑1 (MQ=255)                                                                                                                                            aTTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGc   <  1:147393/151‑1 (MQ=255)                                                                                                                                              tGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTCTTGATGCTCACCTGCATATTGAATCCAGCATGGTGACGCCGGTCACTTTTGAAACCGCTACCCTTcc  >  1:64177/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                TAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCC  >  NC_000913/3843999‑3844289

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