Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,801,836 G→A 51.5% pseudogene (60/882 nt) ygaY → ECK2675:JW5428:b2681; predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,801,8360GA51.5% ‑3.7 / 59.0 66pseudogene (60/882 nt)ygaYECK2675:JW5428:b2681; predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (25/7);  new base A (12/22);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 5.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATG  >  minE/1801781‑1801893
                                                       |                                                         
tCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGc                                              >  1:1795932/1‑69 (MQ=255)
tCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGc                                              >  1:1792003/1‑69 (MQ=255)
 ccTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGc                                              >  1:2430538/1‑68 (MQ=255)
  cTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTTTGGTTCCGCTGGCAGCGACGCTg                                            <  1:2239142/69‑1 (MQ=255)
  cTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTTTGGTTCCGCTGGCAGCGACGCTg                                            <  1:2236341/69‑1 (MQ=255)
  cTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTTTGGTTCCGCTGGCAGCGACGCTg                                            <  1:3072587/69‑1 (MQ=255)
  cTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTTTGGTTCCGCTGGCAGCGACGCTg                                            <  1:1723951/69‑1 (MQ=255)
  cTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTTTGGTTCCGCTGGCAGCGACGCTg                                            <  1:3192217/69‑1 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTTGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:92693/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGc                                                          >  1:2617511/1‑43 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:1890403/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:2626726/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:2168024/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:3130967/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:3411765/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:1003479/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:474063/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:458035/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:1777268/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:3194572/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:1765622/1‑68 (MQ=255)
              ttAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:1772915/1‑68 (MQ=255)
              ttAACCGGTTTATTATCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGa                                 >  1:3055930/1‑68 (MQ=255)
                aaCCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATaa                              >  1:2942535/1‑69 (MQ=255)
                aaCCGGTTTATTCACAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATaa                              >  1:716373/1‑69 (MQ=255)
                        tATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCg                                   >  1:2212413/1‑56 (MQ=255)
                        tATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCg                                   >  1:3415434/1‑56 (MQ=255)
                        tATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCg                                   >  1:2903339/1‑56 (MQ=255)
                            ctcAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGt                  >  1:2128145/1‑69 (MQ=255)
                            ctcAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGt                  >  1:3175326/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:3256070/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:1185044/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:2867646/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:2772185/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:2488034/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:1663479/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:242108/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:1666536/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:1944569/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:2084965/1‑69 (MQ=255)
                             tcAGTCGTGGCACAAATTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTg                 >  1:1996968/1‑69 (MQ=255)
                             tcAGTCGTGGCAAAAATTCTGGTTCCACTGGAAGCGAc                                                >  1:3187530/1‑38 (MQ=37)
                                     ggCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCAc         <  1:1238583/69‑1 (MQ=255)
                                        acaAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACtat      <  1:1426788/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGGGGCAAAGTGGTTGGCACCATTATg  <  1:2539747/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2691172/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:936019/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2528653/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1518807/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:468152/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:370835/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:3244070/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:3220015/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:286221/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2411668/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1942151/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2248677/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGAACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:3366063/69‑1 (MQ=255)
                                            aTTCTGGTTCCACTGGCAGCGACGCTGGCCTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1536295/69‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:325080/68‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:325271/68‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1636301/68‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1532214/68‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:1855796/68‑1 (MQ=255)
                                             ttCTGGTTCCACTGGCAGCGACGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2434563/68‑1 (MQ=255)
                                              tCTGGTTCCACTGGCAGCGTCGCTGGCTTCACCGGACAAACGCGGCAAAGTGGTTGGCACCATTATg  <  1:2069092/67‑1 (MQ=255)
                                                       |                                                         
TCCTCGGTACGGCATTAACCGGTTTATTCTCAGTCGTGGCACAAATTCTGGTTCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATG  >  minE/1801781‑1801893

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: