Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,000 T→C 16.7% A284A (GCA→GCG)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,0000TC16.7% 128.3 / 19.1 54A284A (GCA→GCG) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (3/6);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGA  >  minE/2093931‑2094070
                                                                     |                                                                      
tACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTg                                                                       >  1:651059/1‑71 (MQ=255)
    ggACGACGGAACTTATCGTAATAAGCCAGCACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                   <  1:4042543/71‑1 (MQ=255)
       cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                                                >  1:1670121/1‑71 (MQ=255)
          cgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAc                                                             >  1:4040782/1‑71 (MQ=255)
           taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                             <  1:2549982/68‑1 (MQ=255)
               cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                        >  1:3169939/1‑71 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                     >  1:1998707/1‑71 (MQ=255)
                  aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                     >  1:3507484/1‑71 (MQ=255)
                       aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCa                                                   >  1:30859/1‑68 (MQ=255)
                       aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                >  1:2630179/1‑71 (MQ=255)
                        ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                >  1:868048/1‑70 (MQ=255)
                        ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTg                                               >  1:1189465/1‑71 (MQ=255)
                               aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                        <  1:262743/71‑1 (MQ=255)
                                 tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTTCGTTttctt                                      <  1:682388/71‑1 (MQ=255)
                                 tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                                             <  1:4121814/64‑1 (MQ=255)
                                 tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                      <  1:291785/71‑1 (MQ=255)
                                 tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                      <  1:4096712/71‑1 (MQ=255)
                                  aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                      >  1:224093/1‑70 (MQ=255)
                                   cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                    <  1:327164/71‑2 (MQ=255)
                                   cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                    <  1:2716984/71‑2 (MQ=255)
                                   cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGt                                                                        >  1:3280387/1‑35 (MQ=255)
                                   cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCa                                    >  1:2928687/1‑71 (MQ=255)
                                   cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACAACGCGCCCGCCAGCTGCGTTTTCTTCa                                    >  1:1124931/1‑71 (MQ=255)
                                     gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                   >  1:361529/1‑70 (MQ=255)
                                     gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                   >  1:3214810/1‑70 (MQ=255)
                                      ccGGAATAATGCCTCACGCACCGTTAGTCGGTGCCGTTACCa                                                              >  1:3215485/1‑42 (MQ=255)
                                      ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACt                                  >  1:1250371/1‑70 (MQ=255)
                                      ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                 <  1:139694/71‑1 (MQ=255)
                                       cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                <  1:2098593/71‑1 (MQ=255)
                                           ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                            >  1:3190474/1‑71 (MQ=255)
                                           ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                            >  1:3465786/1‑71 (MQ=255)
                                               tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAcc                              >  1:3695737/1‑65 (MQ=255)
                                               tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGTCAgc                        <  1:424534/71‑1 (MQ=255)
                                                 ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                       >  1:915853/1‑68 (MQ=255)
                                                 ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                       >  1:2807830/1‑68 (MQ=255)
                                                  cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                     <  1:3406803/71‑2 (MQ=255)
                                                  cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                     <  1:1148017/71‑2 (MQ=255)
                                                  cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgc                      >  1:3438212/1‑70 (MQ=255)
                                                   gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                     >  1:2686903/1‑70 (MQ=255)
                                                      cgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACa                 >  1:3285942/1‑71 (MQ=255)
                                                       gcacCGTTAGTCGGCGCAGGCACCACGCGACCGCCAGCAGCGTTctct                                       <  1:3950309/48‑1 (MQ=16)
                                                        cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg               <  1:520067/71‑1 (MQ=255)
                                                        cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg               <  1:2201729/71‑1 (MQ=255)
                                                         acCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGt              <  1:3827713/71‑1 (MQ=255)
                                                          ccGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                       <  1:4082012/45‑1 (MQ=255)
                                                             ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt          >  1:1401770/1‑71 (MQ=255)
                                                             ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt          <  1:2463316/71‑1 (MQ=255)
                                                              tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt          >  1:3042194/1‑70 (MQ=255)
                                                               aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                 <  1:2193282/46‑1 (MQ=255)
                                                                gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATcc        <  1:2943441/70‑1 (MQ=255)
                                                                gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa       <  1:3478794/71‑1 (MQ=255)
                                                                 tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGcc                           >  1:4110362/1‑50 (MQ=255)
                                                                 tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAg      <  1:2559309/71‑1 (MQ=255)
                                                                    gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt          <  1:2510980/64‑1 (MQ=255)
                                                                    gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCg   <  1:3659028/71‑1 (MQ=255)
                                                                     tGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGa  >  1:1423232/1‑71 (MQ=255)
                                                                     |                                                                      
TACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGA  >  minE/2093931‑2094070

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