Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,866,110 T→A 100% intergenic (‑20/+160) rpoS ← / ← nlpD RNA polymerase, sigma S (sigma 38) factor/predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,866,1100TA91.0% 111.8 / 2.3 43intergenic (‑20/+160)rpoS/nlpDRNA polymerase, sigma S (sigma 38) factor/predicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base A (20/19);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATT  >  W3110S.gb/2866043‑2866174
                                                                   |                                                                
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                               >  1:948072/1‑71 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                                <  1:2676749/70‑1 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                                <  1:1935041/70‑1 (MQ=255)
  acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTct                                                             <  1:3480717/71‑1 (MQ=255)
         cTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAgg                                                      <  1:3019376/71‑1 (MQ=255)
            tgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTaa                                                        >  1:1537329/1‑66 (MQ=255)
            tgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGc                                                     <  1:1784875/69‑1 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                 >  1:2923924/1‑71 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                                  <  1:3432547/70‑1 (MQ=255)
                    ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTc                                            <  1:1421958/70‑1 (MQ=255)
                      tGCGATAACAGTTCCTCTTCGGCCAAATAGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                          >  1:3022819/1‑70 (MQ=255)
                       gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         >  1:1834188/1‑70 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAgg                                                      >  1:3812486/1‑56 (MQ=255)
                            aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                    >  1:3372575/1‑70 (MQ=255)
                            aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                    >  1:3146062/1‑70 (MQ=255)
                             aCAGTTCCTCTTCGGCCAAATCGTTTTCACTGGGTTCCAATTCTACTAAGGCCTTTTc                                               >  1:3253444/1‑58 (MQ=255)
                             aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                       >  1:746689/1‑66 (MQ=255)
                             aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   <  1:3137724/70‑1 (MQ=255)
                             aCAGTTACTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   >  1:21924/1‑70 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTCTTCTACTAAGGCCtttt                                                <  1:2334027/55‑1 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTCTTCTACTAAGGCCtttt                                                <  1:2425689/55‑1 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTCTTCTACTAAGGCCtttt                                                <  1:3186482/55‑1 (MQ=255)
                                gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                               <  1:1156016/71‑1 (MQ=255)
                                  tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                       >  1:3067895/1‑61 (MQ=255)
                                  tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                       >  1:3278172/1‑61 (MQ=255)
                                   ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                             >  1:2263235/1‑70 (MQ=255)
                                    ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                            <  1:1895835/70‑1 (MQ=255)
                                             cAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                  <  1:2534720/71‑1 (MQ=255)
                                               aaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                             <  1:2028913/58‑1 (MQ=255)
                                               aaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                <  1:1629927/71‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCTCTTCTACTAAGGCCtttt                                                <  1:2077683/38‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCa                                  <  1:670129/52‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCa                                  <  1:1559796/52‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt               >  1:2944762/1‑71 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt               <  1:1686035/71‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                >  1:109112/1‑70 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                >  1:3397135/1‑70 (MQ=255)
                                                 tCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                    >  1:3490240/1‑49 (MQ=255)
                                                 tCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   >  1:1677563/1‑50 (MQ=255)
                                                             ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                 >  1:2483527/1‑40 (MQ=255)
                                                             ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt  <  1:1554089/71‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                    <  1:1200216/51‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                    <  1:885081/51‑1 (MQ=255)
                                                                   |                                                                
CAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATT  >  W3110S.gb/2866043‑2866174

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: