Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,844,470	A→G	E169E (GAA→GAG)	adeD →	cryptic adenine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,844,470	0	A	G	100.0%	126.6 / NA	34	E169E (GAA→GAG)	adeD	cryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/20);  total (14/20)

CAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCAC  >  NC_000913/3844332‑3844605                                                                                                                                           |                                                                                                                                        cAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGcc                                                                                                                             >  1:195888/1‑151 (MQ=255)  aaCGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCt                                                                                                                            <  2:541355/151‑1 (MQ=255)    cGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTgg                                                                                                                          <  1:143927/151‑1 (MQ=255)    cGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTgg                                                                                                                          >  1:554066/1‑151 (MQ=255)            gCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCa                                                                                                                  <  1:469269/151‑1 (MQ=255)                   cGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGc                                                                                                             >  1:225185/1‑149 (MQ=255)                    ggATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGGGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGc                                                                                                             >  2:391841/1‑148 (MQ=255)                      aTTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGtt                                                                                                        >  2:77585/1‑151 (MQ=255)                       ttCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGtt                                                                                                        >  1:330179/1‑150 (MQ=255)                       ttCGCCTGGTTTGCCCGCTGTGCCGAACAGGAAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGtt                                                                                                        <  2:502102/150‑1 (MQ=255)                         cGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTAcc                                                                                                     <  1:107380/151‑1 (MQ=255)                               gTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCtt                                                                                               >  1:300958/1‑151 (MQ=255)                                 ttGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTg                                                                                              <  1:30982/150‑1 (MQ=255)                                      cGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGaaa                                                                                        <  2:534928/151‑1 (MQ=255)                                          gtgCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAatga                                                                                     <  1:244372/150‑1 (MQ=255)                                           tgCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAatgatg                                                                                   >  1:405787/1‑151 (MQ=255)                                            gCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATgg                                                                                  <  2:37771/151‑1 (MQ=255)                                              cGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGAc                                                                                >  1:116877/1‑151 (MQ=255)                                                 aCAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTAc                                                                             <  2:625412/151‑1 (MQ=255)                                                    ggcaaggcaaAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCt                                                                          <  1:443061/151‑1 (MQ=255)                                                           caaAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTa                                                                    >  1:87830/1‑150 (MQ=255)                                                                  aGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGGTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTATCgg                                                            <  2:165494/151‑1 (MQ=255)                                                                     aCTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGc                                                          >  1:478357/1‑150 (MQ=255)                                                                          cAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATg                                                    <  2:347721/151‑1 (MQ=255)                                                                               cAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCgct                                                <  1:214793/150‑1 (MQ=255)                                                                                aGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACTATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCgctgc                                              <  2:39164/151‑1 (MQ=255)                                                                                     ttGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGAt                                         >  1:323328/1‑151 (MQ=255)                                                                                              cGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGAt                                <  1:162620/151‑1 (MQ=255)                                                                                               gCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATg                               >  2:399333/1‑151 (MQ=255)                                                                                                    ggAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCAttt                          >  1:124909/1‑151 (MQ=255)                                                                                                            gCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCAcc                   <  1:384253/150‑1 (MQ=255)                                                                                                            gCGATGTTAACGGTGCCAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCAc                    <  1:281245/149‑1 (MQ=255)                                                                                                                            cAGTTTTACCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCAc  <  1:83435/151‑1 (MQ=255)                                                                                                                             aGTTTTGCCCTTGAGCAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCAc  <  1:344991/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        CAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCAGTTTTACCCTTGAACAGATGCTCGCCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCAC  >  NC_000913/3844332‑3844605

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