Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,380 2 bp→AA 100% pseudogene (513‑514/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,3800GA100.0% 186.2 / NA 61pseudogene (513/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/37);  total (24/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*W3110S.gb3,824,3810GA100.0% 186.9 / NA 61pseudogene (514/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/37);  total (24/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824318‑3824439
                                                              ||                                                          
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     <  1:2498491/71‑1 (MQ=255)
aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                                     <  1:2082027/71‑1 (MQ=255)
    ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                 >  1:1497214/1‑71 (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)
     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                                <  1:889236/71‑1 (MQ=255)
     cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                                 >  1:739639/1‑70 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:1704871/1‑64 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:611907/1‑64 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:807298/1‑64 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:5152/1‑64 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:3248934/1‑64 (MQ=255)
              gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                                              >  1:1356871/1‑64 (MQ=255)
                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:202902/71‑1 (MQ=255)
                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:2061161/71‑1 (MQ=255)
                aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:3147992/71‑1 (MQ=255)
                 aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                       >  1:3534871/1‑68 (MQ=255)
                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                   >  1:1287817/1‑65 (MQ=255)
                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                   >  1:2798881/1‑65 (MQ=255)
                        gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGa                                   >  1:3394476/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:1341868/1‑71 (MQ=255)
                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:3406218/1‑71 (MQ=255)
                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                 >  1:286195/1‑71 (MQ=255)
                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTAAGCCGCAGAGACCGa                 >  1:2129953/1‑71 (MQ=255)
                                    tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGA‑‑CAGCCGCAGAGACCGACa               >  1:130714/1‑71 (MQ=255)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:3018751/49‑1 (MQ=25)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:2874335/49‑1 (MQ=25)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:3500341/49‑1 (MQ=25)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:3511756/49‑1 (MQ=25)
                                      tGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                                     <  1:3277415/49‑1 (MQ=25)
                                        tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGAc                <  1:3161387/68‑1 (MQ=255)
                                              ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga                       <  1:1416030/55‑1 (MQ=255)
                                                gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAg      <  1:621635/67‑1 (MQ=255)
                                                gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     >  1:745587/4‑71 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:1142581/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:3644116/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:987745/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:3410120/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:3246209/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:2900071/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:2886979/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:2581344/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:2333084/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:1590206/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:1560313/68‑1 (MQ=255)
                                                 taaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt     <  1:1493848/68‑1 (MQ=255)
                                                   aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat  <  1:2422168/71‑1 (MQ=255)
                                                       gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAgcagca                         <  1:2259053/42‑1 (MQ=25)
                                                       gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:3456046/42‑1 (MQ=25)
                                                       gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa                         <  1:1398661/42‑1 (MQ=25)
                                                              ||                                                          
AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGAT  >  W3110S.gb/3824318‑3824439

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: