Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 153,700 A→G 5.9% L576P (CTT→CCT)  htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb153,7000AG 5.9% 251.2 / 6.3 85L576P (CTT→CCT) htrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base G (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGCT  >  W3110S.gb/153633‑153764
                                                                   |                                                                
aaTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGa                                                                <  1:879964/70‑1 (MQ=255)
aaTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGa                                                                <  1:2856313/70‑1 (MQ=255)
aaTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGa                                                                <  1:3444460/70‑1 (MQ=255)
      ctGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAGGATAAACGc                                                         >  1:4566739/1‑71 (MQ=255)
      ctGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAGGATAAACGc                                                         >  1:139722/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGc                                                         <  1:2847915/68‑1 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGc                                                         <  1:3328586/68‑1 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGc                                                         <  1:1867200/68‑1 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:75992/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:2469796/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:5048092/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:1584684/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:4757365/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:711147/1‑71 (MQ=255)
         gAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAACGGTGAAACTAAGATAAACGCTAt                                                      >  1:202978/1‑71 (MQ=255)
          aaaCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATc                                                     <  1:222221/71‑1 (MQ=255)
          aaaCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATc                                                     <  1:1408326/71‑1 (MQ=255)
            aCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAGGATAAACGCTATCGt                                                   <  1:3887918/71‑1 (MQ=255)
            aCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAGGATAAACGCTATCGt                                                   <  1:3881221/71‑1 (MQ=255)
            aCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAGGATAAACGCTATCGt                                                   <  1:1492145/71‑1 (MQ=255)
                gAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAgt                                               >  1:356260/1‑71 (MQ=255)
                gAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAgt                                               >  1:5147838/1‑71 (MQ=255)
                 aaGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGAt                                                               >  1:957086/1‑54 (MQ=255)
                 aaGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGAt                                                               >  1:4644040/1‑54 (MQ=255)
                          gttCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCt                                     >  1:5303661/1‑71 (MQ=255)
                          gttCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCt                                     >  1:3195847/1‑71 (MQ=255)
                          gttCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCt                                     >  1:246275/1‑71 (MQ=255)
                           ttCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCt                                     <  1:4887037/70‑1 (MQ=255)
                           ttCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCt                                     <  1:1866744/70‑1 (MQ=255)
                              aGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCAt                                 <  1:774049/71‑1 (MQ=255)
                              aGTGCCAAGTAATTTTTCAATCGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCAt                                 <  1:5130619/71‑1 (MQ=255)
                                tGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTc                               <  1:4274537/71‑1 (MQ=255)
                                 gCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTc                               >  1:3587404/1‑70 (MQ=255)
                                 gCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTc                               >  1:1502793/1‑70 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:657878/1‑46 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:5285967/1‑46 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:2988434/1‑46 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:714421/1‑46 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTAt                                                      >  1:5201478/1‑46 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTc                                   >  1:4400424/1‑65 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTc                                   >  1:697352/1‑65 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTc                                   >  1:3757545/1‑65 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTc                                   >  1:690598/1‑65 (MQ=255)
                                  ccAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCa                             >  1:731165/1‑71 (MQ=255)
                                    aaGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCAt                                 <  1:2594276/65‑1 (MQ=255)
                                    aaGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATg                           <  1:676950/71‑1 (MQ=255)
                                    aaGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATg                           <  1:3115544/71‑1 (MQ=255)
                                    aaGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATg                           <  1:851668/71‑1 (MQ=255)
                                    aaGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATg                           <  1:4799389/71‑1 (MQ=255)
                                     aGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGa                          <  1:4232504/71‑1 (MQ=255)
                                     aGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGa                          <  1:1130637/71‑1 (MQ=255)
                                     aGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGa                          <  1:3349848/71‑1 (MQ=255)
                                       tAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:5251369/45‑1 (MQ=255)
                                       tAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:4409139/45‑1 (MQ=255)
                                        aaTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:1773543/44‑1 (MQ=255)
                                        aaTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACg                       >  1:4241193/1‑71 (MQ=255)
                                        aaTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACg                       >  1:1199993/1‑71 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:2342771/1‑45 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:4820960/1‑45 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:5323076/1‑45 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGc                      <  1:4886409/71‑1 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGc                      <  1:2452926/71‑1 (MQ=255)
                                         aTTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGc                      <  1:2363743/71‑1 (MQ=255)
                                            tttCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTg                    <  1:526855/70‑1 (MQ=255)
                                                 aTTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:5577994/35‑1 (MQ=255)
                                                 aTTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:4454770/35‑1 (MQ=255)
                                                 aTTGGAATGGTGAAACTAAGATAAACGCTATCGTc                                                  <  1:3922976/35‑1 (MQ=255)
                                                   tGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:5182771/1‑35 (MQ=255)
                                                   tGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:2842542/1‑35 (MQ=255)
                                                   tGGAATGGTGAAACTAAGATAAACGCTATCGTCAg                                                >  1:1051853/1‑35 (MQ=255)
                                                     gAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGa            <  1:5340263/69‑1 (MQ=255)
                                                        tGGTGAAACTAAGATAAACGATATCGTCAGTGTCGCCGTCTTCATTCCATGAACg                       <  1:2244897/55‑1 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGa            <  1:5030863/64‑1 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGa            <  1:4534324/64‑1 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGa            <  1:3382808/64‑1 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTg       >  1:3630634/1‑69 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTa     <  1:1629395/71‑1 (MQ=255)
                                                          gTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTa     <  1:3241350/71‑1 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGaca          >  1:3087025/1‑64 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGaca          >  1:2594587/1‑64 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGc   >  1:122942/1‑71 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGc   >  1:5148981/1‑71 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGc   >  1:3961791/1‑71 (MQ=255)
                                                            gAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGc   >  1:3198038/1‑71 (MQ=255)
                                                             aaaCTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGCt  >  1:1511179/1‑71 (MQ=255)
                                                                   |                                                                
AATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTGTCGCCGTCTTCATTCCATGAACGCTGGGCACTGACACTGTAGCT  >  W3110S.gb/153633‑153764

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: