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%	231.8 / NA	77	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (35/42);  total (35/42)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  AP009048/4563520‑4563654                                                                    |                                                                    ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCCGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:3643087/71‑1 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:2631229/71‑1 (MQ=255)   tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                <  1:1971851/71‑1 (MQ=255)   tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                <  1:4461443/71‑1 (MQ=255)        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                           >  1:5177399/1‑71 (MQ=255)         gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                          >  1:2383474/1‑71 (MQ=255)         gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                          >  1:4167065/1‑71 (MQ=255)         gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                          >  1:4077649/1‑71 (MQ=255)            gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                       <  1:3194855/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:3037691/59‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:450503/59‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:1221269/59‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                      <  1:2939037/71‑1 (MQ=255)               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    >  1:4153166/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   >  1:3574557/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   >  1:300591/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                     <  1:1057305/69‑1 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                     <  1:2379188/69‑1 (MQ=255)                  tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  >  1:795347/1‑54 (MQ=255)                   gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    >  1:334198/1‑67 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               >  1:4052421/1‑71 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                            <  1:4880249/71‑1 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             <  1:150460/70‑1 (MQ=255)                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           <  1:1009382/71‑1 (MQ=255)                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           >  1:1005436/1‑71 (MQ=255)                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           >  1:2335532/1‑71 (MQ=255)                         cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   <  1:4855251/62‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:2511247/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:4247327/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:2858026/71‑1 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        <  1:531847/71‑1 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        <  1:5140643/71‑1 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        <  1:3081156/71‑1 (MQ=255)                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       <  1:2562849/71‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCttttt                                                                 <  1:3854316/42‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                         >  1:3950149/1‑50 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                     >  1:1360463/1‑70 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    <  1:363568/71‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    >  1:450892/1‑71 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:1339571/71‑1 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                 <  1:2849208/39‑1 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 <  1:5134689/71‑1 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 >  1:4588450/1‑71 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 <  1:3510805/71‑1 (MQ=255)                                    tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         >  1:1562004/1‑61 (MQ=255)                                    tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                               >  1:64586/1‑71 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:2499251/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:4610734/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:4771816/70‑1 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                     >  1:4232700/1‑45 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                            <  1:324437/70‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               >  1:1040148/1‑49 (MQ=255)                                           cagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        >  1:3931585/1‑55 (MQ=255)                                            agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                       <  1:3947792/71‑1 (MQ=255)                                            agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                       <  1:4242545/71‑1 (MQ=255)                                             ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                      <  1:1703750/71‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               >  1:3444798/1‑44 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                    >  1:3991045/1‑70 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                    >  1:3985165/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         >  1:2122427/1‑48 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:1199523/39‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:4661408/39‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                   >  1:2155809/1‑69 (MQ=255)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                    >  1:2711534/1‑67 (MQ=255)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa                 <  1:1101391/70‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                <  1:1246004/70‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                <  1:4314630/70‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:2147335/1‑71 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc              <  1:1682069/71‑1 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc              <  1:1911065/71‑1 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                  <  1:1327177/49‑1 (MQ=255)                                                         gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAAGCCAGGTTTACCGCCGAGCAAACCGcc           <  1:3684238/70‑1 (MQ=255)                                                           tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAc        >  1:729734/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      >  1:4650710/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      >  1:4813746/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      <  1:505997/71‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      <  1:1216155/71‑1 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:2532583/1‑37 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg     >  1:328708/1‑71 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg     >  1:4545222/1‑71 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg     >  1:2720557/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    <  1:2726763/71‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    <  1:718003/71‑1 (MQ=255)                                                                 gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  <  1:2261127/71‑1 (MQ=255)                                                                 gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  >  1:993941/1‑71 (MQ=255)                                                                    |                                                                    GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  AP009048/4563520‑4563654

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: ‑