Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,187,308	T→C	P183P (CCA→CCG)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,187,308	0	T	C	100.0%	176.6 / NA	52	P183P (CCA→CCG)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/30);  total (22/30)

GTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT  >  minE/2187245‑2187372                                                                |                                                                 gtgtTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaaaa                                                           <  1:723565/71‑1 (MQ=255)    tttCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATg                                                         <  1:489000/70‑1 (MQ=255)    tttCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATg                                                         <  1:1861208/70‑1 (MQ=255)       ctctTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCg                                                      <  1:780301/70‑1 (MQ=255)       ctctTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGa                                                     >  1:361641/1‑71 (MQ=255)        tctTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaaaa                                                           >  1:1533783/1‑64 (MQ=255)          tttCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGAt                                                    <  1:758184/69‑1 (MQ=255)           ttCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaa                                                             <  1:939504/59‑1 (MQ=255)            tCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATa                                                   >  1:923615/1‑68 (MQ=255)            tCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATa                                                   <  1:1059045/68‑1 (MQ=255)             ccAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATaa                                                  <  1:551089/68‑1 (MQ=255)              cAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGaa                                              <  1:173715/71‑1 (MQ=255)                gTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGcc                                                       <  1:1386211/60‑1 (MQ=255)                gTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATaaa                                                 <  1:167670/66‑1 (MQ=255)                 tCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAg                                           >  1:1540689/1‑71 (MQ=255)                 tCTACCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGc                                                        >  1:143193/1‑58 (MQ=255)                   ttccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAg                                                <  1:2447525/64‑1 (MQ=255)                     ccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAgc                                          >  1:1398444/1‑68 (MQ=255)                     ccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCa                                       >  1:1534518/1‑71 (MQ=255)                      cgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGaa                                              >  1:48249/1‑63 (MQ=255)                      cgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGaa                                              >  1:3035157/1‑63 (MQ=255)                         ttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGc                                    >  1:3049974/1‑70 (MQ=255)                              ggATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaa                                                             <  1:293749/40‑1 (MQ=255)                                aTGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGccc                              >  1:449055/1‑69 (MQ=255)                                 tGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGaa                                              <  1:768181/52‑1 (MQ=255)                                 tGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa                            <  1:1558896/70‑1 (MQ=255)                                  gTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAg                           <  1:1172555/70‑1 (MQ=255)                                   tttttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCg                                                      <  1:2559677/42‑1 (MQ=255)                                   tttttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAgc                                          <  1:1730159/54‑1 (MQ=255)                                     tttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCg                       >  1:1015661/1‑71 (MQ=255)                                      ttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCg                       <  1:7820/70‑1 (MQ=255)                                      ttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCg                       <  1:409522/70‑1 (MQ=255)                                         cGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGAcc                        <  1:3026810/66‑1 (MQ=255)                                          gCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAgc                                          >  1:3057149/1‑47 (MQ=255)                                           cAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAacca                 >  1:1033708/1‑71 (MQ=255)                                           cAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAacca                 >  1:197257/1‑71 (MQ=255)                                            aGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccac                <  1:2339212/71‑1 (MQ=255)                                               cAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAAc                                             >  1:2002626/1‑39 (MQ=255)                                                 gagagCAAAGATGATCGGCAAAATGCCGATAAAGAACAgc                                          <  1:1008547/40‑1 (MQ=255)                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGc                                    >  1:2073430/1‑44 (MQ=255)                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGc                                    >  1:1949524/1‑44 (MQ=255)                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta         <  1:490474/71‑1 (MQ=255)                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta         >  1:1760615/1‑71 (MQ=255)                                                   gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta         <  1:2112096/71‑1 (MQ=255)                                                     gCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCtata       <  1:1069630/71‑1 (MQ=255)                                                      cAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccacc               >  1:2075632/1‑62 (MQ=255)                                                       aaaGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGAc                         <  1:1612598/51‑1 (MQ=255)                                                        aaGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAc    <  1:2035820/71‑1 (MQ=255)                                                         aGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCACAGAcc                        <  1:3871/50‑1 (MQ=38)                                                          gatgatCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGa                          >  1:1177806/1‑47 (MQ=255)                                                          gatgatCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCt  <  1:2451354/71‑1 (MQ=255)                                                            tgatCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCt  >  1:189699/1‑69 (MQ=255)                                                                |                                                                 GTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT  >  minE/2187245‑2187372

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