Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,368 G→A 100% pseudogene (501/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,3680GA100.0% 156.8 / NA 56pseudogene (501/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/26);  total (30/26)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGT  >  W3110S.gb/3824301‑3824433
                                                                   |                                                                 
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                                <  1:3377709/71‑1 (MQ=255)
  gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTCTAAACg                                                               <  1:474490/70‑1 (MQ=255)
      aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                            >  1:3492240/1‑67 (MQ=255)
           aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                     <  1:1612877/71‑4 (MQ=255)
            aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                                              >  1:1611622/1‑61 (MQ=255)
                gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                 >  1:35930/1‑70 (MQ=255)
                  aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                              >  1:2383292/1‑71 (MQ=255)
                   aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                 >  1:1870317/1‑51 (MQ=255)
                   aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                             >  1:561224/1‑71 (MQ=255)
                   aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                             <  1:1369396/71‑1 (MQ=255)
                   aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGCAAACGGTGAAGCAATTTGCGATc                                             <  1:3040459/71‑1 (MQ=255)
                     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                        >  1:1145129/1‑58 (MQ=255)
                     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                        >  1:2886787/1‑58 (MQ=255)
                     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                           >  1:982169/1‑71 (MQ=255)
                        acccGATGAAAGGGATGGTAGCCTCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        <  1:2997136/69‑1 (MQ=255)
                         cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        >  1:1442181/2‑70 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                   >  1:3180695/1‑71 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                                   >  1:2976517/1‑71 (MQ=255)
                              tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                            <  1:2149288/45‑3 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                          >  1:1421669/1‑46 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                              >  1:3188225/1‑58 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                              >  1:3111530/1‑58 (MQ=255)
                                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 <  1:2008313/69‑1 (MQ=255)
                                  aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGat                               >  1:1902413/1‑69 (MQ=255)
                                    ggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                             <  1:3107361/70‑1 (MQ=255)
                                     gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                              >  1:2606752/1‑52 (MQ=255)
                                       tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        <  1:1200285/56‑1 (MQ=255)
                                       tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:3053789/65‑1 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                          <  1:1056643/52‑1 (MQ=255)
                                           aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     >  1:3464116/1‑71 (MQ=255)
                                              cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     >  1:803879/1‑68 (MQ=255)
                                                  gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                    <  1:1465297/49‑1 (MQ=25)
                                                  gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                    <  1:1328060/49‑1 (MQ=25)
                                                   tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCag                  <  1:877753/66‑1 (MQ=255)
                                                    ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               >  1:1482834/1‑52 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                          <  1:1654586/40‑1 (MQ=16)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        >  1:3590080/1‑42 (MQ=19)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                                      >  1:3547008/1‑44 (MQ=21)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                                >  1:45078/1‑50 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa           >  1:515962/1‑71 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa           >  1:2270233/1‑71 (MQ=255)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACa         >  1:418140/1‑71 (MQ=255)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACa         >  1:3711375/1‑71 (MQ=255)
                                                        gTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   <  1:2964650/60‑1 (MQ=255)
                                                         tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg       <  1:3533300/71‑1 (MQ=255)
                                                          cGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                        >  1:948157/1‑37 (MQ=12)
                                                              aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAc    >  1:1109122/1‑69 (MQ=255)
                                                              aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAc    >  1:910156/1‑69 (MQ=255)
                                                              aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAc    >  1:662375/1‑69 (MQ=255)
                                                              aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg   <  1:1591543/70‑1 (MQ=255)
                                                              aTTGTAAACGGTGAAGCAATTTGCAATCTGGAATACGTTg                                 <  1:706472/40‑1 (MQ=19)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                  <  1:2181588/38‑1 (MQ=21)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                               <  1:2172751/41‑1 (MQ=25)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     <  1:3385335/51‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                     <  1:1670337/51‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                   <  1:1188343/53‑1 (MQ=255)
                                                               ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt  <  1:458719/70‑1 (MQ=255)
                                                                tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                 <  1:229250/38‑1 (MQ=21)
                                                                   |                                                                 
TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGT  >  W3110S.gb/3824301‑3824433

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: