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,1100TA100.0% 181.5 / NA 61intergenic (‑20/+160)rpoS/nlpDRNA polymerase, sigma S (sigma 38) factor/predicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (23/38);  total (23/38)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCT  >  W3110S.gb/2866044‑2866169
                                                                  |                                                           
aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                         <  1:1093356/71‑1 (MQ=255)
  cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                       >  1:2377693/1‑71 (MQ=255)
  cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                       >  1:2361674/1‑71 (MQ=255)
  cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTct                                                        >  1:1786065/1‑70 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                     <  1:4753443/71‑1 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                     <  1:4857934/71‑1 (MQ=255)
     acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                    <  1:3677922/71‑1 (MQ=255)
     acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                    <  1:3069552/71‑1 (MQ=255)
       gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGGCAAATCGTTATCACTGGGTTCCAATTc                                                         >  1:4882749/1‑64 (MQ=255)
           tgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCt                                              <  1:3928749/71‑1 (MQ=255)
            gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCt                                              >  1:114964/1‑70 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                      >  1:2468860/1‑61 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                      >  1:2724185/1‑61 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                            <  1:3364084/71‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             <  1:5152937/70‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             <  1:296076/70‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCATCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             <  1:4228079/70‑1 (MQ=255)
                cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCg                                         >  1:2203804/1‑71 (MQ=255)
                cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCg                                         >  1:2995499/1‑71 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             <  1:4901025/62‑1 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                     >  1:1690737/1‑70 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                   <  1:3618779/70‑1 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                   <  1:725070/70‑1 (MQ=255)
                          tAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                     <  1:407099/49‑1 (MQ=255)
                          tAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                     <  1:5007232/49‑1 (MQ=255)
                          tAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                >  1:2553854/1‑70 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                       >  1:798699/1‑45 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                        >  1:893687/1‑60 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                        >  1:4928483/1‑60 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                              <  1:1565391/70‑1 (MQ=255)
                              aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTc                                       >  1:4507395/1‑59 (MQ=255)
                               gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                         <  1:1944141/40‑1 (MQ=255)
                                ttCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                      <  1:4671416/42‑1 (MQ=255)
                                ttCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                            <  1:3696822/68‑1 (MQ=255)
                                 tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                       >  1:3979054/1‑40 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                           <  1:3978063/35‑1 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                           <  1:4286152/35‑1 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGc                                                <  1:4564678/46‑1 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             >  1:3470666/1‑49 (MQ=255)
                                   ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                       <  1:4789111/70‑1 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                       <  1:1085035/69‑1 (MQ=255)
                                     ctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                    >  1:1270736/1‑71 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCt                                              >  1:1240654/1‑44 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                     <  1:1258699/53‑1 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                     <  1:761103/69‑1 (MQ=255)
                                       tCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                   >  1:3598037/1‑70 (MQ=255)
                                           ccAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCa                                      <  1:3592898/47‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                             <  1:2534376/38‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                    <  1:4516096/47‑1 (MQ=255)
                                              aaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                   <  1:2732785/63‑1 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt          <  1:4622267/71‑1 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt          <  1:2154072/71‑1 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt          <  1:5172091/71‑1 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt          <  1:645207/68‑1 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg       >  1:3527400/1‑71 (MQ=255)
                                                   ttATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGc      <  1:4361072/71‑1 (MQ=255)
                                                      tCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                    <  1:1352870/38‑1 (MQ=255)
                                                       cACTGGGTTCCAATTCTACTTAGGCCTTTTCGTCAAAAACATCAACTCCGTTCTCATCAAATTCCGCATCt  >  1:221207/1‑71 (MQ=255)
                                                        aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                  >  1:1102528/1‑38 (MQ=255)
                                                           gggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca              <  1:2982835/55‑1 (MQ=255)
                                                           gggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca              <  1:38743/55‑1 (MQ=255)
                                                                  |                                                           
AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCT  >  W3110S.gb/2866044‑2866169

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: