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

AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTC  >  W3110S.gb/2866044‑2866171
                                                                  |                                                             
aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                           <  1:3416006/71‑1 (MQ=255)
 acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                             <  1:439922/68‑1 (MQ=255)
   acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                            >  1:1422547/1‑67 (MQ=255)
   acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATAACTGGGTTCCAAtt                                                            >  1:1210498/1‑67 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                       <  1:1044/71‑1 (MQ=255)
    cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                       <  1:1160240/71‑1 (MQ=255)
     acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                       <  1:3705582/70‑1 (MQ=255)
      cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                       <  1:3665297/69‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                              <  1:1051256/71‑1 (MQ=255)
             tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                               <  1:1092705/70‑1 (MQ=255)
              ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                              <  1:2156554/70‑1 (MQ=255)
               gCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGcc                                                 >  1:1026529/1‑66 (MQ=255)
               gCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtttt                                             <  1:3539177/70‑1 (MQ=255)
                   ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                          <  1:4045248/69‑1 (MQ=255)
                   ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTc                                         >  1:54945/1‑70 (MQ=255)
                    cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                           <  1:490152/51‑1 (MQ=255)
                     tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                            <  1:4235829/49‑1 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                      <  1:4267760/69‑1 (MQ=255)
                       cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                      <  1:1990115/69‑1 (MQ=255)
                        gATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAc                                   <  1:3814456/71‑1 (MQ=255)
                         aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                  <  1:1175823/71‑1 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                    >  1:1852655/1‑67 (MQ=255)
                           aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                      <  1:2263770/65‑1 (MQ=255)
                            aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                <  1:559813/70‑1 (MQ=255)
                            aCAGTTACTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                         >  1:874763/1‑45 (MQ=255)
                               gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                    >  1:1585889/1‑63 (MQ=255)
                               gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                            >  1:3392841/1‑71 (MQ=255)
                               gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                            >  1:3518295/1‑71 (MQ=255)
                                 tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                             >  1:2200665/1‑36 (MQ=255)
                                 tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                          >  1:779462/1‑71 (MQ=255)
                                 tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                          >  1:2018988/1‑71 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTaa                                                     >  1:333346/1‑43 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                               >  1:104960/1‑49 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCt                                                >  1:2022077/1‑48 (MQ=255)
                                  ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                    >  1:4143862/1‑60 (MQ=255)
                                    tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGc                                                  <  1:4136930/44‑1 (MQ=255)
                                     ctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                    >  1:1073468/1‑57 (MQ=255)
                                     ctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                      >  1:119676/1‑71 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTc                                            <  1:2964301/48‑1 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                          <  1:3766533/66‑1 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                     >  1:2293133/1‑71 (MQ=255)
                                      ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                     >  1:4104626/1‑71 (MQ=255)
                                         ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                   <  1:2825819/70‑1 (MQ=255)
                                             aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                  <  1:3612806/51‑1 (MQ=255)
                                                  gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt            >  1:3591470/1‑68 (MQ=255)
                                                      tCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc     >  1:3201305/1‑71 (MQ=255)
                                                       cACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt    >  1:1526565/1‑71 (MQ=255)
                                                        aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                 >  1:2424854/1‑41 (MQ=255)
                                                        aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                 >  1:3739416/1‑41 (MQ=255)
                                                         cTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt   >  1:2673139/1‑70 (MQ=255)
                                                         cTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTc  <  1:573771/71‑1 (MQ=255)
                                                         cTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTc  <  1:579129/71‑1 (MQ=255)
                                                              ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                   <  1:239982/49‑1 (MQ=255)
                                                               tCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                <  1:1465761/51‑1 (MQ=255)
                                                               tCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                <  1:1381046/51‑1 (MQ=255)
                                                                  |                                                             
AACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTC  >  W3110S.gb/2866044‑2866171

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: