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

TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACA  >  W3110S.gb/3824301‑3824426
                                                             |                                                                
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                         >  1:1059108/1‑71 (MQ=255)
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                         >  1:973517/1‑71 (MQ=255)
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                         >  1:2577962/1‑71 (MQ=255)
tggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc                                                         >  1:1925007/1‑71 (MQ=255)
          cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                               <  1:1195975/71‑3 (MQ=255)
          cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                               <  1:1228974/71‑3 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                        <  1:2498491/71‑1 (MQ=255)
                 aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                        <  1:2082027/71‑1 (MQ=255)
                     ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                    >  1:1497214/1‑71 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                     >  1:1195597/1‑51 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                   <  1:889236/71‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                   <  1:3434670/71‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                   <  1:2305458/71‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                   <  1:2534936/71‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                   <  1:2539083/71‑1 (MQ=255)
                      cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                    >  1:739639/1‑70 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                              <  1:3234656/37‑1 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                              <  1:3284277/37‑1 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                              <  1:515578/37‑1 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                              <  1:2231456/37‑1 (MQ=255)
                             aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                              <  1:614850/37‑1 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:3248934/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:5152/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:611907/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:807298/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:1356871/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:1704871/1‑64 (MQ=255)
                               gAAAGGGATGGTAGCCGCAGTTTCTGTCGCTATTGTaaa                                                          >  1:3236766/1‑39 (MQ=25)
                                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:2061161/71‑1 (MQ=255)
                                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:3147992/71‑1 (MQ=255)
                                 aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:202902/71‑1 (MQ=255)
                                  aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                          >  1:3534871/1‑68 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                      >  1:3394476/1‑65 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                      >  1:1287817/1‑65 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                      >  1:2798881/1‑65 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                >  1:92894/1‑71 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                >  1:590250/1‑71 (MQ=255)
                                         gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                >  1:585965/1‑71 (MQ=255)
                                            gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg              <  1:1955912/70‑1 (MQ=255)
                                             ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg              <  1:3257094/69‑1 (MQ=255)
                                                cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:2921041/56‑1 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                 >  1:2101907/1‑42 (MQ=19)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                      >  1:3261085/1‑53 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                      >  1:1947003/1‑53 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa    >  1:286195/1‑71 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa    >  1:1341868/1‑71 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa    >  1:3406218/1‑71 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTAAGCCGCAGAGACCGa    >  1:2129953/1‑71 (MQ=255)
                                                     tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGA‑‑CAGCCGCAGAGACCGACa  >  1:130714/1‑71 (MQ=255)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:2874335/49‑1 (MQ=25)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:3511756/49‑1 (MQ=25)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:3500341/49‑1 (MQ=25)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:3018751/49‑1 (MQ=25)
                                                       tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                        <  1:3277415/49‑1 (MQ=25)
                                                         tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc   <  1:3161387/68‑1 (MQ=255)
                                                             |                                                                
TGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACA  >  W3110S.gb/3824301‑3824426

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: