Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,366 G→A 100% pseudogene (501/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,3660GA100.0% 216.5 / NA 72pseudogene (501/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/41);  total (31/41)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTA  >  AP009048/3824300‑3824432
                                                                  |                                                                  
ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                                >  1:5218398/1‑71 (MQ=255)
 gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                               >  1:2909608/1‑71 (MQ=255)
 gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                               <  1:4858929/71‑1 (MQ=255)
  aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                  >  1:2634112/1‑67 (MQ=255)
   aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                              <  1:4445981/70‑2 (MQ=255)
   aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                              >  1:1718140/1‑69 (MQ=255)
   aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                                <  1:3969534/68‑1 (MQ=255)
    aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                               >  1:1505201/1‑68 (MQ=255)
    aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                            <  1:5118812/71‑1 (MQ=255)
     aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                  >  1:1355205/1‑64 (MQ=255)
      cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                             <  1:4834124/68‑3 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                         >  1:5367070/1‑70 (MQ=255)
        gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGCTATTGTAAACGGTGAAGc                                                         >  1:1179099/1‑70 (MQ=255)
         cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                                 <  1:1067627/61‑1 (MQ=255)
             cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                   <  1:664473/71‑1 (MQ=255)
             cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGAGAAGCAATTTg                                                   >  1:2355962/1‑71 (MQ=255)
               gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                  >  1:3037079/1‑70 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                >  1:3451413/1‑71 (MQ=255)
                aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                >  1:180746/1‑71 (MQ=255)
                 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                               >  1:4660370/1‑55 (MQ=255)
                 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                               <  1:4163764/71‑1 (MQ=255)
                 aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                               >  1:1900461/1‑71 (MQ=255)
                  aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                              <  1:4235171/71‑1 (MQ=255)
                      aacccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                 <  1:1276449/61‑1 (MQ=255)
                        cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                         >  1:202895/2‑70 (MQ=255)
                         ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAAt                                        <  1:2003655/70‑1 (MQ=255)
                         ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                       <  1:1468140/71‑1 (MQ=255)
                         ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATATGGa                                          >  1:1635690/1‑68 (MQ=255)
                           gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                          >  1:494115/1‑50 (MQ=255)
                           gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                 <  1:2871004/59‑1 (MQ=255)
                           gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                                 <  1:2454485/59‑1 (MQ=255)
                           gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                     <  1:2633645/71‑1 (MQ=255)
                           gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                     <  1:1206914/71‑1 (MQ=255)
                             tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                   <  1:4219743/71‑1 (MQ=255)
                              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                           >  1:265434/1‑46 (MQ=255)
                              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                   >  1:272002/1‑70 (MQ=255)
                              gAAAGGGATGGTAGACGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                   >  1:2116264/1‑70 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                  <  1:968946/70‑1 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                 <  1:4534448/71‑1 (MQ=255)
                               aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                 <  1:3418086/71‑1 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                   >  1:2561458/1‑68 (MQ=255)
                                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                <  1:15194/71‑1 (MQ=255)
                                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg                               <  1:1766644/71‑1 (MQ=255)
                                   ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                             >  1:2867265/1‑71 (MQ=255)
                                    gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                            <  1:1594620/71‑1 (MQ=255)
                                     aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                           <  1:4927624/71‑1 (MQ=255)
                                      tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                    >  1:1411444/1‑61 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGACGCAATTTGCGATCt                                             <  1:4867931/50‑1 (MQ=255)
                                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                <  1:1698921/63‑1 (MQ=255)
                                            ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                           <  1:2303058/48‑1 (MQ=25)
                                              gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                  <  1:4479737/71‑1 (MQ=255)
                                              gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                    >  1:2370119/1‑69 (MQ=255)
                                              gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                    >  1:1520678/1‑69 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                    <  1:2260693/52‑1 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                 <  1:3135069/71‑1 (MQ=255)
                                               cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                 <  1:1388429/71‑1 (MQ=255)
                                                 gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc               <  1:1488457/71‑1 (MQ=255)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                   <  1:1866857/50‑1 (MQ=255)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                      <  1:5304609/63‑1 (MQ=255)
                                                  tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc              >  1:1019996/1‑71 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                    >  1:4177099/1‑63 (MQ=255)
                                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa            >  1:3415845/1‑71 (MQ=255)
                                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                     >  1:2099090/1‑44 (MQ=21)
                                                       gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc                                      <  1:5008720/42‑1 (MQ=19)
                                                       gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt         >  1:4176042/1‑71 (MQ=255)
                                                        tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                    <  1:4175433/59‑1 (MQ=255)
                                                        tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                    <  1:5133620/59‑1 (MQ=255)
                                                        tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCTACa          <  1:3334613/69‑1 (MQ=255)
                                                         cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg        >  1:867262/1‑70 (MQ=255)
                                                          ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa      >  1:954892/4‑71 (MQ=255)
                                                          ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa       <  1:3322657/67‑1 (MQ=255)
                                                            taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                >  1:3462160/2‑43 (MQ=25)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg    <  1:1757383/69‑1 (MQ=255)
                                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt   <  1:2453736/70‑1 (MQ=255)
                                                               tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                              >  1:4081326/1‑42 (MQ=25)
                                                               tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa       >  1:3591983/1‑65 (MQ=255)
                                                               tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa  <  1:4759384/70‑1 (MQ=255)
                                                                  |                                                                  
GGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTA  >  AP009048/3824300‑3824432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: