Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,528,353	A→G	S189P (TCT→CCT)	sgcQ ←	KpLE2 phage‑like element; putative nucleoside triphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,528,353	0	A	G	100.0%	116.8 / NA	33	S189P (TCT→CCT)	sgcQ	KpLE2 phage‑like element; putative nucleoside triphosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/12);  total (21/12)

CGACGGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGAAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATCGCGATTGCCCAGATACA  >  NC_000913/4528278‑4528430                                                                            |                                                                              cgacgGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAacac                                                                          <  1:835000/81‑1 (MQ=255)  gacgGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAaca                                                                           >  1:835040/1‑79 (MQ=255)   acgGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACag                                                                        <  2:835001/81‑1 (MQ=255)    cgGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACaga                                                                       <  2:835002/81‑1 (MQ=255)     ggTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACagag                                                                      >  2:835049/1‑81 (MQ=255)         tCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAgcgc                                                                   <  2:835003/80‑1 (MQ=255)          cAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCAt                                                                 >  1:835045/1‑81 (MQ=255)           aGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCAt                                                                 <  2:835005/80‑1 (MQ=255)           aGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATc                                                                <  2:835004/81‑1 (MQ=255)             gCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAg                                                              >  1:835006/1‑81 (MQ=255)                ccGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGAt                                                           >  2:835038/1‑81 (MQ=255)                 cGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATg                                                          <  1:835007/81‑1 (MQ=255)                        ttGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCTAGGAACACAGAGCGCATCAGGATGgttgtt                                                    >  2:835046/1‑80 (MQ=255)                          gACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTaaa                                                 >  1:835036/1‑81 (MQ=255)                           aCCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTaaa                                                 >  2:835075/1‑80 (MQ=255)                           aCCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAAt                                                <  1:835008/81‑1 (MQ=255)                              cGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACt                                             <  2:835009/81‑1 (MQ=255)                              cGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACt                                             >  2:835059/1‑81 (MQ=255)                                 tttAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTc                                          >  2:835053/1‑81 (MQ=255)                                    aaCAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGAt                                       >  1:835047/1‑81 (MQ=255)                                     aCAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGAtt                                      <  1:835010/81‑1 (MQ=255)                                      cAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGAtt                                      >  2:835011/1‑80 (MQ=255)                                       aGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGAtt                                      >  1:835043/1‑79 (MQ=255)                                            gcTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCaa                                >  1:835065/1‑80 (MQ=255)                                            gcTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAAt                               <  1:835012/81‑1 (MQ=255)                                             cTGTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCaa                                >  2:835074/1‑79 (MQ=255)                                               gTCGGTACGCGTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAg                             >  1:835041/1‑80 (MQ=255)                                                      cgcgTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGata                     <  1:835013/81‑1 (MQ=255)                                                       gcgTCCCGG‑GGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAG‑AATAGAACAGATATcg                  >  2:835060/1‑81 (MQ=255)                                                         gTCCCGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATcg                  >  2:835051/1‑81 (MQ=255)                                                           cccGGCGGTCAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATcgc                 >  2:835061/1‑80 (MQ=255)                                                                    cAGGCCAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATCGCGATTGcc          >  2:835067/1‑78 (MQ=255)                                                                         cAGGAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATCGCGATTGCCCAGATaca  >  2:835057/1‑81 (MQ=255)                                                                            |                                                                              CGACGGTGTCAGGCACCGTTTCTTTGACCCGTTTTAACAGCGCGCTGTCGGTACGCGTCCCGGCGGTCAGGCCAGAAACACAGAGCGCATCAGGATGGTTGTTAAATACTGTCGATTTAGCAATAGAACAGATATCGCGATTGCCCAGATACA  >  NC_000913/4528278‑4528430

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: ‑