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

CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824308‑3824439
                                                                  |                                                                 
cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                               >  1:1081814/1‑71 (MQ=255)
   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                            <  1:2038869/71‑3 (MQ=255)
       cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                        <  1:616804/71‑1 (MQ=255)
        tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                        <  1:961848/70‑1 (MQ=255)
        tGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                                        <  1:1046261/70‑1 (MQ=255)
            aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                                                   <  1:1924906/71‑1 (MQ=255)
                  cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                             <  1:2042826/54‑2 (MQ=255)
                  cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              <  1:831918/69‑1 (MQ=255)
                       tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa                                                                 <  1:1575113/46‑1 (MQ=255)
                       tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                                        >  1:466770/1‑71 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:382458/70‑1 (MQ=255)
                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:1139789/70‑1 (MQ=255)
                           aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                             <  1:395934/46‑2 (MQ=255)
                              gATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                                 >  1:1973572/1‑71 (MQ=255)
                               aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                                <  1:2066570/71‑1 (MQ=255)
                                     gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                           >  1:1493688/1‑70 (MQ=255)
                                     gCCGCAGTTTCTGACGGTATTGTAAACGGTGAAGca                                                             <  1:1599791/36‑2 (MQ=16)
                                      ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGGGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:564530/71‑1 (MQ=255)
                                      ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:1524158/71‑1 (MQ=255)
                                          aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:528276/67‑1 (MQ=255)
                                           gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                    <  1:1811802/71‑1 (MQ=255)
                                            tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       <  1:894375/51‑1 (MQ=255)
                                              tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:394915/1‑42 (MQ=19)
                                                  tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                <  1:1762627/68‑1 (MQ=255)
                                                  tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                <  1:361856/68‑1 (MQ=255)
                                                      taTTGTAAACGGTAAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg         >  1:106682/2‑71 (MQ=255)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:443284/41‑1 (MQ=25)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:2122885/41‑1 (MQ=25)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:183882/41‑1 (MQ=25)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:1632251/41‑1 (MQ=25)
                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGt        <  1:2214843/70‑1 (MQ=255)
                                                          gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                >  1:1675030/4‑60 (MQ=255)
                                                             aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  >  1:140281/1‑71 (MQ=255)
                                                             aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  >  1:868350/1‑71 (MQ=255)
                                                               cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtga   >  1:2172271/1‑68 (MQ=255)
                                                                  |                                                                 
CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824308‑3824439

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: