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%	199.8 / NA	66	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/36);  total (30/36)
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:3430441/71‑1 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2008176/71‑1 (MQ=255)  gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:4657627/71‑1 (MQ=255)    cGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:975977/1‑71 (MQ=255)     ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGTCTGTAACGCCATTTTACt                                                            <  1:3835234/71‑1 (MQ=255)      gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           <  1:2883089/71‑1 (MQ=255)      gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:4310863/1‑71 (MQ=255)      gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:5157758/1‑71 (MQ=255)        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                          >  1:4756294/1‑70 (MQ=255)        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                         >  1:1357193/1‑71 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                         <  1:1537302/66‑1 (MQ=255)             cTTGATGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:2483149/1‑71 (MQ=255)               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2509356/57‑1 (MQ=255)               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:4301366/1‑60 (MQ=255)               tGCTGATCGGCAAATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:5249447/1‑60 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:5280630/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:1890754/70‑1 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:3278744/70‑1 (MQ=255)                  tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                        >  1:5145468/1‑62 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:5032646/1‑57 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:4385363/1‑71 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:646252/1‑71 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:1213591/1‑71 (MQ=255)                     tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:3158308/1‑70 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:811347/1‑71 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCAATGTGGAaca                                           >  1:3005712/1‑71 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:3133407/49‑1 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:3304253/70‑1 (MQ=255)                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:4526513/48‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:2941755/71‑1 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:1388904/1‑70 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:1394972/1‑70 (MQ=255)                              cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   >  1:2797666/1‑71 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  <  1:2133707/71‑1 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   <  1:2808602/69‑1 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:101433/71‑1 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:3509158/71‑1 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               <  1:1053830/71‑1 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:3624935/70‑1 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:5157124/70‑1 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                >  1:2299113/1‑66 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:3157831/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:2156882/47‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:4401335/70‑1 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:2776586/70‑1 (MQ=255)                                         agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:387585/52‑1 (MQ=255)                                         agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAcc                        >  1:2406693/1‑71 (MQ=255)                                            agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:3083485/71‑1 (MQ=255)                                             ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                    <  1:4697869/71‑1 (MQ=255)                                              gTCATAAATAGTCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                   <  1:2753145/71‑1 (MQ=255)                                              gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:4257992/43‑1 (MQ=255)                                              gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                   <  1:2696687/71‑1 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:4559813/41‑1 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:5194684/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 <  1:4569205/70‑1 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 >  1:4261508/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                <  1:4121054/71‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:4092872/71‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:5265134/71‑1 (MQ=255)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              >  1:15087/1‑71 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                      >  1:3226949/1‑46 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            >  1:2096062/1‑71 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:2045263/71‑1 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:1208744/71‑1 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAACAAACCGc          >  1:2398400/1‑71 (MQ=255)                                                         gCTGTAACGCCATTTTACTGTCGCCCATGTGGAacact                                         <  1:4673831/38‑2 (MQ=38)                                                          cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:5129192/1‑47 (MQ=255)                                                            gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGACAAcg     >  1:2580784/1‑71 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:4157268/42‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcg     <  1:4289977/70‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:933043/71‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGACAAcgcg   <  1:1929918/70‑1 (MQ=255)                                                                    |                                                                  GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  AP009048/4563520‑4563652

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