Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,408 T→A 100% pseudogene (543/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,4080TA100.0% 192.8 / NA 65pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (33/32);  total (33/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAC  >  AP009048/3824344‑3824473
                                                                |                                                                 
ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                             <  1:3503062/71‑1 (MQ=255)
     gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                        >  1:4060063/1‑71 (MQ=255)
      tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                       <  1:116376/71‑1 (MQ=255)
       ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                       >  1:2135402/1‑70 (MQ=255)
        tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                               <  1:585111/61‑1 (MQ=255)
        tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                               <  1:3258081/61‑1 (MQ=255)
         cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                     <  1:4021500/70‑1 (MQ=255)
         cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                     <  1:2209524/70‑1 (MQ=255)
           gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                             >  1:78123/1‑60 (MQ=255)
            tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                 <  1:2989136/71‑1 (MQ=255)
             cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                 >  1:1787282/1‑70 (MQ=255)
                taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:1139677/2‑71 (MQ=255)
                 aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:3243315/1‑70 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  <  1:3321201/64‑1 (MQ=255)
                  ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:175693/1‑69 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGCCCGACATGAACGTAGtgat                                      >  1:498858/1‑71 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                <  1:2565892/61‑1 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      <  1:4065390/71‑1 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      <  1:3992509/71‑1 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      <  1:231623/71‑1 (MQ=255)
                       aaCGGTGAAGCAATTTGCGATCTGGAATACGTCGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      <  1:3842845/71‑1 (MQ=255)
                        aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                             >  1:4531103/1‑47 (MQ=37)
                        aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                                                    >  1:684489/1‑56 (MQ=255)
                        aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg                                     <  1:1045025/71‑1 (MQ=255)
                        aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      >  1:2686611/1‑70 (MQ=255)
                        aCGGTGAAGCAATTTGAGATCTGGAATACGTTGAAGACTCAGCCGc                                                              >  1:4247706/1‑46 (MQ=255)
                         cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTTAGCCGCAGACACCGACATGAACGTAGtgatga                                    >  1:233203/1‑71 (MQ=255)
                          ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                   <  1:3639576/68‑1 (MQ=255)
                          ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc                                   >  1:3672291/4‑71 (MQ=255)
                            tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:1286563/58‑1 (MQ=255)
                                 caatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                                  <  1:1441365/62‑1 (MQ=255)
                                  aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      >  1:615873/3‑60 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                                               >  1:2831525/1‑47 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                       >  1:4123882/1‑71 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                       >  1:3435450/1‑71 (MQ=255)
                                      tGCGATCTGGAATACGTTGAAGACTCAGCCCCAGAGACCGACATGAACGTAGtgatga                                    >  1:1159782/1‑58 (MQ=255)
                                       gCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                       <  1:4344509/70‑1 (MQ=255)
                                         gATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                               <  1:2801426/60‑1 (MQ=255)
                                          aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                    >  1:65114/1‑70 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:4245895/44‑1 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                      >  1:4235304/1‑51 (MQ=255)
                                           tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                 >  1:2099041/1‑56 (MQ=255)
                                              ggAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                 <  1:2624170/37‑1 (MQ=255)
                                               gAATACTTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa              <  1:1136861/71‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACTTGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                >  1:3782438/1‑69 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                  <  1:185978/67‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                    <  1:1889430/65‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg               <  1:4449611/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg               <  1:3142885/70‑1 (MQ=255)
                                               gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg               <  1:3024441/70‑1 (MQ=255)
                                                aaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg               >  1:3044923/1‑69 (MQ=255)
                                                 aTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg               >  1:4198928/1‑68 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGCCCGACATGAACGTAGtga                                       >  1:3249133/1‑41 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg          >  1:1127815/1‑70 (MQ=255)
                                                    cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg          <  1:1067396/70‑1 (MQ=255)
                                                    agTTGAAGACTCAGCCGCAGAGACCGCCATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg          >  1:4127940/2‑70 (MQ=255)
                                                      ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAg       <  1:1914866/71‑1 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                               <  1:3548499/45‑1 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                        >  1:2444480/1‑52 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCAtt                >  1:1981943/1‑60 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg     >  1:3490087/1‑71 (MQ=255)
                                                        gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg     >  1:1046392/1‑71 (MQ=255)
                                                           gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 <  1:1500314/56‑1 (MQ=255)
                                                           gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc  >  1:4119092/1‑71 (MQ=255)
                                                            aCTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg          <  1:631507/62‑1 (MQ=255)
                                                                |                                                                 
CCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAC  >  AP009048/3824344‑3824473

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: