Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	4,563,587	T→A	100%	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	4,563,587	0	T	A	100.0%	178.2 / NA	59	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/29);  total (30/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  AP009048/4563520‑4563652                                                                    |                                                                  ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:199030/71‑1 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:4266933/71‑1 (MQ=255)  gtaggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1426843/68‑1 (MQ=255)  gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:3684262/71‑1 (MQ=255)   tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:2222115/71‑1 (MQ=255)     ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:2908896/70‑1 (MQ=255)       cagcagCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                          >  1:4494414/1‑71 (MQ=255)            gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                      >  1:2838317/1‑70 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                     >  1:676889/1‑70 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:865899/1‑71 (MQ=255)                 cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:891339/55‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1971273/53‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                >  1:2789754/1‑68 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                >  1:3148754/1‑68 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             <  1:2836812/71‑1 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:3117901/1‑71 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:327546/1‑71 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:862796/1‑71 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:3180743/1‑71 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               >  1:1262044/1‑66 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:2058666/1‑71 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1752898/71‑1 (MQ=255)                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:3635989/1‑63 (MQ=255)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                      >  1:692819/1‑53 (MQ=255)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGTTGAc                                    <  1:2211458/71‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                <  1:3443822/57‑1 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:3243068/1‑70 (MQ=255)                                 aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                <  1:1672440/71‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:2136953/40‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:1018004/40‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                <  1:3263267/53‑1 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                        >  1:339434/1‑59 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:97542/1‑68 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:2257499/1‑71 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:1523504/70‑1 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:4046253/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:2883340/70‑1 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:2487251/70‑1 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                         <  1:2377061/71‑1 (MQ=255)                                           cagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgc                      >  1:4353286/1‑71 (MQ=255)                                            agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:4013018/71‑1 (MQ=255)                                             ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1349257/39‑1 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   >  1:4178144/1‑37 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:793561/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 >  1:2930223/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                <  1:1259729/71‑1 (MQ=255)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               <  1:3269300/54‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                          >  1:2843663/1‑42 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                     >  1:1034177/1‑47 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGCCGCCAGGTTTACCGCCGAGCAAAc             >  1:2794532/1‑71 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              <  1:1335121/70‑1 (MQ=255)                                                    aaTAGGCGGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             <  1:2484094/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCTCCGAGCAAACCGcc         <  1:1299799/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:2623948/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:867430/71‑1 (MQ=255)                                                           tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:3437929/50‑1 (MQ=255)                                                           tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAc      >  1:4282352/1‑71 (MQ=255)                                                            gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:3906858/67‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                     >  1:3439726/1‑38 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                              <  1:4115695/45‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:1582778/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:2410313/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:3409162/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             <  1:2489404/60‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:1160938/71‑1 (MQ=255)                                                                    |                                                                  GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  AP009048/4563520‑4563652

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