Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,906,143 T→A 5.1% N152I (AAC→ATC)  eno ← enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,906,1430TA 5.1% 270.8 / 5.5 97N152I (AAC→ATC) enoenolase
Reads supporting (aligned to +/- strand):  ref base T (38/54);  new base A (2/3);  total (40/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGCTCGTACAG  >  AP009048/2906074‑2906210
                                                                     |                                                                   
ccAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGtt                                                                    <  1:469472/71‑1 (MQ=255)
 cAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGtt                                                                    >  1:5913877/1‑70 (MQ=255)
   aCCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                                                 <  1:7367406/71‑1 (MQ=255)
    ccGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACGACCGTTGATGATGTtcatc                                                                <  1:2186032/71‑1 (MQ=255)
    ccGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatc                                                                <  1:4985401/71‑1 (MQ=255)
    ccGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatc                                                                <  1:5635515/71‑1 (MQ=255)
    ccGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatc                                                                <  1:3653942/71‑1 (MQ=255)
     cGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:3016361/71‑1 (MQ=255)
     cGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:3628695/71‑1 (MQ=255)
      ggCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatc                                                                <  1:2628620/69‑1 (MQ=255)
      ggCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTcatcat                                                              <  1:7567408/71‑1 (MQ=255)
      ggCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTcatcat                                                              <  1:4710398/71‑1 (MQ=255)
      ggCCGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTcatcat                                                              <  1:1536388/71‑1 (MQ=255)
       gCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCatcatc                                                             >  1:7410857/1‑71 (MQ=255)
        cTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtca                                                                  <  1:5272065/65‑1 (MQ=255)
        cTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtca                                                                  <  1:3764865/65‑1 (MQ=255)
        cTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCg                                                            <  1:7164800/71‑1 (MQ=255)
        cTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCg                                                            <  1:6306166/71‑1 (MQ=255)
          gAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCgg                                                           >  1:4751231/1‑70 (MQ=255)
           aaTCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCgg                                                           <  1:4567978/69‑1 (MQ=255)
           aaTCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGa                                                          <  1:6968298/70‑1 (MQ=255)
             tCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAAcc                                                       >  1:5452252/1‑71 (MQ=255)
             tCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAAcc                                                       <  1:4960075/71‑1 (MQ=255)
                  aaTTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCAt                                                  <  1:7699522/71‑1 (MQ=255)
                    ttCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCa                                                   >  1:6295524/1‑68 (MQ=255)
                    ttCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCa                                                   >  1:3846966/1‑68 (MQ=255)
                     tCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 >  1:5443707/1‑69 (MQ=255)
                          gATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGc                                                    >  1:7420493/1‑61 (MQ=255)
                          gATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGc                                                    >  1:5192292/1‑61 (MQ=255)
                          gATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAtt                                          >  1:6970550/1‑71 (MQ=255)
                          gATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                                           <  1:7803346/70‑1 (MQ=255)
                           atatCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:5544116/1‑71 (MQ=255)
                           atatCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                                           <  1:6488848/69‑1 (MQ=255)
                            tatCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTg                                        <  1:6383556/71‑1 (MQ=255)
                            tatCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTg                                        <  1:6282047/71‑1 (MQ=255)
                              tCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGcc                                      >  1:7864531/1‑71 (MQ=255)
                              tCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGcc                                      <  1:657388/71‑1 (MQ=255)
                              tCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGcc                                      >  1:5503348/1‑71 (MQ=255)
                               cAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGccc                                     >  1:2429090/1‑71 (MQ=255)
                               cAACGTTGTTGTCAGCGTGCTCACCACCGATGATGAACATCATCAACgg                                                           <  1:7244931/49‑1 (MQ=18)
                               cAACGTTGTTGTCAGCGTGCTCACCACCGATGATGAACATCATCAACgg                                                           <  1:1050225/49‑1 (MQ=18)
                               cAACGTTGTTGTCAGCGTGCTCACCACCGATGATGAACATCATCAACgg                                                           <  1:2634446/49‑1 (MQ=18)
                                aaCGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                                           >  1:524397/1‑64 (MQ=255)
                                aaCGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCg                                    <  1:7266411/71‑1 (MQ=255)
                                aaCGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCg                                    >  1:2999564/1‑71 (MQ=255)
                                 aCGTTGTTGTCAGCGTGCTAACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCg                                    >  1:2236950/1‑70 (MQ=255)
                                  cGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:4692987/1‑64 (MQ=255)
                                  cGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:3080132/1‑64 (MQ=255)
                                  cGTTGTTGTCAGCGTGCTCACCACCGATGATGAACATCATCAACg                                                            >  1:6926872/1‑45 (MQ=17)
                                    ttgttgTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:7736259/40‑1 (MQ=255)
                                    ttgttgTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:1687341/40‑1 (MQ=255)
                                    ttgttgTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:6033695/40‑1 (MQ=255)
                                    ttgttgTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:2791170/40‑1 (MQ=255)
                                    ttgttgTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                                               <  1:4201473/40‑1 (MQ=255)
                                     tgttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTa                               >  1:4069573/1‑71 (MQ=255)
                                      gttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTAc                              >  1:433932/1‑71 (MQ=255)
                                       ttgtCAGCGTGCTCTCCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               <  1:4133991/53‑1 (MQ=255)
                                       ttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               <  1:1856632/53‑1 (MQ=255)
                                       ttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               <  1:7889062/53‑1 (MQ=255)
                                       ttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               <  1:4844875/53‑1 (MQ=255)
                                       ttgtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTAcc                             <  1:145377/71‑1 (MQ=255)
                                         gtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGt                                             >  1:3856133/1‑53 (MQ=255)
                                         gtCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCg                            <  1:3699178/70‑1 (MQ=255)
                                          tCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgtt                          <  1:2783854/71‑1 (MQ=255)
                                          tCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgtt                          >  1:6730006/1‑71 (MQ=255)
                                           cAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 <  1:3791502/47‑1 (MQ=255)
                                           cAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 <  1:6340608/47‑1 (MQ=255)
                                             gCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               >  1:2039377/1‑47 (MQ=255)
                                             gCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 <  1:7819377/45‑1 (MQ=255)
                                             gCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcag                       <  1:4115960/71‑1 (MQ=255)
                                              cGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCAt                                                  >  1:1081137/1‑43 (MQ=255)
                                              cGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 >  1:6415824/1‑44 (MQ=255)
                                              cGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                                                 >  1:3864018/1‑44 (MQ=255)
                                                tGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagttc                    <  1:5852612/71‑1 (MQ=255)
                                                 gCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagtt                     >  1:6713282/1‑69 (MQ=255)
                                                 gCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagtt                     >  1:6160750/1‑69 (MQ=255)
                                                  cTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                               >  1:868494/1‑42 (MQ=255)
                                                    caccacCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGccc                                     >  1:1551638/1‑50 (MQ=255)
                                                    caccacCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGc                 >  1:138805/1‑70 (MQ=255)
                                                     accaccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGc                 <  1:1440260/69‑1 (MQ=255)
                                                     accaccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCg                <  1:1991545/70‑1 (MQ=255)
                                                       caccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:7812264/1‑43 (MQ=255)
                                                       caccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:7821570/1‑43 (MQ=255)
                                                       caccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAtt                                          >  1:7037690/1‑42 (MQ=255)
                                                         ccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                                           <  1:6941832/39‑1 (MQ=255)
                                                         ccGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                                           <  1:6139813/39‑1 (MQ=255)
                                                          cGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGc          <  1:4084079/71‑1 (MQ=255)
                                                           gTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCtttcagt                      <  1:1026060/58‑1 (MQ=255)
                                                            ttGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGc                 <  1:3266791/62‑1 (MQ=255)
                                                            ccgatgatGAACATCATCGGAACCGGCATAGAGTATTTGCCAGGAGTACCGTTCAGTTCAGCGATGTGCTc        >  1:5006076/3‑71 (MQ=255)
                                                              gatgatGTTCATCATCGGAACCGGCATAGAGTAttt                                         >  1:6104600/1‑36 (MQ=255)
                                                              gatgatGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGa               >  1:2345275/1‑62 (MQ=255)
                                                              gatgatGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGCTCg       >  1:7869817/1‑70 (MQ=255)
                                                                tgatgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTa                               >  1:2306978/1‑44 (MQ=255)
                                                                 gatgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCATCg                <  1:6069501/58‑1 (MQ=255)
                                                                 gatgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCg                <  1:6812205/58‑1 (MQ=255)
                                                                 gatgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCg                <  1:1140565/58‑1 (MQ=255)
                                                                 gatgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCg                <  1:930395/58‑1 (MQ=255)
                                                                  atgTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGCTCGTACAg  <  1:1961123/71‑1 (MQ=255)
                                                                     |                                                                   
CCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGCTCGTACAG  >  AP009048/2906074‑2906210

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: ‑