Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,014,276	G→A	G286D (GGT→GAT)	ECB_02815 →	Capsule polysaccharide export inner‑membrane protein kpsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,014,276	0	G	A	100.0%	118.0 / NA	33	G286D (GGT→GAT)	ECB_02815	Capsule polysaccharide export inner‑membrane protein kpsE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/18);  total (15/18)

TGCGCGTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGGTGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGCTGT  >  REL606/3014205‑3014347                                                                        |                                                                        tGCGCGTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGAc                                                                     <  1:453860/76‑1 (MQ=255)  gcgcGTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACa                                                                    <  2:583740/76‑1 (MQ=255)  gcgcGTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACa                                                                    <  1:554065/76‑1 (MQ=255)      gTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCt                                                                >  1:851538/1‑76 (MQ=255)         aTGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTaaa                                                             <  2:1203730/76‑1 (MQ=255)          tGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGCCAAGCTAAAc                                                            <  1:1117261/76‑1 (MQ=255)            cGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCg                                                          <  1:895314/76‑1 (MQ=255)               ttCAGTCATTGCAGGCACAAATTGACGNAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTAt                                                       >  1:1204183/1‑76 (MQ=255)                tCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATg                                                      >  1:383037/1‑76 (MQ=255)                  aGTCAATGCAGGCACAAATTGACGCNGAAAAAAGCAAAATNACAGCGCTGCNGNATGACAAGCGAAACCGTATGGc                                                    >  2:20584/1‑76 (MQ=255)                   gTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCg                                                   <  1:423023/76‑1 (MQ=255)                        tGCAGGCACAAATTGACGAAGAAAAAAGAAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGa                                              >  2:599309/1‑76 (MQ=255)                         gCAGGCTCAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGAt                                             <  1:926369/76‑1 (MQ=255)                              cacaAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTga                                        <  2:346113/76‑1 (MQ=255)                              cacaAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGCCGGTGGATTTTga                                        <  1:86397/76‑1 (MQ=255)                               acaAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTgaa                                       <  2:49981/76‑1 (MQ=255)                                    ttGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGaaat                                  >  2:149948/1‑76 (MQ=255)                                    ttGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGaaat                                  >  1:184779/1‑76 (MQ=255)                                        cGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGacatcaaa                              >  2:909439/1‑76 (MQ=255)                                        cGAAGAAAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGaaatcaaa                              >  2:909434/1‑76 (MQ=255)                                          aagaaAAAAGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGaaatcaaatc                            <  2:52608/76‑1 (MQ=255)                                             aaaaaaGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAAtcaaatcaaa                         <  1:477088/76‑1 (MQ=255)                                                  aGCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTag                    <  1:236849/76‑1 (MQ=255)                                                   gCAAAATCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTaga                   <  1:389438/76‑1 (MQ=255)                                                     aaaaTCACAGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGt                 >  2:712579/1‑76 (MQ=255)                                                           acaGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAaca           >  2:479627/1‑76 (MQ=255)                                                            caGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAacac          >  2:415948/1‑76 (MQ=255)                                                             aGCGCCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACAcc         >  1:1161116/1‑76 (MQ=255)                                                              gcgcCGCAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCg        <  1:446450/76‑1 (MQ=255)                                                               cgccgcAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGa       <  2:505533/76‑1 (MQ=255)                                                               cgccgcAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGa       <  2:294982/76‑1 (MQ=255)                                                                 ccgcAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGc     <  1:82379/76‑1 (MQ=255)                                                                  cgcAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGCt    >  2:734485/1‑76 (MQ=255)                                                                   gcAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGCTg   >  2:471011/1‑76 (MQ=255)                                                                    cAGGATGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGCTGt  >  2:926884/1‑76 (MQ=255)                                                                        |                                                                        TGCGCGTAATGCGATTCAGTCATTGCAGGCACAAATTGACGAAGAAAAAAGCAAAATCACAGCGCCGCAGGGTGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAGCTGT  >  REL606/3014205‑3014347

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