Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,451,762 T→A I144F (ATC→TTC)  acpD ← acyl carrier protein phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,451,7620TA97.9% 166.2 / NA 48I144F (ATC→TTC) acpDacyl carrier protein phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (15/32);  minor base C (0/1);  total (15/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAG  >  REL606/1451689‑1451832
                                                                         |                                                                      
tCGGTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAAcc                                                                      <  1:1278012/76‑1 (MQ=255)
 cGGTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAAccc                                                                     >  1:258965/1‑76 (MQ=255)
 cGGTAATGCCGATAAAGCCGAGGAACGTGGACAGANACGGCGTCACCAGGTCCGTTGGNCCATCTTTGTGGAAccc                                                                     >  1:798665/1‑76 (MQ=255)
  ggTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccg                                                                    <  2:1129782/76‑1 (MQ=255)
   gTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgc                                                                   >  2:262146/1‑76 (MQ=255)
    tAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgcc                                                                  <  2:865519/76‑1 (MQ=255)
    tAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgcc                                                                  >  1:797888/1‑76 (MQ=255)
    tAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgcc                                                                  >  1:708778/1‑76 (MQ=255)
    tAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgcc                                                                  <  1:430547/76‑1 (MQ=255)
     aaTGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgccg                                                                 <  1:1015548/76‑1 (MQ=255)
     aaTGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgccg                                                                 <  2:1028007/76‑1 (MQ=255)
     aaTGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACccgccg                                                                 <  1:913178/76‑1 (MQ=255)
       tGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCcgcg                                                               <  1:1537269/76‑1 (MQ=255)
       tGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCcgcg                                                               <  1:88214/76‑1 (MQ=255)
          cGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCt                                                            <  2:1557575/76‑1 (MQ=255)
             gaaaGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGt                                                         <  2:778435/75‑1 (MQ=255)
              aaaGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGta                                                        <  2:343234/76‑1 (MQ=255)
              aaaGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGta                                                        <  2:169479/76‑1 (MQ=255)
              aaaGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGta                                                        <  2:1252209/76‑1 (MQ=255)
              aaaGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGta                                                        <  1:217222/76‑1 (MQ=255)
                aGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGtaat                                                      <  1:1184317/76‑1 (MQ=255)
                 gCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGtaata                                                     <  1:2811/76‑1 (MQ=255)
                       gAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATg                                               <  1:426178/76‑1 (MQ=255)
                               aCAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCtttttt                                       <  2:1356460/76‑1 (MQ=255)
                               aCAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCtttttt                                       <  1:253838/76‑1 (MQ=255)
                               aCAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCtttttt                                       <  2:1332198/76‑1 (MQ=255)
                               aCAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCtttttt                                       >  1:592485/1‑76 (MQ=255)
                                    tACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCg                                  <  2:436012/76‑1 (MQ=255)
                                     aCGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGt                                 >  2:221513/1‑76 (MQ=255)
                                     aCGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTAACCGt                                 >  2:392297/1‑76 (MQ=255)
                                      cGGCGTCCCCAGGTCCGTTGGTCCTTCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGtt                                <  1:418806/76‑1 (MQ=255)
                                      cGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGtt                                >  1:56948/1‑76 (MQ=255)
                                         cGTCCCCAGGTCCGTTGGTCCATCTTTGTGGACCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTAcc                             <  1:723883/76‑1 (MQ=255)
                                         cGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTAcc                             <  2:520165/76‑1 (MQ=255)
                                         cGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTAcc                             <  1:418795/76‑1 (MQ=255)
                                         cGTCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTAcc                             <  1:1010361/76‑1 (MQ=255)
                                           tCACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAg                           >  1:937063/1‑76 (MQ=255)
                                            cACCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGa                          <  1:465978/76‑1 (MQ=255)
                                             aCCAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGAc                         >  1:347924/1‑76 (MQ=255)
                                              ccAGGTCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGAcc                        >  2:1369461/1‑76 (MQ=255)
                                                  ntCCGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGACGGCTTTTTTACCCGTTACCAGACCATcc                    >  1:1579925/2‑76 (MQ=255)
                                                    ccGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCgg                  <  1:44251/76‑1 (MQ=255)
                                                    ccGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCgg                  >  2:245430/1‑76 (MQ=255)
                                                     cGTTGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGCTGGCTTTTTTACCCCTTACCAGACCTCCCGGa                 >  1:764110/1‑76 (MQ=255)
                                                        tGGTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCg              <  1:1320721/76‑1 (MQ=255)
                                                         ggTCCATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGt             >  1:109484/1‑76 (MQ=255)
                                                            ccATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTtct          <  2:975276/76‑1 (MQ=255)
                                                             cATCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTtctc         >  1:738465/1‑76 (MQ=255)
                                                              aTCTTTGTGGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCg        <  1:921559/76‑1 (MQ=255)
                                                                  ttGTGGACCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGtat    <  2:341421/76‑1 (MQ=255)
                                                                   tgtgGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGtata   >  2:1088845/1‑76 (MQ=255)
                                                                    gtgGAACCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAg  <  1:697785/76‑1 (MQ=255)
                                                                         |                                                                      
TCGGTAATGCCGATAAAGCCGAGGAACGTGGACAGATACGGCGTCACCAGGTCCGTTGGTCCATCTTTGTGGATCCCGCCGCGGCTGGTAATAACGATGGCTTTTTTACCCGTTACCAGACCTTCCGGACCGTTCTCGGTATAG  >  REL606/1451689‑1451832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: