Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 464,946 A→G 10.0% F531L (TTC→CTC)  ybaE ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048464,9460AG10.0% 130.0 / 7.2 50F531L (TTC→CTC) ybaEpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (2/3);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTC  >  AP009048/464886‑465015
                                                            |                                                                     
gcgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                                                                       >  1:3948246/1‑61 (MQ=255)
   gtgtCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGt                                                              <  1:1448725/67‑1 (MQ=255)
     gtCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                                                                       >  1:4250924/1‑56 (MQ=255)
        aGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCg                                                      >  1:2191488/1‑70 (MQ=255)
          tCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCAt                                                   >  1:3427904/1‑71 (MQ=255)
          tCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCAt                                                   >  1:2215012/1‑71 (MQ=255)
           cGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATc                                                  >  1:1553797/1‑71 (MQ=255)
            gCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGGGGCGTCAGCGTCGCATc                                                  >  1:170120/1‑70 (MQ=255)
              acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACa                                                                     >  1:3018934/1‑49 (MQ=255)
              acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACa                                                                     >  1:1840639/1‑49 (MQ=255)
              acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGt                                               <  1:161303/71‑1 (MQ=255)
                acCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCa                                             >  1:4100418/1‑71 (MQ=255)
                acCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCa                                             >  1:4062482/1‑71 (MQ=255)
                 ccGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                                                                       >  1:127906/1‑44 (MQ=255)
                       aCGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACt                                       <  1:1392153/70‑1 (MQ=255)
                       aCGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACt                                       <  1:4147932/70‑1 (MQ=255)
                        cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGt                                               >  1:4132945/1‑61 (MQ=255)
                        cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGt                                               <  1:1042891/61‑1 (MQ=255)
                        cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTgg                                     >  1:3778025/1‑71 (MQ=255)
                        cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTgg                                     >  1:1567176/1‑71 (MQ=255)
                        cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTgg                                     >  1:3529347/1‑71 (MQ=255)
                          ccGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCt                                   <  1:2771105/71‑1 (MQ=255)
                          ccGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGc                                    >  1:2440603/1‑70 (MQ=255)
                          ccGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGc                                    >  1:3543150/1‑70 (MQ=255)
                          ccGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGGCATTAACTGGc                                    >  1:3642827/1‑70 (MQ=255)
                           cGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATc                                                  >  1:865471/1‑55 (MQ=255)
                           cGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCt                                   >  1:2928240/1‑70 (MQ=255)
                            gggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaa                                 <  1:1191473/71‑1 (MQ=255)
                            gggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaa                                 <  1:903520/71‑1 (MQ=255)
                             ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                <  1:2168975/71‑1 (MQ=255)
                             ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                <  1:1391227/71‑1 (MQ=255)
                             ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaa                                 >  1:2279010/1‑70 (MQ=255)
                                    cACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                <  1:982544/64‑1 (MQ=255)
                                     aCTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                              >  1:1368337/1‑65 (MQ=255)
                                       tAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCa                                                    <  1:1538343/41‑1 (MQ=255)
                                        aaTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCAtt                                           <  1:4174503/49‑1 (MQ=255)
                                        aaTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCAtt                                           <  1:3116040/49‑1 (MQ=255)
                                         aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAg                    >  1:4556739/1‑71 (MQ=255)
                                           gCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACt                                       <  1:461070/50‑1 (MQ=255)
                                           gCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCtt                        <  1:1097777/65‑1 (MQ=255)
                                            cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCa                     >  1:3597203/1‑67 (MQ=255)
                                              atagtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCAt               <  1:2615810/71‑1 (MQ=255)
                                              atagtgatagtTGAACAGCGGCGGCAGCGTCGCATCTg                                                <  1:4047125/38‑1 (MQ=255)
                                                agtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCAtt                                           >  1:2339126/1‑41 (MQ=255)
                                                 gtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTa             >  1:3739687/1‑70 (MQ=255)
                                                  tgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTc                      <  1:2142985/60‑1 (MQ=255)
                                                   gatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTaa            >  1:3831331/1‑69 (MQ=255)
                                                    atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTaaa           <  1:1969776/69‑1 (MQ=255)
                                                    atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATc         <  1:4346635/71‑1 (MQ=255)
                                                     tagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCg        >  1:1822023/1‑71 (MQ=255)
                                                        ttGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACt                                       <  1:1054382/37‑1 (MQ=255)
                                                         tGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGtt     <  1:100247/70‑1 (MQ=255)
                                                          gAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGttt    <  1:1167103/70‑1 (MQ=255)
                                                           aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttc  >  1:169419/1‑71 (MQ=255)
                                                            |                                                                     
GCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTC  >  AP009048/464886‑465015

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: