Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 981,683 Δ1 bp 33.3% coding (11/533 nt) ydfQ ← predicted lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE981,6830C.33.3% 50.4 / 23.0 24coding (11/533 nt)ydfQpredicted lysozyme
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base . (5/3);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCT  >  minE/981616‑981749
                                                                   |                                                                  
gTATCTGAGGAGCAGATGCTCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATtttt                                                                 <  1:3648653/71‑1 (MQ=255)
gTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATtttt                                                                 >  1:3251480/1‑71 (MQ=255)
  aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATttttg                                                                <  1:938423/70‑2 (MQ=255)
   tCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTtgt                                                              >  1:984255/1‑70 (MQ=255)
   tCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTtgt                                                              >  1:5252930/1‑70 (MQ=255)
     tGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGtt                                                           >  1:207669/1‑71 (MQ=255)
       aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGTTc                                                          <  1:4797198/70‑1 (MQ=255)
       aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGTTCa                                                         >  1:5373979/1‑71 (MQ=255)
         gagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATttttg                                                                >  1:2617870/1‑62 (MQ=255)
         gagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATttttg                                                                >  1:2996416/1‑62 (MQ=255)
                 gcgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTtgtg                                                             >  1:2555640/1‑57 (MQ=255)
                  cgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCt                                               >  1:1070278/1‑71 (MQ=255)
                     cAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGTTCa                                                         <  1:4332644/57‑1 (MQ=255)
                     cAGCACCAATCAGAGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGtt                                                           <  1:4041225/55‑1 (MQ=255)
                       gCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATttttg                                                                >  1:4892089/1‑48 (MQ=255)
                         aCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCt                                         >  1:1143069/1‑70 (MQ=255)
                            aaTCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTAc                                      >  1:3255114/1‑70 (MQ=255)
                            aaTCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTAc                                      >  1:5145841/1‑70 (MQ=255)
                              tCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTAc                                      <  1:3082781/68‑1 (MQ=255)
                                  cgcCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATtc                               <  1:297876/71‑1 (MQ=255)
                                   gcCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtg                                                               <  1:2838583/38‑1 (MQ=255)
                                        aaCGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCa                                                         >  1:892158/1‑39 (MQ=255)
                                          cGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACccgcccgc                                                >  1:2273123/1‑46 (MQ=255)
                                          cGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACccgcccgc                                                >  1:959034/1‑46 (MQ=255)
                                                ccGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTgg                 <  1:2635366/71‑1 (MQ=255)
                                                     aGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCa            <  1:1152428/71‑1 (MQ=255)
                                                     aGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCa            <  1:785446/71‑1 (MQ=255)
                                                       gCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGa                                    <  1:748868/45‑1 (MQ=255)
                                                       gCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGa                                    <  1:2446459/45‑1 (MQ=255)
                                                                 aTCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCt  <  1:1288189/69‑1 (MQ=255)
                                                                   cTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACt                     >  1:2735436/1‑48 (MQ=255)
                                                                   |                                                                  
GTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCT  >  minE/981616‑981749

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