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

GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  AP009048/3824298‑3824426
                                                              |                                                                  
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgtaa                                                             <  1:1414674/70‑3 (MQ=255)
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:2351791/69‑2 (MQ=255)
gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              <  1:760980/69‑2 (MQ=255)
  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)
                                gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                                >  1:2752718/1‑35 (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:1520678/1‑69 (MQ=255)
                                                gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa              >  1:2370119/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)
                                                              |                                                                  
GTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  AP009048/3824298‑3824426

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: