Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,378 2 bp→AA 100% pseudogene (513‑514/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,3780GA100.0% 212.9 / NA 71pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/37);  total (34/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*AP0090483,824,3790GA100.0% 214.3 / NA 71pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/37);  total (34/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGAC  >  AP009048/3824313‑3824440
                                                                 ||                                                             
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                           <  1:664473/71‑1 (MQ=255)
cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGAGAAGCAATTTg                                                           >  1:2355962/1‑71 (MQ=255)
  gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                          >  1:3037079/1‑70 (MQ=255)
   aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                        >  1:180746/1‑71 (MQ=255)
   aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                        >  1:3451413/1‑71 (MQ=255)
    aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                       >  1:1900461/1‑71 (MQ=255)
    aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                                       <  1:4163764/71‑1 (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)
              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)
                 gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                           >  1:272002/1‑70 (MQ=255)
                 gAAAGGGATGGTAGACGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                           >  1:2116264/1‑70 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                          <  1:968946/70‑1 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                         <  1:3418086/71‑1 (MQ=255)
                  aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                         <  1:4534448/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:5133620/59‑1 (MQ=255)
                                           tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                            <  1:4175433/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)
                                                     aaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                       >  1:3706357/2‑54 (MQ=255)
                                                      aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                       >  1:4059965/1‑53 (MQ=255)
                                                      aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat     >  1:3176559/1‑71 (MQ=255)
                                                      aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat     >  1:170804/1‑71 (MQ=255)
                                                       aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg       >  1:4758057/1‑68 (MQ=255)
                                                        cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga   >  1:3872962/1‑71 (MQ=255)
                                                         ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTa          >  1:472274/4‑63 (MQ=255)
                                                         ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAc  >  1:2319739/4‑71 (MQ=255)
                                                           tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa              >  1:912953/2‑57 (MQ=255)
                                                           tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt        >  1:1121125/2‑63 (MQ=255)
                                                            gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg       <  1:2866328/63‑1 (MQ=255)
                                                                 ||                                                             
CTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGAC  >  AP009048/3824313‑3824440

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: