Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,326 T→C 100% pseudogene (459/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,3260TC100.0% 231.2 / NA 67pseudogene (459/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/34);  total (33/34)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824265‑3824392
                                                             |                                                                  
ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaa                                                            <  1:1909233/70‑1 (MQ=255)
 cTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                              <  1:4110753/67‑1 (MQ=255)
   gCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                               <  1:2950411/64‑1 (MQ=255)
    cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:3349747/69‑1 (MQ=255)
    cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:3004055/69‑1 (MQ=255)
    cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                         <  1:2271140/69‑1 (MQ=255)
     ggggCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt                                                               <  1:4482454/60‑1 (MQ=255)
     gTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGAt                                                      >  1:2745132/1‑71 (MQ=255)
        gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                                    <  1:2485265/70‑1 (MQ=255)
              gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                                                                >  1:3521905/3‑52 (MQ=255)
              gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa                                                       <  1:438745/59‑1 (MQ=255)
              gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                              >  1:3834748/3‑70 (MQ=255)
              gcaGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg                                              >  1:3671070/3‑70 (MQ=255)
                  aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                             <  1:4069352/51‑1 (MQ=255)
                  aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                                   <  1:978507/61‑1 (MQ=255)
                  aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:2204145/1‑71 (MQ=255)
                  aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:2829136/1‑71 (MQ=255)
                  aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                         >  1:2537910/1‑71 (MQ=255)
                   tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                         <  1:3121187/70‑1 (MQ=255)
                    gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa                                                           >  1:3946117/1‑51 (MQ=255)
                    gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                       >  1:2570402/1‑71 (MQ=255)
                    gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTc                                       >  1:1635792/1‑71 (MQ=255)
                     cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                      <  1:383704/71‑1 (MQ=255)
                     cgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                      >  1:416032/1‑71 (MQ=255)
                         aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg                                  >  1:91958/1‑71 (MQ=255)
                          cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                   <  1:880646/69‑1 (MQ=255)
                          cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                 >  1:3055816/1‑71 (MQ=255)
                          cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                 <  1:1854961/71‑1 (MQ=255)
                          cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                 <  1:2685767/71‑1 (MQ=255)
                          cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                 >  1:4527166/1‑71 (MQ=255)
                           aGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg                                 >  1:629633/1‑70 (MQ=255)
                            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATg                                                              >  1:2833382/1‑40 (MQ=255)
                                 tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                               <  1:890226/66‑3 (MQ=255)
                                  ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          >  1:1788243/1‑69 (MQ=255)
                                   gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGc                                             <  1:3112847/50‑1 (MQ=255)
                                   gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:4078934/69‑2 (MQ=255)
                                   gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:776216/69‑2 (MQ=255)
                                    tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                       >  1:921742/1‑71 (MQ=255)
                                    tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                       >  1:4231053/1‑71 (MQ=255)
                                      gAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg                                                        >  1:546303/1‑36 (MQ=255)
                                      gAAAACGGAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                    >  1:1034303/1‑70 (MQ=255)
                                        aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                   <  1:3696557/71‑3 (MQ=255)
                                        aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                      >  1:3706427/1‑68 (MQ=255)
                                         aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta                               >  1:366348/1‑56 (MQ=255)
                                          aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                   >  1:3643771/1‑67 (MQ=255)
                                           cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                >  1:3068957/1‑71 (MQ=255)
                                           cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                <  1:3718068/71‑1 (MQ=255)
                                             gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca              <  1:1335597/71‑2 (MQ=255)
                                             gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca              <  1:1318419/71‑2 (MQ=255)
                                               aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat            <  1:407023/71‑4 (MQ=255)
                                               aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat            <  1:3805220/71‑4 (MQ=255)
                                                 gCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc               <  1:1938037/66‑1 (MQ=255)
                                                   tGAAAACCAACCCGTTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                          <  1:1213812/53‑2 (MQ=38)
                                                   tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTaaa                        >  1:1742302/1‑55 (MQ=255)
                                                    gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc               <  1:1153514/63‑1 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTGTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa      >  1:3548455/1‑71 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                                    >  1:699053/1‑41 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc                                   <  1:687596/42‑1 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgtaa                         <  1:2578056/52‑3 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg       <  1:3990401/70‑1 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg       <  1:571632/70‑1 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa      >  1:321600/1‑71 (MQ=255)
                                                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGGTATTGTAAACGGTGAAGCAATTTGCg       <  1:3861988/70‑1 (MQ=255)
                                                      aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt     >  1:3728049/1‑71 (MQ=255)
                                                       aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc    >  1:4313045/1‑71 (MQ=255)
                                                        aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt   >  1:1343977/1‑71 (MQ=255)
                                                         ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg  >  1:4496380/1‑71 (MQ=255)
                                                          cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt                    >  1:4513134/1‑51 (MQ=255)
                                                             |                                                                  
CCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTG  >  W3110S.gb/3824265‑3824392

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: