Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,233,561	A→G	E228E (GAA→GAG)	sanA →	DUF218 domain‑containing protein SanA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,233,561	0	A	G	100.0%	82.8 / NA	25	E228E (GAA→GAG)	sanA	DUF218 domain‑containing protein SanA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/17);  total (8/17)

GCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTT  >  NC_000913/2233489‑2233635                                                                         |                                                                           gCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACt                                                                    >  2:395223/1‑81 (MQ=255)  ccGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACtt                                                                   <  2:395198/81‑1 (MQ=255)    gCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGa                                                                 <  1:395199/81‑1 (MQ=255)     cTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGaa                                                                <  2:395200/81‑1 (MQ=255)       ggTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAAt                                                               >  1:395231/1‑80 (MQ=255)        gTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTa                                                             <  2:395201/81‑1 (MQ=255)          cccTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTAc                                                            >  2:395233/1‑80 (MQ=255)           ccTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGAGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACa                                                           <  1:395202/80‑1 (MQ=255)            cTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACaaa                                                         <  2:395203/81‑1 (MQ=255)               ttCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGa                                                      <  1:395204/81‑1 (MQ=255)                  cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAAc                                                   <  1:395205/81‑1 (MQ=255)                     cTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAg                                                <  2:395206/81‑1 (MQ=255)                        tGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGaa                                             <  2:395207/81‑1 (MQ=255)                           aCCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAag                                           >  2:395228/1‑80 (MQ=255)                                     ggAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGAt                                >  1:395227/1‑81 (MQ=255)                                      gAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATg                               <  2:395208/81‑1 (MQ=255)                                                  gggTATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGtttt                   <  2:395209/81‑1 (MQ=255)                                                     tATCCCGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTg                >  1:395230/1‑81 (MQ=255)                                                        cccGCCGTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCg             >  2:395232/1‑81 (MQ=255)                                                          cgccgTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGtt           <  1:395210/81‑1 (MQ=255)                                                            ccgTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGc         <  1:395211/81‑1 (MQ=255)                                                             cgTCACACCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGcc        <  1:395212/81‑1 (MQ=255)                                                                cacaCCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGtt     >  2:395237/1‑81 (MQ=255)                                                                 acacCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGtt     <  2:395213/80‑1 (MQ=255)                                                                   acCCGAGCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGttttt  <  1:395214/81‑1 (MQ=255)                                                                         |                                                                           GCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTT  >  NC_000913/2233489‑2233635

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