Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,533:1 +G 100% pseudogene (666/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,5331.G100.0% 258.6 / NA 75pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (31/44);  total (31/44)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597
                                                                       |                                                                
aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                   <  1:6395146/71‑1 (MQ=255)
 aGTGCAGGGGACGGGAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                 <  1:4604734/71‑1 (MQ=255)
     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:1947380/1‑70 (MQ=255)
        gggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                             <  1:7788002/68‑2 (MQ=255)
           aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                                        <  1:1829609/71‑1 (MQ=255)
                aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                                   <  1:874521/71‑1 (MQ=255)
                aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                                   <  1:2354629/71‑1 (MQ=255)
                aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                 >  1:4164726/1‑56 (MQ=255)
                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:2139122/1‑71 (MQ=255)
                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                       >  1:4936076/1‑63 (MQ=255)
                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:2317236/1‑55 (MQ=255)
                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:2707366/1‑55 (MQ=255)
                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:3953243/1‑55 (MQ=255)
                         gCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:587410/1‑71 (MQ=255)
                          ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                         <  1:5469895/71‑1 (MQ=255)
                          ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:4413152/1‑49 (MQ=255)
                          ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                                              >  1:5145938/1‑49 (MQ=255)
                           cGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                                        <  1:434113/71‑1 (MQ=255)
                            gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:2184134/53‑1 (MQ=255)
                            gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGc                                       <  1:754001/71‑1 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:563260/70‑1 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:939873/1‑70 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:1133733/1‑70 (MQ=255)
                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:118785/70‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:1362021/69‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:5987580/69‑1 (MQ=255)
                                aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:6390437/69‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                                        <  1:2411538/64‑1 (MQ=255)
                                  ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg                                 >  1:4978768/1‑71 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:5028531/66‑1 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:4405816/66‑1 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:3974716/66‑1 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:3828196/66‑1 (MQ=255)
                                   cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                                >  1:1497411/1‑71 (MQ=255)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                                   >  1:2885762/1‑35 (MQ=37)
                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                               >  1:613779/1‑71 (MQ=255)
                                      gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                                                 <  1:5185575/51‑1 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  >  1:925024/1‑49 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            <  1:1942244/71‑1 (MQ=255)
                                       aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            <  1:1932923/71‑1 (MQ=255)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                                       >  1:4814406/1‑43 (MQ=25)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                           <  1:5166287/71‑1 (MQ=255)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                           >  1:535655/1‑71 (MQ=255)
                                        agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc                           <  1:477527/71‑1 (MQ=255)
                                         gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCt                          <  1:7861937/71‑1 (MQ=255)
                                          agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:3347950/39‑1 (MQ=25)
                                          agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                         <  1:1156647/39‑1 (MQ=25)
                                           gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     >  1:2011530/1‑42 (MQ=25)
                                           gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     >  1:6919679/1‑42 (MQ=25)
                                           gCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:3442749/1‑53 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                                      >  1:2789926/1‑56 (MQ=255)
                                            cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:3301187/57‑1 (MQ=255)
                                              aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:3092102/1‑55 (MQ=255)
                                               cTCACCTTGTTGGCGCTGTCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa                     <  1:614989/70‑1 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:1214244/1‑43 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:2628075/1‑43 (MQ=25)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                                   >  1:106136/1‑56 (MQ=255)
                                               cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa                     >  1:3010704/1‑70 (MQ=255)
                                                tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                                >  1:2177203/1‑58 (MQ=255)
                                                  accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:2396534/39‑1 (MQ=25)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:5535833/39‑1 (MQ=37)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:6116878/39‑1 (MQ=37)
                                                   ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:615981/39‑1 (MQ=37)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:5628676/39‑1 (MQ=37)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:6974619/39‑1 (MQ=37)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:4175538/1‑49 (MQ=255)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt               >  1:7760573/1‑71 (MQ=255)
                                                    cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt               <  1:1814349/71‑1 (MQ=255)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:3111647/1‑43 (MQ=37)
                                                     ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt               <  1:4593070/70‑1 (MQ=255)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGATTCCATTGt                                               >  1:3879258/1‑37 (MQ=25)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               >  1:2667718/1‑37 (MQ=37)
                                                      tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt               >  1:6185108/1‑69 (MQ=255)
                                                         tgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                                <  1:1812781/49‑1 (MQ=255)
                                                         tgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                                <  1:1122502/49‑1 (MQ=255)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                                      <  1:2834532/42‑1 (MQ=37)
                                                          ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa          <  1:2535501/70‑1 (MQ=255)
                                                           gcgctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTaa         <  1:6710025/67‑1 (MQ=255)
                                                           gcgctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAg        >  1:1623123/4‑71 (MQ=255)
                                                              ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg                                  >  1:2493049/1‑42 (MQ=38)
                                                              ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg                                  >  1:84577/1‑42 (MQ=38)
                                                              ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg                                  >  1:93850/1‑42 (MQ=38)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:5224802/36‑1 (MQ=37)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:72309/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1615528/71‑1 (MQ=255)
                                                                 gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:1663000/71‑1 (MQ=255)
                                                                       |                                                                
AAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824463‑3824597

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: