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

GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAG  >  W3110S.gb/3824447‑3824562
                                                                   |                                                 
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                <  1:3207199/71‑1 (MQ=255)
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                <  1:767542/71‑1 (MQ=255)
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                <  1:1533328/71‑1 (MQ=255)
 aCGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                <  1:2057004/70‑1 (MQ=255)
          tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                      >  1:1590700/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGGAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2438004/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGCGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:1456949/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                     >  1:473050/1‑70 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:480426/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:294331/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2069770/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2186837/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2266652/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:3358698/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2667965/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2742406/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:2889033/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:126985/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:1740766/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:1622673/1‑71 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                    >  1:889800/1‑71 (MQ=255)
                aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                 <  1:2658127/70‑1 (MQ=255)
                aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                 <  1:458283/70‑1 (MQ=255)
                aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                 <  1:2617115/70‑1 (MQ=255)
                aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                 <  1:1301582/70‑1 (MQ=255)
                aaGTGCAGGGGACGCCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                 <  1:1391437/70‑1 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:2013276/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:318099/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:1055746/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:2371808/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:561725/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:634294/71‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:1771748/71‑2 (MQ=255)
                      aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                          <  1:1257018/70‑2 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                           >  1:1405682/1‑53 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:26038/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:361133/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:722025/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:2868907/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:2617039/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:2608057/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:1460614/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:2442274/1‑68 (MQ=255)
                       ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                        >  1:167921/1‑68 (MQ=255)
                                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                <  1:3279084/36‑1 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCCCCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:1786002/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:2140511/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:1145726/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:801605/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:1532600/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:1802966/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:2035978/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:3113384/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:3288917/1‑69 (MQ=255)
                                             ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg     >  1:3130371/1‑69 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATGGAATCCATTGTAGCGACGCAg  <  1:343108/70‑1 (MQ=255)
                                                                   |                                                 
GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAG  >  W3110S.gb/3824447‑3824562

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: