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

AAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824305‑3824432
                                                               |                                                                
aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           <  1:973809/71‑3 (MQ=255)
 aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           >  1:687221/1‑68 (MQ=255)
  aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           >  1:2065231/1‑67 (MQ=255)
  aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           >  1:2223891/1‑67 (MQ=255)
   cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                        >  1:1117007/1‑71 (MQ=255)
    ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                                       <  1:1114151/71‑1 (MQ=255)
     gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      <  1:459703/71‑2 (MQ=255)
      cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                     <  1:2483399/71‑3 (MQ=255)
      cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                     <  1:1545927/71‑3 (MQ=255)
       aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                                    <  1:618545/71‑4 (MQ=255)
        aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                           <  1:394980/63‑3 (MQ=255)
          cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                                         >  1:3342808/1‑63 (MQ=255)
           tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                >  1:3640684/1‑71 (MQ=255)
            gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                >  1:1784038/1‑70 (MQ=255)
             aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                                               <  1:2658297/70‑1 (MQ=255)
             aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                              >  1:249979/1‑71 (MQ=255)
              aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                             >  1:650640/1‑71 (MQ=255)
               aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                >  1:3282768/1‑67 (MQ=255)
               aaccaaccAGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                          >  1:2950482/1‑57 (MQ=255)
                aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                                           >  1:583756/1‑71 (MQ=255)
                  cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                          >  1:1385818/1‑70 (MQ=255)
                    acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                       <  1:1293884/69‑1 (MQ=255)
                         aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                              <  1:1270416/59‑1 (MQ=255)
                           gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAATCAATTTGCGa                                              >  1:1381513/1‑57 (MQ=255)
                             aaGGGATTGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:2974998/71‑1 (MQ=255)
                              aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                >  1:3555599/1‑68 (MQ=255)
                                 gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                          <  1:1559881/71‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:2940540/63‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                         <  1:1253981/54‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                         <  1:1247329/54‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:3700645/63‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:2170792/63‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:1212511/63‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:1081231/65‑1 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                          >  1:518440/1‑69 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                          >  1:131747/1‑69 (MQ=255)
                                   tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                          >  1:651598/1‑69 (MQ=255)
                                    ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                         <  1:3609772/69‑1 (MQ=255)
                                    ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                         <  1:3636939/69‑1 (MQ=255)
                                     gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                      >  1:630595/1‑71 (MQ=255)
                                        gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                                >  1:2512161/1‑42 (MQ=19)
                                          cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                    >  1:1853139/1‑68 (MQ=255)
                                                ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              >  1:3583759/1‑52 (MQ=255)
                                                ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc            >  1:554003/1‑70 (MQ=255)
                                                ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAcc            >  1:225258/1‑70 (MQ=255)
                                                 tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              >  1:1870366/1‑51 (MQ=255)
                                                 tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                  >  1:3279704/1‑63 (MQ=255)
                                                 tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCg           <  1:3199485/70‑1 (MQ=255)
                                                 tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa          >  1:1063349/1‑71 (MQ=255)
                                                     tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg      <  1:2718677/71‑1 (MQ=255)
                                                       ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGaa    <  1:3193241/68‑1 (MQ=255)
                                                          aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              >  1:1198428/1‑42 (MQ=25)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   <  1:2736220/36‑1 (MQ=16)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   <  1:669687/36‑1 (MQ=16)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   <  1:2845961/36‑1 (MQ=16)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:3810711/41‑1 (MQ=25)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                               <  1:3037496/40‑1 (MQ=25)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg  <  1:2457143/69‑1 (MQ=255)
                                                           ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg  <  1:762581/69‑1 (MQ=255)
                                                            tgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa                               <  1:3401806/39‑1 (MQ=21)
                                                               |                                                                
AAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824305‑3824432

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: