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% 202.6 / NA 67intergenic (‑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 (33/34);  total (33/34)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATT  >  W3110S.gb/2866044‑2866174
                                                                  |                                                                
aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                              >  1:3793725/1‑71 (MQ=255)
   acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                           >  1:4525412/1‑71 (MQ=255)
   acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                           >  1:1940110/1‑71 (MQ=255)
   acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGACAAATCGTTATCACTGGGTTCCAATTctac                                                           >  1:119839/1‑71 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                          <  1:4245069/71‑1 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                          <  1:586190/71‑1 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                          <  1:3723001/71‑1 (MQ=255)
       gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                       <  1:1680936/71‑1 (MQ=255)
       gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                       <  1:3759402/71‑1 (MQ=255)
        cTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                       <  1:1449433/70‑1 (MQ=255)
         tgtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGTCAAATCGTTATCACTGGGTTCCAATTctacta                                                         <  1:3423621/67‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                 >  1:1657167/1‑71 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                 >  1:2858496/1‑71 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                 >  1:2819021/1‑71 (MQ=255)
                 tCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                             >  1:4434721/1‑71 (MQ=255)
                   ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                          <  1:629314/56‑1 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                          >  1:3570846/1‑70 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                          >  1:3475342/1‑70 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTAGTCaa                                          >  1:2791186/1‑70 (MQ=255)
                      gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         >  1:3664535/1‑70 (MQ=255)
                      gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         >  1:3170659/1‑70 (MQ=255)
                      gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         >  1:702422/1‑70 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                           <  1:2419187/51‑1 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         <  1:3611194/69‑1 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                         >  1:4235016/1‑49 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                            >  1:105073/1‑46 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                            >  1:402250/1‑46 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                       >  1:3185955/1‑51 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         <  1:2857475/64‑1 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         <  1:2810484/64‑1 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   <  1:1064952/70‑1 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   <  1:695120/70‑1 (MQ=255)
                             cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                            >  1:394342/1‑44 (MQ=255)
                              aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                         >  1:4169525/1‑46 (MQ=255)
                              aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                >  1:2233058/1‑71 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                               >  1:4296107/1‑36 (MQ=255)
                                   ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                            <  1:1084278/70‑1 (MQ=255)
                                   ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGt                           >  1:4795687/1‑71 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                          <  1:2808535/39‑1 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                              <  1:2285093/35‑1 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                              <  1:3283497/35‑1 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                       <  1:1624404/42‑1 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtttt                                                <  1:2992929/47‑1 (MQ=255)
                                         ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                      <  1:2328095/70‑1 (MQ=255)
                                         ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                      <  1:4462195/70‑1 (MQ=255)
                                           ccAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                   <  1:760512/71‑1 (MQ=255)
                                            cAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                   >  1:2104003/1‑70 (MQ=255)
                                            cAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                    >  1:1866515/1‑69 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                    <  1:421173/52‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                    <  1:67795/52‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                 <  1:4246752/71‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                 <  1:1052425/71‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                  >  1:702288/1‑70 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                         <  1:3736796/45‑1 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                            >  1:4035109/1‑58 (MQ=255)
                                               aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt               <  1:2608152/71‑1 (MQ=255)
                                                 cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                   <  1:801721/49‑1 (MQ=255)
                                                 cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCGGTTCtcatca                   <  1:3538769/65‑1 (MQ=255)
                                                 cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc             >  1:4035394/1‑71 (MQ=255)
                                                 cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc              >  1:3303645/1‑70 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                 >  1:3047614/1‑50 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg            <  1:39833/71‑1 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg            <  1:2213476/71‑1 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg            >  1:1117934/1‑71 (MQ=255)
                                                  gTTATCACTGGGTTCAAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                 >  1:1280618/1‑50 (MQ=255)
                                                   ttATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                   <  1:643776/63‑1 (MQ=255)
                                                            ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt  >  1:3924079/1‑71 (MQ=255)
                                                                  |                                                                
AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATT  >  W3110S.gb/2866044‑2866174

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: