Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,430,136 C→T 23.8% pseudogene (965/1527 nt)
R102H (CGC→CAC) 
lomR →
insH ←
ECK1366:JW5884+JW5904:b4570; Rac prophage region; hypothetical protein
IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,430,1360CT23.8% 78.8 / 22.5 42pseudogene (965/1527 nt)
R102H (CGC→CAC) 
lomR
insH
ECK1366:JW5884+JW5904:b4570; Rac prophage region; hypothetical protein
IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (10/22);  new base T (8/2);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTA  >  W3110S.gb/1430069‑1430197
                                                                   |                                                             
ggTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGAt                                                            >  1:3702095/1‑71 (MQ=255)
    ttGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCgg                                                        >  1:1699079/1‑71 (MQ=255)
      gAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGc                                                       <  1:1456991/70‑1 (MQ=255)
      gAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGc                                                       <  1:3732947/70‑1 (MQ=255)
           aTTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGgcg                                                 <  1:710188/71‑1 (MQ=255)
            ttGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGa                                                             <  1:1852234/58‑1 (MQ=255)
             tGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGa                                                             <  1:1869542/57‑1 (MQ=255)
                 gggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATcc                                           <  1:3847468/71‑1 (MQ=255)
                 gggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATcc                                           <  1:826829/71‑1 (MQ=255)
                        ttGATGCTGCTCGAGGAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATcc                                           <  1:417210/64‑1 (MQ=255)
                        ttGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCa                                          >  1:4112795/1‑65 (MQ=255)
                        ttGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCa                                          >  1:1569449/1‑65 (MQ=255)
                         tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCg                                                         >  1:5106398/1‑49 (MQ=255)
                         tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCTGGGCGCTATCCAGGGa                                      >  1:296196/1‑68 (MQ=255)
                         tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATaa                                   >  1:363701/1‑71 (MQ=255)
                          gATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCt                                               >  1:1395635/1‑58 (MQ=255)
                          gATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAAt                                  >  1:1197852/1‑71 (MQ=255)
                           aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATc                                 <  1:126602/71‑1 (MQ=255)
                                gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGa                                      >  1:1956090/1‑61 (MQ=255)
                                gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGAt                                     >  1:3750387/1‑62 (MQ=255)
                                 ctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa                            <  1:3809067/70‑1 (MQ=255)
                                    gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa                          <  1:1187399/69‑1 (MQ=255)
                                    gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa                          <  1:749821/69‑1 (MQ=255)
                                    gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa                          <  1:2499032/69‑1 (MQ=255)
                                    gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa                          <  1:4130611/69‑1 (MQ=255)
                                           tGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCa                  <  1:147786/70‑1 (MQ=255)
                                            ggCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg                >  1:2922513/1‑71 (MQ=255)
                                             gCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg                >  1:3984114/1‑70 (MQ=255)
                                             gCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCa                  <  1:2287258/68‑1 (MQ=255)
                                              cGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATgg               <  1:3571987/70‑1 (MQ=255)
                                                 aaaTTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg                <  1:1532003/66‑1 (MQ=255)
                                                   aTTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGa         >  1:4114107/1‑71 (MQ=255)
                                                      cATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg                      >  1:2240305/1‑55 (MQ=255)
                                                      cATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGAttt      >  1:2445668/1‑71 (MQ=255)
                                                        tgatgGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg    <  1:4818995/71‑1 (MQ=255)
                                                        tgatgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg    <  1:2483508/71‑1 (MQ=255)
                                                           tggtggTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGAATGGAGGCGATTTCGTa  >  1:677199/1‑70 (MQ=255)
                                                                gtgtGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg    <  1:2650496/63‑1 (MQ=255)
                                                                gtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGc           <  1:4069969/56‑1 (MQ=255)
                                                                gtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGc           <  1:4383714/56‑1 (MQ=255)
                                                                 tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc                   <  1:1265850/47‑1 (MQ=255)
                                                                 tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc                   <  1:2787648/47‑1 (MQ=255)
                                                                  gCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGa         >  1:3004385/1‑56 (MQ=255)
                                                                   |                                                             
GGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTA  >  W3110S.gb/1430069‑1430197

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: