Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,522 T→G 100% pseudogene (657/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,5220TG100.0% 287.5 / NA 81pseudogene (657/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/47);  total (34/47)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  AP009048/3824454‑3824586
                                                                    |                                                                   
atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                                   >  1:1649746/1‑71 (MQ=255)
atcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgct                                                                   >  1:1766467/1‑71 (MQ=255)
 tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                  >  1:3638616/1‑71 (MQ=255)
 tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                  >  1:1686129/1‑71 (MQ=255)
  catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                                 >  1:2127348/1‑71 (MQ=255)
    ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                               >  1:808876/1‑71 (MQ=255)
      gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              <  1:4011340/70‑1 (MQ=255)
      gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                              >  1:1543959/1‑70 (MQ=255)
        aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:3998775/71‑1 (MQ=255)
        aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:1450519/71‑1 (MQ=255)
        aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:3091346/71‑1 (MQ=255)
         gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                        <  1:1215560/71‑1 (MQ=255)
         gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                        <  1:2937648/71‑1 (MQ=255)
          tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                  <  1:4100916/62‑1 (MQ=255)
          tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                       <  1:3200147/71‑1 (MQ=255)
             aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                    <  1:1625377/71‑3 (MQ=255)
             aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                        >  1:401086/1‑67 (MQ=255)
              ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                                                   >  1:3472279/1‑68 (MQ=255)
               gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                      <  1:388939/67‑1 (MQ=255)
                ggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGTGGgaatc                                                 <  1:2204984/71‑1 (MQ=255)
                 gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                <  1:2059803/71‑1 (MQ=255)
                 gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                <  1:4456662/71‑1 (MQ=255)
                    ggCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                          <  1:2738779/59‑1 (MQ=255)
                     gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                            >  1:4274375/1‑71 (MQ=255)
                      cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                            >  1:3612343/1‑70 (MQ=255)
                        gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:2156009/53‑1 (MQ=255)
                         aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                                        >  1:2757799/1‑71 (MQ=255)
                         aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                                        <  1:3388543/71‑1 (MQ=255)
                          aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                       <  1:2909006/71‑1 (MQ=255)
                          aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                    <  1:633878/58‑3 (MQ=255)
                          aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                    <  1:4119877/58‑3 (MQ=255)
                              gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                  >  1:2572589/1‑42 (MQ=38)
                              gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   >  1:1702098/1‑71 (MQ=255)
                               aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                  <  1:1448893/71‑1 (MQ=255)
                                gCCGTTCACGCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                                  >  1:1712121/1‑70 (MQ=255)
                                gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                       >  1:726700/1‑49 (MQ=255)
                                gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                       >  1:1947211/1‑49 (MQ=255)
                                    ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:2498737/41‑1 (MQ=38)
                                    ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:814858/41‑1 (MQ=38)
                                    ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                             <  1:349125/41‑1 (MQ=38)
                                     tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               >  1:2140511/1‑68 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAatccaatcca                                          >  1:1904658/1‑56 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                          >  1:3302379/1‑56 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:856093/70‑1 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:5323230/70‑1 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:5293998/70‑1 (MQ=255)
                                      cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGAGACGCAg                            <  1:2262496/70‑1 (MQ=255)
                                       aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                            <  1:5300224/69‑1 (MQ=255)
                                       aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          <  1:1001932/71‑1 (MQ=255)
                                       aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                           <  1:950119/70‑1 (MQ=255)
                                         ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                        >  1:2778095/1‑71 (MQ=255)
                                         ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                        <  1:4097046/71‑1 (MQ=255)
                                          cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                 >  1:2043063/1‑61 (MQ=255)
                                           aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                          >  1:132913/1‑51 (MQ=255)
                                           aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               <  1:85282/62‑1 (MQ=255)
                                           aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                      >  1:3040977/1‑71 (MQ=255)
                                           aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                      <  1:2774330/71‑1 (MQ=255)
                                           aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACACAGAAggcg                      >  1:3598421/1‑71 (MQ=255)
                                             gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               <  1:2852550/60‑1 (MQ=255)
                                              aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          >  1:1324126/1‑64 (MQ=255)
                                              aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCAACGCAGAAGgcggcg                   >  1:2717052/1‑71 (MQ=255)
                                                gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGAAATCGAATCCATTGTAGc                                   >  1:2390187/1‑53 (MQ=255)
                                                   cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                                     >  1:4981970/1‑48 (MQ=255)
                                                   cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGc                              <  1:2534944/55‑1 (MQ=255)
                                                     aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                          >  1:4265161/1‑57 (MQ=255)
                                                      cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa            <  1:4369559/70‑1 (MQ=255)
                                                         accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                       <  1:1532321/54‑1 (MQ=255)
                                                         accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         >  1:3446913/3‑70 (MQ=255)
                                                         accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt         <  1:4733650/68‑1 (MQ=255)
                                                          ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   <  1:3793334/42‑1 (MQ=37)
                                                          ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   <  1:4647054/42‑1 (MQ=37)
                                                          ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   <  1:4544717/42‑1 (MQ=37)
                                                           cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                   <  1:3439919/42‑1 (MQ=37)
                                                           cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt      >  1:1539683/1‑71 (MQ=255)
                                                             tgttgGCGCTGTCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               <  1:713087/44‑1 (MQ=25)
                                                             tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                               <  1:5249932/44‑1 (MQ=38)
                                                             tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt    >  1:1303322/1‑71 (MQ=255)
                                                             tgttgGCGCTGGACCGAGGGGGAATCGAATCCATTGTAGCGAc                                >  1:3093496/1‑43 (MQ=25)
                                                              gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt     <  1:4712786/69‑1 (MQ=255)
                                                              gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt      >  1:5183624/1‑68 (MQ=255)
                                                              gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGAAGCGACGCAg                            <  1:4577174/46‑1 (MQ=37)
                                                               ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt  <  1:2921643/71‑1 (MQ=255)
                                                                 ggCGCTGCCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt    >  1:4531883/1‑67 (MQ=255)
                                                                    |                                                                   
ATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  AP009048/3824454‑3824586

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: