Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,790 T→C 7.1% L354L (TTA→TTG) ‡ tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,7900TC 7.1% 159.3 / 2.5 56L354L (TTA→TTG) ‡tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

AAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGA  >  minE/2093721‑2093855
                                                                     |                                                                 
aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                                  <  1:1068627/71‑1 (MQ=255)
aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                                  <  1:1061863/71‑1 (MQ=255)
 aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                                 <  1:3242207/71‑1 (MQ=255)
   gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               <  1:4757352/71‑1 (MQ=255)
   gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               <  1:829348/71‑1 (MQ=255)
            cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                                                      >  1:5469700/1‑71 (MQ=255)
             aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     >  1:3856781/1‑71 (MQ=255)
              tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     >  1:2843138/1‑70 (MQ=255)
                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                                   >  1:682396/1‑54 (MQ=255)
                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:5362224/1‑58 (MQ=255)
                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:2489800/1‑58 (MQ=255)
                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:1375060/1‑58 (MQ=255)
                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                                  >  1:849638/1‑71 (MQ=255)
                  gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAtt                                                 >  1:4173286/1‑70 (MQ=255)
                  gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTg                                                <  1:4048256/71‑1 (MQ=255)
                    tttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     <  1:2904729/64‑1 (MQ=255)
                       ccgccgCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGagc                                             <  1:4323993/69‑3 (MQ=255)
                       ccgccgCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACa                                            <  1:161852/70‑1 (MQ=255)
                        cgccgcGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacg                                          >  1:4570746/1‑69 (MQ=255)
                             cATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCa                                                   >  1:2495466/1‑57 (MQ=255)
                               ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:4963297/1‑70 (MQ=255)
                                tGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCAcgccacgcca                                  <  1:543722/71‑1 (MQ=255)
                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAAc                                             >  1:1229003/1‑59 (MQ=255)
                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4730031/70‑1 (MQ=255)
                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4668797/70‑1 (MQ=255)
                                     aCCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCACGCCACGCCaactt                              >  1:3027797/1‑67 (MQ=255)
                                       cTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTccc                           >  1:4067880/1‑71 (MQ=255)
                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAtc                               >  1:5441428/1‑64 (MQ=255)
                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt                        <  1:4555033/71‑1 (MQ=255)
                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt                        <  1:756842/71‑1 (MQ=255)
                                            ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa                      <  1:344475/71‑1 (MQ=255)
                                            ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa                      <  1:663082/71‑1 (MQ=255)
                                              ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                                              >  1:3686379/1‑45 (MQ=255)
                                              ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACa                    <  1:1509830/71‑1 (MQ=255)
                                               gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAg                   >  1:258567/1‑71 (MQ=255)
                                                   gccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCg                >  1:2350007/1‑70 (MQ=255)
                                                    ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               >  1:2655589/1‑70 (MQ=255)
                                                    ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               <  1:384699/70‑1 (MQ=255)
                                                        ccAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTcc          <  1:995537/71‑1 (MQ=255)
                                                         cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:3114831/1‑44 (MQ=255)
                                                         cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:552915/1‑44 (MQ=255)
                                                             agtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4597969/42‑1 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:998413/1‑37 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:1670102/1‑37 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:3106350/1‑37 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                >  1:5150474/1‑43 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                >  1:756109/1‑43 (MQ=255)
                                                              gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCgcgc      >  1:3478369/1‑69 (MQ=255)
                                                              gtTCAGTTAACACTGCCGCCGCCATTGAACAGGCCAc                                      >  1:5085282/1‑37 (MQ=255)
                                                               ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               <  1:2961925/59‑1 (MQ=255)
                                                               ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTc           <  1:1246004/63‑1 (MQ=255)
                                                                tCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACa                    <  1:4327679/53‑1 (MQ=255)
                                                                 cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCaga  <  1:155535/70‑1 (MQ=255)
                                                                 cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:1812151/1‑68 (MQ=255)
                                                                 cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:5442082/1‑68 (MQ=255)
                                                                 cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:4780405/1‑68 (MQ=255)
                                                                 cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:4659274/1‑68 (MQ=255)
                                                                     |                                                                 
AAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGA  >  minE/2093721‑2093855

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