Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 241,172 C→T 15.4% G711D (GGT→GAT)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048241,1720CT15.4% 134.1 / 15.2 52G711D (GGT→GAT) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base T (6/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCTTGTC  >  AP009048/241103‑241242
                                                                     |                                                                      
gcAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                       >  1:1275095/1‑71 (MQ=255)
  aaTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATcgc                                                                     <  1:4341535/71‑1 (MQ=255)
  aaTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                     <  1:403139/71‑1 (MQ=255)
    tCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGa                                                                   >  1:2183975/1‑71 (MQ=255)
      acaTCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGAcc                                                                 <  1:1022474/71‑1 (MQ=255)
        aTCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                       <  1:4198598/63‑1 (MQ=255)
         tCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAAt                                                              <  1:3036567/71‑1 (MQ=255)
         tCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAAt                                                              >  1:2961293/1‑71 (MQ=255)
          ccaccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                       >  1:260375/1‑61 (MQ=255)
          ccaccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATg                                                             >  1:2103973/1‑71 (MQ=255)
           caccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACGAATGc                                                            <  1:1691889/71‑1 (MQ=255)
             ccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCTATGCgg                                                          <  1:601000/71‑1 (MQ=255)
              caGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                       >  1:2222463/1‑57 (MQ=255)
              caGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCgg                                                          >  1:4002690/1‑70 (MQ=255)
               aGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCaa                                                               >  1:2821696/1‑64 (MQ=255)
                gcgcCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGaa                                                       >  1:3878250/1‑71 (MQ=255)
                    ctctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACggg                                                   <  1:760761/71‑1 (MQ=255)
                    ctctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACggg                                                   <  1:2252072/71‑1 (MQ=255)
                     tctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGaa                                                       >  1:3747475/1‑66 (MQ=255)
                       ttCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                                  <  1:580434/69‑1 (MQ=255)
                       ttCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTgg                                                >  1:4161398/1‑71 (MQ=255)
                       ttCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTgg                                                >  1:3534928/1‑71 (MQ=255)
                         ccAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAAc                                                      >  1:1131906/1‑63 (MQ=255)
                            gCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCaa                                                               >  1:1140078/1‑51 (MQ=255)
                            gCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTc                                           >  1:2212768/1‑71 (MQ=255)
                            gCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTc                                           >  1:1910955/1‑71 (MQ=255)
                            gCAAGCCAACCGGATTATGCTCGCGCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGtt                                            <  1:2508483/70‑1 (MQ=255)
                             cAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCg                                          >  1:2549212/1‑71 (MQ=255)
                              aaGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAAt                                                              >  1:1211153/1‑50 (MQ=255)
                               aGCGAACCGGATTATGCTCGCTCGGGGTCAGGTACTGACCGCGACCAATGCGGGaa                                                       >  1:1308094/1‑56 (MQ=255)
                                gCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCa                                       >  1:853995/1‑71 (MQ=255)
                                 ccAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                       >  1:1611396/1‑38 (MQ=255)
                                   aaCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACtt                                    <  1:2624839/71‑1 (MQ=255)
                                   aaCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACtt                                    <  1:2630377/71‑1 (MQ=255)
                                    aCCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTg                                   >  1:1637039/1‑71 (MQ=255)
                                    aCCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTg                                   >  1:1642476/1‑71 (MQ=255)
                                    aCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTg                                   >  1:9117/1‑71 (MQ=255)
                                     ccGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTg                                   >  1:1582464/1‑70 (MQ=255)
                                      cGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTa                                 <  1:3836192/71‑1 (MQ=255)
                                        gATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTaa                                >  1:1971847/1‑70 (MQ=255)
                                           tATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATc                            >  1:4458834/1‑71 (MQ=255)
                                           tATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATc                            >  1:1293496/1‑71 (MQ=255)
                                             ttctcgctcgGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCtt                          <  1:87493/69‑1 (MQ=255)
                                               ctcgctcgGCGTCAGGTACTGACCGCGACCAATGCGGGAACgg                                                    >  1:3949595/1‑43 (MQ=255)
                                               ctcgctcgGCGTCAGGTACTGACCGCGACCAATGCGGGAACg                                                     <  1:109438/42‑1 (MQ=255)
                                                tcgctcgGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGt                                                  >  1:248527/1‑44 (MQ=255)
                                                     cgGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTaaa                               <  1:3983107/58‑1 (MQ=255)
                                                      ggCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACtt                  >  1:400376/1‑70 (MQ=255)
                                                       gCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTa                                 <  1:798817/54‑1 (MQ=255)
                                                          tCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATg              <  1:2752804/70‑1 (MQ=255)
                                                          tCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGa             >  1:177602/1‑71 (MQ=255)
                                                             ggTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATcc          <  1:1719165/71‑1 (MQ=255)
                                                              gTACTGACCGCGACCAATGCGGGAACGGGTGGCGt                                             <  1:963341/35‑1 (MQ=255)
                                                                aCTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACtt                                    >  1:2346220/1‑42 (MQ=255)
                                                                  tGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCtt     >  1:81313/1‑71 (MQ=255)
                                                                   gACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCTTg    >  1:4449904/1‑71 (MQ=255)
                                                                    aCCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCt                           <  1:2916210/47‑1 (MQ=255)
                                                                     ccGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCTTGt   >  1:3835641/1‑70 (MQ=255)
                                                                     ccGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCTTGTc  <  1:4331572/71‑1 (MQ=255)
                                                                     |                                                                      
GCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGCTTGTC  >  AP009048/241103‑241242

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: