Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP028101	761,062	A→G	K479K (AAA→AAG)	C7Y58_03755 →	hemolysin activator protein precursor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP028101	761,062	0	A	G	100.0%	124.6 / NA	33	K479K (AAA→AAG)	C7Y58_03755	hemolysin activator protein precursor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

TGGAGTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATAGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAATTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAAAATTAAATTTTAGAT  >  CP028101/760957‑761295                                                                                                          |                                                                                                                                                                                                                                          tGGAGTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCt                                                                                           <  1:269326/250‑1 (MQ=37)   gAGTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATttagatttagat                                                                                                 <  2:293235/242‑1 (MQ=21)   gAGTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATttagatttagat                                                                                                 >  1:293235/1‑242 (MQ=21)     gTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATg                                                                                       >  2:132044/1‑250 (MQ=255)            gtgtTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAAc                                                                                 >  1:463422/1‑249 (MQ=255)                 ggtggtTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAg                                                                           >  1:190185/1‑250 (MQ=255)                   tggtTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCt                                                                         <  1:679675/250‑1 (MQ=255)                         tCATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAGAATATAGAAAAGGATATATGTCAGGTGCTATATTAGCTTTAAGATATAATATGAACTATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCaa                                                                  >  1:602606/1‑251 (MQ=255)                          cATAGAACTGGAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGGAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTTGTTCaaa                                                                 <  1:698477/251‑1 (MQ=255)                                   ggAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATAAAATAGTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTaaa                                                        >  1:502284/1‑251 (MQ=255)                                    gAAATATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTaaaa                                                       >  2:567565/1‑251 (MQ=255)                                     aaaTATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAAc                                                      <  1:132044/251‑1 (MQ=255)                                     aaaTATACAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGAAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGGGCTATATTATGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAAc                                                      <  2:707728/251‑1 (MQ=255)                                        tataCAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTa                                                   >  2:260666/1‑251 (MQ=255)                                         ataCAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTaa                                                  >  2:461198/1‑251 (MQ=255)                                          taCAAGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAg                                                 <  2:302586/251‑1 (MQ=255)                                             aaGGTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAAt                                              <  1:341969/251‑1 (MQ=255)                                                gTGATAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATaga                                           >  1:362884/1‑251 (MQ=255)                                                   aTAAAGCAATGGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGGAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAgaga                                         <  2:350425/250‑1 (MQ=255)                                                            ttgaaatagaaaATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATAtttt                               <  2:518647/249‑1 (MQ=255)                                                            tggaaatagaaaATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATAtttt                               <  2:573483/249‑1 (MQ=255)                                                                  tagaaaATGAATTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCa                         >  2:400874/1‑251 (MQ=255)                                                                            aTTATCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTaaaaa               <  1:150530/251‑1 (MQ=255)                                                                               aacATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACATTTAGCTCGTTCAATTTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAaaatta            >  2:37269/3‑251 (MQ=255)                                                                               aTCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAaaatta            >  1:84758/1‑251 (MQ=255)                                                                               aTCATATAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAaaatta            >  1:158564/1‑251 (MQ=255)                                                                                  atatAGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAaaattaaat         >  2:622401/1‑251 (MQ=255)                                                                                    ataGAGTGTTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAAAATTAAAttt       >  1:30252/1‑251 (MQ=255)                                                                                          tgtTAGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAAAATTAAATTTTAGAt  <  2:463422/250‑1 (MQ=255)                                                                                              aGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAAc                                                      >  2:190925/1‑194 (MQ=255)                                                                                              aGACTCAGAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAAc                                                      <  1:190925/194‑1 (MQ=255)                                                                                                     gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATAtttt                               <  1:66310/210‑1 (MQ=255)                                                                                                     gAAAAGTTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATAtttt                               >  2:66310/1‑210 (MQ=255)                                                                                                          |                                                                                                                                                                                                                                          TGGAGTAGGAAGTGTTGGTGGTTATCATAGAACTGGAAATATACAAGGTGATAAAGCAATAGAAATAGAAAATGAATTATCATATAGAGTGTTAGACTCAGAAAAATTTGGAAAAATAACTCCTTATTTAAGTTATTCATATGGAAAAGTAAGAAATAACAAAAATAATTCAAAATATAGAAAAGGATATATGTCAGGTGCTATATTAGGTTTAAGATATAATATGAAATATTTAGATTTAGATGTAGCTTATGCAAAACCTTTAGCTCGTTCAAATTACTTAAAACCTAAGAATAGAGAAATATATTTTAGTGCAACATTAAAAATTAAATTTTAGAT  >  CP028101/760957‑761295

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