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% 183.7 / NA 62pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/36);  total (26/36)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  AP009048/3824293‑3824426
                                                                   |                                                                  
aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                                                                 >  1:3427834/1‑71 (MQ=255)
    ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              >  1:2577600/1‑69 (MQ=255)
    ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                              >  1:4498110/1‑69 (MQ=255)
     gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:4207383/1‑71 (MQ=255)
         aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                        <  1:4312769/71‑2 (MQ=255)
            aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       >  1:438505/1‑67 (MQ=255)
               gCAAGCTGAAAACCAACCCGATGAAAGGTATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:834686/71‑2 (MQ=255)
               gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCACAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                  <  1:1203430/71‑2 (MQ=255)
                 aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                <  1:140021/71‑4 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                <  1:1717685/68‑4 (MQ=255)
                    cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                             >  1:4091924/1‑71 (MQ=255)
                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                                                                 >  1:96296/1‑48 (MQ=39)
                       aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                            >  1:2805649/1‑53 (MQ=255)
                        aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         >  1:2765275/1‑71 (MQ=255)
                          aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                                                        <  1:3230139/54‑2 (MQ=255)
                            cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                     <  1:2354774/71‑1 (MQ=255)
                             aacccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGa                                    <  1:760545/68‑1 (MQ=255)
                              acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                   <  1:3916880/69‑1 (MQ=255)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       >  1:1889710/1‑47 (MQ=38)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                         <  1:3476162/47‑1 (MQ=39)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                     >  1:2769127/1‑51 (MQ=255)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                         <  1:415317/63‑1 (MQ=255)
                                ccGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGGTATTGTAAACggt                                                        <  1:2182992/48‑2 (MQ=37)
                                   aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                              >  1:2788350/1‑71 (MQ=255)
                                   aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                              >  1:286995/1‑71 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                             >  1:2777044/1‑55 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:1625172/71‑1 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:2209204/71‑1 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             <  1:1810266/71‑1 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             >  1:309222/1‑71 (MQ=255)
                                    tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                              >  1:2818997/1‑70 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                       <  1:823771/43‑3 (MQ=38)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:502093/70‑1 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            >  1:2778388/1‑70 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:32826/70‑1 (MQ=255)
                                      aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                            <  1:4200508/70‑1 (MQ=255)
                                       aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:689048/71‑1 (MQ=255)
                                        aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                      <  1:1104164/42‑1 (MQ=37)
                                        aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                           >  1:2159003/1‑69 (MQ=255)
                                          ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc                       <  1:2703001/71‑1 (MQ=255)
                                           gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                      <  1:4224753/71‑1 (MQ=255)
                                             tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:1463095/65‑1 (MQ=255)
                                              ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                           <  1:2484374/63‑1 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                 <  1:2699879/40‑3 (MQ=255)
                                               gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                     <  1:4472966/52‑1 (MQ=255)
                                                tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             >  1:2286174/1‑59 (MQ=255)
                                                   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)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                   >  1:443394/1‑42 (MQ=19)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                   <  1:1477228/42‑1 (MQ=19)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                             >  1:727257/1‑48 (MQ=25)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          >  1:3504584/1‑51 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                <  1:3258081/61‑1 (MQ=255)
                                                           tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                <  1:585111/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)
                                                                cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                          <  1:1820818/46‑1 (MQ=25)
                                                                cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg  >  1:1787282/1‑70 (MQ=255)
                                                                   |                                                                  
AGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATG  >  AP009048/3824293‑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: