Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 905,152 Δ1 bp 10.9% coding (138/324 nt) ybjQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb905,1520T.10.9% 114.1 / 11.9 46coding (138/324 nt)ybjQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base . (3/2);  total (23/23)
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

GCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAA  >  W3110S.gb/905082‑905206
                                                                      |                                                      
gCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGt                                                        >  1:3431039/1‑71 (MQ=255)
   aaGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTcc                                                     <  1:168129/71‑1 (MQ=255)
         ttttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               <  1:2879723/71‑1 (MQ=255)
         ttttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               <  1:1421223/71‑1 (MQ=255)
         ttttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               <  1:3710712/71‑1 (MQ=255)
          tttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               >  1:1488688/1‑70 (MQ=255)
          tttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               >  1:3799902/1‑70 (MQ=255)
          tttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               <  1:683291/70‑1 (MQ=255)
          tttAGGTGCCAATATTTTCCGTGATTTCTTTCCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCg                                               >  1:3527298/1‑70 (MQ=255)
              ggTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATg                                           <  1:3908055/70‑1 (MQ=255)
              ggTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGa                                          <  1:1238570/71‑1 (MQ=255)
                 gCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTCCGGTGCGTa                                             >  1:2667592/1‑64 (MQ=255)
                 gCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTCCGGTGCGTa                                             >  1:281983/1‑64 (MQ=255)
                 gCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTCCGGTGCGTa                                             >  1:695164/1‑64 (MQ=255)
                      tatTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTcc                                                     <  1:3793189/52‑1 (MQ=255)
                      tatTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTcc                                                     <  1:1341507/52‑1 (MQ=255)
                       atTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGtt                                                       >  1:572034/1‑49 (MQ=255)
                       atTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGAGGACGtt                                                       >  1:3337325/1‑49 (MQ=255)
                        ttttCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTg                                >  1:3034553/1‑71 (MQ=255)
                        ttttCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTg                                >  1:485761/1‑71 (MQ=255)
                        ttttCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTg                                >  1:2857674/1‑71 (MQ=255)
                             cGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTa                            <  1:3459039/70‑1 (MQ=255)
                             cGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTa                            <  1:1693055/70‑1 (MQ=255)
                                 aTTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAg                                     <  1:151622/57‑1 (MQ=255)
                                     cTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACgg                    <  1:1779480/70‑1 (MQ=255)
                                          ccGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGgag                 <  1:2940660/68‑1 (MQ=255)
                                            ggTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCAc                      <  1:3627414/61‑1 (MQ=255)
                                            ggTATCCGAGATATCGTTGGCGGACGTCCGGTGCGTATGAAAAAg                                     <  1:2259880/45‑1 (MQ=38)
                                                tCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAg                                     <  1:3309324/42‑1 (MQ=255)
                                                tCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGgag                 <  1:3640317/62‑1 (MQ=255)
                                                tCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGgag                 <  1:2097124/62‑1 (MQ=255)
                                                tCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCtt         <  1:3867008/70‑1 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCTGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:564233/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGAt               >  1:2048707/1‑62 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGAt               >  1:3644304/1‑62 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:3164587/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:2968876/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:3554048/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:1488524/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCttt        >  1:2382120/1‑69 (MQ=255)
                                                  cgcgATATCGTTGGCGGACCTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGgaga                >  1:2316229/1‑61 (MQ=255)
                                                     gATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAgg    <  1:1621506/70‑1 (MQ=255)
                                                      atatCGTTGGCGGACGTCCGGTGCGTATGAAAAAg                                     <  1:166459/35‑1 (MQ=37)
                                                         tCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGgag                 <  1:162062/53‑1 (MQ=255)
                                                         tCGTAGGCGGACGTTCCGGTGCGTATGAAAAAGaa                                   >  1:2947988/1‑35 (MQ=255)
                                                           gTTGGCGGCCGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGaa  >  1:2179410/1‑66 (MQ=255)
                                                           gTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTaaa                          >  1:1025370/1‑42 (MQ=255)
                                                           gTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGaa  >  1:1060524/1‑66 (MQ=255)
                                                           gTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGaa  >  1:1718409/1‑66 (MQ=255)
                                                                      ttCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAgg    <  1:316204/53‑1 (MQ=255)
                                                                      ttCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAgg    <  1:2892888/53‑1 (MQ=255)
                                                                      ttCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAgg    <  1:1582091/53‑1 (MQ=255)
                                                                      |                                                      
GCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAA  >  W3110S.gb/905082‑905206

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: