Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,410 T→A 100% pseudogene (543/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,4100TA100.0% 195.1 / NA 66pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (37/29);  total (37/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGA  >  W3110S.gb/3824345‑3824474
                                                                 |                                                                
gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              <  1:1637724/70‑1 (MQ=255)
gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              >  1:3338537/1‑70 (MQ=255)
 ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                              <  1:3126838/69‑1 (MQ=255)
 ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:268550/71‑1 (MQ=255)
        ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc                                                      >  1:3429196/1‑70 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:769137/1‑71 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:2055310/1‑71 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:3337848/1‑71 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    >  1:297211/1‑71 (MQ=255)
          cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            >  1:2608084/1‑62 (MQ=255)
             tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGc                                                             <  1:3194150/58‑1 (MQ=255)
                 taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:625125/54‑1 (MQ=255)
                 taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:1160668/70‑1 (MQ=255)
                 taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:2951511/70‑1 (MQ=255)
                 taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:1795261/70‑1 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:2365629/53‑1 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            >  1:3573473/1‑69 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:2707807/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:2701837/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:2421184/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        >  1:232639/4‑71 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:2844892/64‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                       <  1:1720416/69‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                       <  1:2218113/69‑1 (MQ=255)
                      taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                        <  1:3650140/68‑1 (MQ=255)
                        aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                                                          >  1:1244153/1‑50 (MQ=255)
                        aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:1083781/48‑1 (MQ=37)
                        aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg                                                >  1:160377/1‑60 (MQ=255)
                        aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                                     >  1:3351123/1‑71 (MQ=255)
                         aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:463985/63‑1 (MQ=255)
                                   aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                               >  1:914381/3‑50 (MQ=255)
                                   aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                               >  1:2683293/3‑50 (MQ=255)
                                   aatttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                               >  1:843218/3‑50 (MQ=255)
                                     tttGCGATCTGGAATACGTTGAAGACTCAGCCGCa                                                            <  1:796084/35‑1 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg                                >  1:909573/1‑61 (MQ=255)
                                       tGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgcg                      >  1:1707406/1‑71 (MQ=255)
                                           aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                   >  1:1678987/1‑70 (MQ=255)
                                            tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                                                    <  1:2983975/36‑1 (MQ=255)
                                            tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                            <  1:3169389/44‑1 (MQ=255)
                                            tCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatca                 >  1:3288813/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                        >  1:1884286/1‑63 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                        >  1:2090870/1‑63 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACAGACATGa                                               >  1:1231359/1‑40 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                   <  1:242523/65‑1 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTg              <  1:2265588/70‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              <  1:2022905/51‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              <  1:3218558/51‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGt          <  1:3456398/71‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGt          <  1:915097/71‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtg                                       <  1:1207914/40‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              >  1:1607506/1‑49 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                              >  1:210645/1‑49 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:3264117/1‑51 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:3209851/1‑51 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:3518825/1‑51 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                            >  1:1003029/1‑51 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcat                   >  1:2862005/1‑60 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGaa            >  1:1660560/1‑67 (MQ=255)
                                                      gTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatc                  <  1:2135707/60‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                       <  1:2366133/54‑1 (MQ=255)
                                                         gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                        >  1:2134151/1‑51 (MQ=255)
                                                         gAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    >  1:2259508/1‑71 (MQ=255)
                                                          aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    <  1:1682526/70‑1 (MQ=255)
                                                          aaGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAggg    <  1:2483961/70‑1 (MQ=255)
                                                           aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa  >  1:3174300/1‑71 (MQ=255)
                                                           aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGa  >  1:3211073/1‑71 (MQ=255)
                                                                 |                                                                
GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGA  >  W3110S.gb/3824345‑3824474

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: