Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,670,644 A→G 20.0% I232T (ATT→ACT)  yfgL ← protein assembly complex, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,670,6440AG20.0% 81.0 / 17.1 40I232T (ATT→ACT) yfgLprotein assembly complex, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base A (12/20);  new base G (3/5);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCC  >  minE/1670574‑1670710
                                                                      |                                                                  
gTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGaaa                                                                    <  1:2932110/71‑1 (MQ=255)
              gAGTCGTGTCAACATCGCTCAGACGGTCAATTTTGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCaa                                                      <  1:2708139/71‑1 (MQ=255)
              gAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCaa                                                      <  1:31849/71‑1 (MQ=255)
                 tCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAatc                                                   <  1:582431/71‑1 (MQ=255)
                    tgtCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACgct                                                              <  1:3140006/57‑1 (MQ=255)
                      tCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTg                                              >  1:23938/1‑71 (MQ=255)
                          cATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGccc                                          <  1:3629419/71‑1 (MQ=255)
                           aTCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCtgctgc                                                         <  1:1910444/55‑1 (MQ=255)
                           aTCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCt                                               >  1:145307/1‑65 (MQ=255)
                            tCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGctgctg                                                          <  1:2450292/53‑1 (MQ=255)
                              gCTCAGACGGTCAATTTCGGTATAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGtt                                      >  1:856422/1‑71 (MQ=255)
                              gCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATAc                                                                 >  1:2382113/1‑44 (MQ=255)
                              gCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGtt                                      >  1:2625904/1‑71 (MQ=255)
                                 cAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCa                                   >  1:3561592/1‑71 (MQ=255)
                                 cAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTcc                                    <  1:2698825/70‑1 (MQ=255)
                                  aGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCTAAATCATCTGGCCCTGTTCCa                                   <  1:1665107/70‑1 (MQ=255)
                                  aGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTcc                                    <  1:1826131/69‑1 (MQ=255)
                                    aCGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAatca                                                  >  1:1449130/1‑53 (MQ=255)
                                       gTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAatcat                                                 <  1:618855/51‑1 (MQ=255)
                                        tCAATTTCGGTAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCAc                            >  1:1373405/1‑71 (MQ=255)
                                         cAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACt                           >  1:200169/1‑71 (MQ=255)
                                          aaTTTCGGTAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATc                                 >  1:1598695/1‑64 (MQ=255)
                                            tttCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCaa                                                      <  1:567899/41‑1 (MQ=255)
                                            tttCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAg                               <  1:3708827/64‑1 (MQ=255)
                                             ttCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCaaa                                                     >  1:478751/1‑41 (MQ=255)
                                                ggTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGa                    <  1:1570755/71‑1 (MQ=255)
                                                 gTAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAg                               <  1:327016/59‑1 (MQ=255)
                                                 gTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAatc                                                   <  1:805068/39‑1 (MQ=255)
                                                  tAGAACCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCa                                   <  1:672792/54‑1 (MQ=255)
                                                  tAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCt                      <  1:704833/67‑1 (MQ=255)
                                                      aCCGGTCGCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGAc              <  1:3078395/71‑1 (MQ=255)
                                                      aCCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGAc              <  1:2849015/71‑1 (MQ=255)
                                                         ggTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCAt           <  1:2813622/71‑1 (MQ=255)
                                                          gTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGa               >  1:3682756/1‑66 (MQ=255)
                                                             gCCTGGGAAGTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCt                      <  1:2599830/56‑1 (MQ=255)
                                                             gCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATc       >  1:1270500/1‑71 (MQ=255)
                                                              ccTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCg      <  1:1260442/71‑1 (MQ=255)
                                                                tGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTgc                         <  1:1913076/50‑1 (MQ=255)
                                                                 gggAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAg                               <  1:3654045/43‑1 (MQ=255)
                                                                  ggAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATGAGCACTGCGCTGACGCGACCATTATCGcccc  >  1:1580169/1‑71 (MQ=255)
                                                                  ggAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGcccc  >  1:2396935/1‑71 (MQ=255)
                                                                      |                                                                  
GTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCC  >  minE/1670574‑1670710

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: