Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 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
*AP0090482,866,1100TA100.0% 274.3 / NA 90intergenic (‑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 (37/53);  total (37/53)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTA  >  AP009048/2866043‑2866176
                                                                   |                                                                  
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 >  1:3400602/1‑71 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 <  1:419405/71‑1 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 <  1:4314520/71‑1 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 <  1:452675/71‑1 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                                  <  1:3398444/70‑1 (MQ=255)
cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAt                                                                  >  1:4762481/1‑70 (MQ=255)
 aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:41625/71‑1 (MQ=255)
 aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:1674190/71‑1 (MQ=255)
  acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                >  1:749150/1‑70 (MQ=255)
  acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:3631036/70‑1 (MQ=255)
   cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTct                                                               >  1:2758409/1‑70 (MQ=255)
     cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            <  1:2663986/71‑1 (MQ=255)
      acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            >  1:2962563/1‑70 (MQ=255)
      acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:868927/66‑1 (MQ=255)
        gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:1649365/64‑1 (MQ=255)
           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTaa                                                          <  1:579132/67‑1 (MQ=255)
           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGcc                                                      >  1:3278980/1‑71 (MQ=255)
             gtgGCTCCCTGCGATAACCGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                              <  1:3141244/61‑1 (MQ=255)
                  tCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   <  1:2839244/67‑1 (MQ=255)
                   cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:5098803/53‑1 (MQ=255)
                   cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               >  1:2906012/1‑70 (MQ=255)
                   cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               >  1:3926137/1‑70 (MQ=255)
                   cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               >  1:2750825/1‑70 (MQ=255)
                   cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               >  1:1718716/1‑70 (MQ=255)
                    ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTc                                              >  1:1895534/1‑70 (MQ=255)
                     cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            <  1:3862598/71‑1 (MQ=255)
                     cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            <  1:5281297/71‑1 (MQ=255)
                     cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            <  1:5296302/71‑1 (MQ=255)
                      tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            >  1:2010323/1‑70 (MQ=255)
                      tGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            >  1:3618070/1‑70 (MQ=255)
                       gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           >  1:4294294/1‑70 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                           <  1:1074656/53‑1 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                         >  1:919490/1‑71 (MQ=255)
                          aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:1329368/67‑1 (MQ=255)
                           tAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:211298/45‑1 (MQ=255)
                            aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                         >  1:4725904/1‑67 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:2684161/62‑1 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                  >  1:2674555/1‑71 (MQ=255)
                               aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                  >  1:520090/1‑71 (MQ=255)
                                gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAc                                        <  1:873601/64‑1 (MQ=255)
                                gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                                 >  1:2591765/1‑71 (MQ=255)
                                gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                                 >  1:2585486/1‑71 (MQ=255)
                                 ttCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                                 >  1:1690096/1‑70 (MQ=255)
                                  tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:1703777/59‑1 (MQ=255)
                                    ctctTCGGCCAAATCGTTATCACTGGGATCCAATTCTACTAAGGCCTTATCGTCaa                                            <  1:3586212/56‑1 (MQ=255)
                                     tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCa                                             <  1:617549/54‑1 (MQ=255)
                                     tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAc                                        <  1:3584623/59‑1 (MQ=255)
                                     tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAc                                        <  1:663898/59‑1 (MQ=255)
                                     tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                   <  1:4496225/64‑1 (MQ=255)
                                      ctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                           >  1:1356640/1‑71 (MQ=255)
                                       ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                          >  1:1796259/1‑55 (MQ=255)
                                       ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                               >  1:3492511/1‑66 (MQ=255)
                                       ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                          >  1:484427/1‑71 (MQ=255)
                                        tCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   >  1:3770176/1‑45 (MQ=255)
                                          ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                        <  1:4704840/70‑1 (MQ=255)
                                          ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                        <  1:1539846/70‑1 (MQ=255)
                                           gCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   <  1:253836/42‑1 (MQ=255)
                                            ccAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           >  1:1305790/1‑49 (MQ=255)
                                             cAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:2217069/48‑1 (MQ=255)
                                             cAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                     >  1:3220218/1‑70 (MQ=255)
                                              aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                   <  1:1474505/71‑1 (MQ=255)
                                              aaaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                   <  1:3323325/71‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:2757255/45‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                   >  1:1428090/1‑69 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                 >  1:2016552/1‑71 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                 <  1:2813618/71‑1 (MQ=255)
                                                  cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc                         >  1:4608346/1‑61 (MQ=255)
                                                  cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc                         >  1:3653531/1‑61 (MQ=255)
                                                  cGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc                >  1:1423712/1‑70 (MQ=255)
                                                   gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               <  1:3899758/38‑1 (MQ=255)
                                                   gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGt                                               <  1:1279038/38‑1 (MQ=255)
                                                   gTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg              <  1:1224376/71‑1 (MQ=255)
                                                    ttATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGc             <  1:2491044/71‑1 (MQ=255)
                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                     >  1:1994053/1‑45 (MQ=255)
                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCa            >  1:3846283/1‑70 (MQ=255)
                                                        cACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:108125/37‑1 (MQ=255)
                                                         aCTTGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt         <  1:4462938/70‑1 (MQ=255)
                                                         aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                    <  1:3794221/59‑1 (MQ=255)
                                                          cTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTc       <  1:2802871/71‑1 (MQ=255)
                                                          cTGAGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt        >  1:3431613/1‑70 (MQ=255)
                                                           tGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc               <  1:469416/62‑1 (MQ=255)
                                                           tGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCa      >  1:1464979/1‑71 (MQ=255)
                                                             ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCa      <  1:4794878/69‑1 (MQ=255)
                                                             ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt    <  1:5221273/71‑1 (MQ=255)
                                                              gTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGt                             <  1:103939/45‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt    <  1:1177973/69‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa  <  1:4812649/71‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa  <  1:4164236/71‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa  <  1:1835342/71‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCATTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa  <  1:4120516/71‑1 (MQ=255)
                                                                   |                                                                  
CAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTA  >  AP009048/2866043‑2866176

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: