Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,792 A→G 7.5% L354L (TTA→CTA) ‡ tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,7920AG 7.5% 136.4 / 2.6 53L354L (TTA→CTA) ‡tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (30/19);  new base G (2/2);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

AGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGG  >  minE/2093722‑2093860
                                                                      |                                                                    
aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                                      <  1:3242207/71‑1 (MQ=255)
  gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                                    <  1:829348/71‑1 (MQ=255)
  gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                                    <  1:4757352/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)
               gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                                    >  1:1375060/1‑58 (MQ=255)
               gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                                    >  1:5362224/1‑58 (MQ=255)
               gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                                    >  1:2489800/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:4668797/70‑1 (MQ=255)
                                gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                       <  1:4730031/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:384699/70‑1 (MQ=255)
                                                   ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa                    >  1:2655589/1‑70 (MQ=255)
                                                       ccAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTcc               <  1:995537/71‑1 (MQ=255)
                                                        cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                         >  1:552915/1‑44 (MQ=255)
                                                        cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                         >  1:3114831/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:756109/1‑43 (MQ=255)
                                                             gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                     >  1:5150474/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:5442082/1‑68 (MQ=255)
                                                                cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa         >  1:1812151/1‑68 (MQ=255)
                                                                cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa         >  1:4780405/1‑68 (MQ=255)
                                                                cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa         >  1:4659274/1‑68 (MQ=255)
                                                                     aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa         >  1:2918784/1‑63 (MQ=255)
                                                                     aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAgaga     >  1:3914397/1‑67 (MQ=255)
                                                                     aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAgaga     >  1:4206147/1‑67 (MQ=255)
                                                                     aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATgg  <  1:2754508/70‑1 (MQ=255)
                                                                     aaCCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATg   >  1:2734480/1‑69 (MQ=255)
                                                                      |                                                                    
AGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGAGATGG  >  minE/2093722‑2093860

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: