Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,476,848 A→G 29.8% Y59Y (TAT→TAC inaA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,476,8480AG29.8% 157.9 / 80.4 164Y59Y (TAT→TACinaAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/77);  new base G (25/24);  total (63/101)
Fisher's exact test for biased strand distribution p-value = 3.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTT  >  minE/1476782‑1476899
                                                                  |                                                   
gACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGAt                                                    >  1:3127325/1‑68 (MQ=255)
gACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATa                                                   >  1:3162076/1‑69 (MQ=255)
   acCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACg                                                <  1:3099026/69‑1 (MQ=255)
   acCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACg                                                <  1:1288778/69‑1 (MQ=255)
        cccGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACgg                                           <  1:383781/69‑1 (MQ=255)
         ccGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGa                                          <  1:2454208/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:2263306/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:2615875/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:722479/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:1449878/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:2414003/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:1888500/69‑1 (MQ=255)
             tCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      <  1:395383/69‑1 (MQ=255)
              ccggTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:1465193/4‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATg                                       >  1:2675778/1‑67 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:2099692/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:3361548/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:3318917/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:3326917/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:551902/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:2418077/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:569429/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:2671850/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:1144766/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:1079498/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGa                                      >  1:1014871/1‑68 (MQ=255)
              ccAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGCCCGAACGGATAACGTACGGa                                          >  1:1956183/1‑64 (MQ=255)
              ccAGTTCTTTAATTACAGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAAc                                                 >  1:1213894/1‑57 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2874918/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2534793/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2548282/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:3270018/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2395734/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2383904/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:3272496/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:2275822/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:211396/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:849239/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:632142/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:1595978/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:609291/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:419298/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:1922779/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:1924772/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:3382582/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:3354258/63‑1 (MQ=255)
                 gTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAAt                                        <  1:333238/63‑1 (MQ=255)
                    cTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:237834/69‑1 (MQ=255)
                    cTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:1201741/69‑1 (MQ=255)
                    cTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:2631103/69‑1 (MQ=255)
                    cTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:3058064/69‑1 (MQ=255)
                    cTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:815436/69‑1 (MQ=255)
                     tttAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAgat                               <  1:1424201/68‑1 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTAc                                             >  1:2201023/1‑50 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:890617/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2923060/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2925852/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:674879/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:3119936/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:488046/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:360816/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:3429892/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2014976/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:1067713/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2753627/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:1181977/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2558846/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:1897697/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:1844466/1‑62 (MQ=255)
                         aTTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGGTAACGTACAGAATGAAACAg                                 >  1:2408989/1‑62 (MQ=255)
                           tACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:515495/69‑1 (MQ=255)
                           tACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:2025590/69‑1 (MQ=255)
                           tACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:2052899/69‑1 (MQ=255)
                           tACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:1869818/69‑1 (MQ=255)
                           tACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:1643023/69‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGaa                                     <  1:1016622/55‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGaa                                     <  1:1919081/55‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGaa                                     <  1:534412/55‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGaa                                     <  1:855674/55‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGaa                                     <  1:3231371/55‑1 (MQ=255)
                            aCCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGt                        <  1:2233755/68‑1 (MQ=255)
                               gCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATg                    <  1:2138206/69‑1 (MQ=255)
                               gCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATg                    <  1:2821991/69‑1 (MQ=255)
                               gCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATg                    <  1:1091680/69‑1 (MQ=255)
                               gCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATg                    <  1:1612426/69‑1 (MQ=255)
                                  aCCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATgcgc                 >  1:1118002/1‑69 (MQ=255)
                                                cGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATATAGCtt   <  1:1749626/69‑1 (MQ=255)
                                                cGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCtt   <  1:2470552/69‑1 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:2309578/60‑1 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:1724869/60‑1 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:1652050/60‑1 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCttt  >  1:2844090/1‑69 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCttt  >  1:3094033/1‑69 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCttt  >  1:1507196/1‑69 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCttt  >  1:493357/1‑69 (MQ=255)
                                                 gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCttt  >  1:394331/1‑69 (MQ=255)
                                                               cGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:2711979/46‑1 (MQ=255)
                                                               cGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:2800662/46‑1 (MQ=255)
                                                               cGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:1999717/46‑1 (MQ=255)
                                                                ggATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTtaca           <  1:264433/45‑1 (MQ=255)
                                                                  |                                                   
GACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTT  >  minE/1476782‑1476899

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: