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,5890TA94.3% 92.0 / ‑1.2 35K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (12/21);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCG  >  W3110S.gb/4563524‑4563653
                                                                 |                                                                
tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                             <  1:2247682/71‑1 (MQ=255)
tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                             <  1:538873/71‑1 (MQ=255)
  ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                            <  1:807386/70‑1 (MQ=255)
      gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                         <  1:3000415/69‑1 (MQ=255)
      gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                         <  1:18947/69‑1 (MQ=255)
      gcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                        >  1:3346412/1‑70 (MQ=255)
       cagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                      <  1:1934241/71‑1 (MQ=255)
        agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                     >  1:700781/1‑71 (MQ=255)
            tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                 >  1:2999134/1‑71 (MQ=255)
             gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                 <  1:2097458/70‑1 (MQ=255)
             gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                 <  1:2414615/70‑1 (MQ=255)
             gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                 <  1:3796678/70‑1 (MQ=255)
                         aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                          <  1:883089/65‑1 (MQ=255)
                         aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                    >  1:2486036/1‑71 (MQ=255)
                         aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                    <  1:2555496/71‑1 (MQ=255)
                             cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                <  1:1470115/71‑1 (MQ=255)
                               gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                              <  1:946538/71‑1 (MQ=255)
                               gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                              <  1:131075/71‑1 (MQ=255)
                                ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                     <  1:116231/63‑1 (MQ=255)
                                 tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                            >  1:3685210/1‑71 (MQ=255)
                                  ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                           >  1:1840567/1‑71 (MQ=255)
                                          ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                   <  1:3067537/71‑1 (MQ=255)
                                          ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                   <  1:32732/71‑1 (MQ=255)
                                             cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                               >  1:2209547/1‑40 (MQ=255)
                                             cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                >  1:3225368/1‑39 (MQ=255)
                                              aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                >  1:2736154/1‑70 (MQ=255)
                                              aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                >  1:2461929/1‑70 (MQ=255)
                                                 gacaGGCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc            <  1:3432970/68‑1 (MQ=255)
                                                 gacaGGCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc            <  1:2932457/68‑1 (MQ=255)
                                                 aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc            >  1:2558964/1‑71 (MQ=255)
                                                   tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                <  1:207976/65‑1 (MQ=255)
                                                   tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                <  1:3207145/65‑1 (MQ=255)
                                                          tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc   <  1:3586822/71‑1 (MQ=255)
                                                           aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg  >  1:3762063/1‑71 (MQ=255)
                                                             cGCCATTTTACTGTCGCCGATGTGGAACACGGTGACGCCAg                              <  1:1900422/41‑1 (MQ=38)
                                                                 |                                                                
TCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCG  >  W3110S.gb/4563524‑4563653

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: