Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,824,362 A→T 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,3620AT98.2% 164.1 / ‑5.7 54pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (31/22);  minor base G (1/0);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGC  >  W3110S.gb/3824297‑3824412
                                                                 |                                                  
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1010909/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1072125/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1188960/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1254352/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1379343/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:1407065/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:686191/1‑68 (MQ=255)
cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  >  1:2033662/1‑68 (MQ=255)
          aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  <  1:110838/58‑1 (MQ=255)
          aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  <  1:3769150/58‑1 (MQ=255)
          aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                  <  1:1700588/58‑1 (MQ=255)
          aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCCAGTTTCTGTCGGTAt                                                  <  1:1965917/57‑1 (MQ=255)
            ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                >  1:60841/1‑58 (MQ=255)
            ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                >  1:2842378/1‑58 (MQ=255)
            ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                >  1:3136947/1‑58 (MQ=255)
            ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtg                                                >  1:546438/1‑58 (MQ=255)
                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                            <  1:3145416/69‑1 (MQ=255)
                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                           <  1:1529465/70‑1 (MQ=255)
                     aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg                           <  1:3906090/70‑1 (MQ=255)
                        aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCt                       >  1:4001186/1‑71 (MQ=255)
                        aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                            >  1:2981635/1‑50 (MQ=255)
                        aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc                        >  1:569867/1‑70 (MQ=255)
                        aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                       >  1:3045409/1‑71 (MQ=255)
                        aCCAACCCGATGAAAGGGATGGTAGCCGAAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt                       >  1:3262449/1‑71 (MQ=255)
                            acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:3228940/69‑1 (MQ=255)
                            acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:3584757/69‑1 (MQ=255)
                            acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:2176847/69‑1 (MQ=255)
                            acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:1955239/69‑1 (MQ=255)
                            acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:838538/69‑1 (MQ=255)
                             cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:154754/69‑1 (MQ=255)
                                gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg               <  1:2215413/71‑1 (MQ=255)
                                gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg               <  1:3075801/71‑1 (MQ=255)
                                gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg               <  1:764257/71‑1 (MQ=255)
                                   gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGaa                                     >  1:792651/1‑46 (MQ=255)
                                   gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt             >  1:3656858/1‑70 (MQ=255)
                                   gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt             >  1:1067936/1‑70 (MQ=255)
                                   gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt             >  1:2905889/1‑70 (MQ=255)
                                   gAAAGGGATGGTAGCCGCAGTTTCTGTCGGGATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt             >  1:3806210/1‑70 (MQ=255)
                                    aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg            <  1:3860395/70‑1 (MQ=255)
                                    aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg            <  1:2534711/70‑1 (MQ=255)
                                    aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg            <  1:2667950/70‑1 (MQ=255)
                                    aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg            <  1:3519229/70‑1 (MQ=255)
                                     aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg            <  1:1148967/69‑1 (MQ=255)
                                             gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc  >  1:560493/1‑71 (MQ=255)
                                             gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc  >  1:939555/1‑71 (MQ=255)
                                             gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc  >  1:1435717/1‑71 (MQ=255)
                                             gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc  >  1:3849542/1‑71 (MQ=255)
                                             gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc  >  1:345405/1‑71 (MQ=255)
                                               aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   <  1:1172387/52‑1 (MQ=255)
                                                  cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCTATCTGGaa                   >  1:3992993/1‑49 (MQ=255)
                                                  cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   >  1:1501502/1‑49 (MQ=25)
                                                  cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   >  1:3339105/1‑49 (MQ=25)
                                                  cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                   >  1:827531/1‑49 (MQ=25)
                                                  cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTAGCGATCTGGaa                   >  1:2921406/1‑49 (MQ=255)
                                                                 |                                                  
CTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGC  >  W3110S.gb/3824297‑3824412

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: