Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,563,589 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
*W3110S.gb4,563,5890TA95.6% 100.1 / ‑5.4 36K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (14/20);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654
                                                                   |                                                                 
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2126652/71‑1 (MQ=255)
  tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:2484687/71‑1 (MQ=255)
    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            <  1:2253625/71‑1 (MQ=255)
    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            <  1:631252/71‑1 (MQ=255)
         cagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1367518/63‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   <  1:706203/69‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   <  1:1304334/69‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   <  1:1771361/69‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   <  1:1728372/69‑1 (MQ=255)
                cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 <  1:2436069/70‑1 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:1044200/1‑71 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:1024239/70‑1 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:1988151/70‑1 (MQ=255)
                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:2275255/70‑1 (MQ=255)
                           aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                        <  1:2295708/52‑1 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatt                                                                 >  1:1018278/1‑40 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatt                                                                 >  1:1181180/1‑40 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            >  1:768467/1‑45 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:608525/1‑71 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:1826603/1‑71 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:885926/1‑71 (MQ=255)
                                 gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1429336/50‑1 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:246899/58‑1 (MQ=255)
                                       cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:1695016/70‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:101829/38‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1050559/38‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:225917/1‑39 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                         >  1:1405207/1‑63 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:1533905/37‑1 (MQ=255)
                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1489097/47‑1 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                >  1:1198263/1‑52 (MQ=255)
                                                     tAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:326126/49‑1 (MQ=255)
                                                           gTAACGCCTTTTTACTGTCGCCCATGTGTAACACGGTGACGCCAGGtt                            <  1:90074/48‑1 (MQ=255)
                                                           gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1297308/37‑1 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           >  1:1904189/1‑46 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          >  1:2066834/1‑47 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          >  1:1432404/1‑47 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            >  1:2354304/1‑61 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:1426065/1‑71 (MQ=255)
                                                                gCCATTTTACTGTCGCCCAGGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              >  1:1666608/1‑57 (MQ=255)
                                                                   |                                                                 
GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654

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