Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,028,040	A→T	L141* (TTA→TAA)	ECB_02828 ←	Polysialic acid transport protein kpsM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,028,040	0	A	T	100.0%	118.0 / NA	32	L141* (TTA→TAA)	ECB_02828	Polysialic acid transport protein kpsM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/18);  total (14/18)

CCAACGACCATGAATATTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTAAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCAACGT  >  REL606/3027968‑3028112                                                                         |                                                                         ccAACAACCATGAATATTAAGCCGACCCCACATNAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTtaaa                                                                       <  1:772396/76‑1 (MQ=255)  cAACGACCATGAATATTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTtaaaa                                                                      <  2:250605/76‑1 (MQ=255)        ccATGAATATTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGt                                                                >  1:490274/1‑76 (MQ=255)               tatTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCg                                                         >  1:412627/1‑76 (MQ=255)               tatTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGATATCTCg                                                         <  1:643367/76‑1 (MQ=255)                 tgaaGCGGACGCCACATGAAAGNATGATCAGCANGCTCCAGGNAGCCACAAGTCGTTAAAATTTTNTTATCTCGaa                                                       <  2:838528/74‑1 (MQ=255)                    aGCCGACGCCACATGAAAGATCGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAata                                                    <  2:810114/76‑1 (MQ=255)                     gCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAAGTTGTTATCTCGAAatat                                                   <  2:40117/76‑1 (MQ=255)                           gCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCt                                             <  1:876258/76‑1 (MQ=255)                            ccACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTg                                            >  2:96684/1‑76 (MQ=255)                                 tGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTACAAATTTGTTATCTCGAAATATTCGCCTGCCATc                                       >  1:175898/1‑76 (MQ=255)                                  gAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATcc                                      <  1:786273/76‑1 (MQ=255)                                  gAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATcc                                      <  1:935492/76‑1 (MQ=255)                                  gAAAGAATAATCAGCAGGCTCCAGGTAGCCACAAATTGTTAAAAATTTGTTATCTCGAAATATTCGCCCGCCATcc                                      <  2:865227/76‑1 (MQ=255)                                     aGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGa                                   <  1:250790/76‑1 (MQ=255)                                      gAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGAc                                  >  2:465483/1‑76 (MQ=255)                                        aTGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACaa                                <  1:518432/76‑1 (MQ=255)                                        aTGATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACaa                                <  1:934325/76‑1 (MQ=255)                                        aTGATCAGCACGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACaa                                >  1:126616/1‑76 (MQ=255)                                          gATCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATg                              <  1:295389/76‑1 (MQ=255)                                           aTCAGCAGGCTCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGa                             >  1:924476/1‑76 (MQ=255)                                                     tCCAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAg                   >  2:545990/1‑76 (MQ=255)                                                      ccAGGTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGa                  >  1:60491/1‑76 (MQ=255)                                                        aGGTAGCNACAAGTTGTTAAACCTNTGTTATCTCGAAATATNCGCCTGCCTTCCAGACAATGAGCATATGCAGAat                <  2:852655/76‑1 (MQ=255)                                                         ggTAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAata               <  2:77301/76‑1 (MQ=255)                                                           tAGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAAtata             >  1:321851/1‑76 (MQ=255)                                                            aGCCACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATaa            >  1:250052/1‑76 (MQ=255)                                                              ccACAAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAAc          <  2:220713/76‑1 (MQ=255)                                                              ccACAAGTTGTTAAAAATTTGTTATCTAGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAAc          <  1:581811/76‑1 (MQ=255)                                                                acaAGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAg        >  2:730821/1‑76 (MQ=255)                                                                  aaGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCa      >  2:334261/1‑76 (MQ=255)                                                                  aaGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCa      <  2:104427/76‑1 (MQ=255)                                                                  aaGTTGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCa      <  1:798022/76‑1 (MQ=255)                                                                      tGTTAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCAACGt  >  1:469370/1‑76 (MQ=255)                                                                         |                                                                         CCAACGACCATGAATATTAAGCCGACGCCACATGAAAGAATGATCAGCAGGCTCCAGGTAGCCACAAGTTGTAAAAAATTTGTTATCTCGAAATATTCGCCTGCCATCCAGACAATGAGCATAAGCAGAATATAAACAGCAACGT  >  REL606/3027968‑3028112

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