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% 135.8 / NA 45intergenic (‑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 (16/29);  total (16/29)
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:4610234/1‑71 (MQ=255)
 aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:2850216/71‑1 (MQ=255)
 aaCACACGCTGTGTGACTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTc                                                                <  1:3847562/71‑1 (MQ=255)
   cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTcta                                                              <  1:1568548/71‑1 (MQ=255)
     cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            <  1:609812/71‑1 (MQ=255)
     cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            <  1:4024814/71‑1 (MQ=255)
      acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            >  1:3088785/1‑70 (MQ=255)
       cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                           <  1:3101482/70‑1 (MQ=255)
       cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctact                                                            >  1:1584391/1‑69 (MQ=255)
        gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAg                                                         <  1:597160/71‑1 (MQ=255)
           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGcc                                                      >  1:2597845/1‑71 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   <  1:901386/71‑1 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   <  1:2079675/71‑1 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                                    <  1:3268291/70‑1 (MQ=255)
              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAAGCGTTATCACTGGGTTCCAATTCTACTAAGGCCttt                                                   <  1:102395/71‑1 (MQ=255)
                    ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTc                                              >  1:3870080/1‑70 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 >  1:2557324/1‑47 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaaa                                         <  1:3686024/71‑1 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                          <  1:2496606/70‑1 (MQ=255)
                        cGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaa                                           <  1:2672608/69‑1 (MQ=255)
                             aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                          >  1:1558619/1‑65 (MQ=255)
                              cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                          <  1:2451581/64‑1 (MQ=255)
                              cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                   <  1:2896762/71‑1 (MQ=255)
                                 ttCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAAtt                                                                 >  1:405985/1‑38 (MQ=255)
                                  tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                >  1:4170811/1‑70 (MQ=255)
                                    ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtttt                                                  >  1:3795003/1‑50 (MQ=255)
                                    ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                            <  1:3902751/56‑1 (MQ=255)
                                    ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                              <  1:1396602/70‑1 (MQ=255)
                                       ttCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTACGtt                            <  1:2349769/69‑1 (MQ=255)
                                         cGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                          <  1:285294/69‑1 (MQ=255)
                                         cGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                          <  1:1101895/69‑1 (MQ=255)
                                               aaTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                  >  1:4304460/1‑71 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                 <  1:2614121/71‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                 <  1:4351658/71‑1 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                 >  1:1211493/1‑71 (MQ=255)
                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                  >  1:2129908/1‑70 (MQ=255)
                                                 tCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCa                                             <  1:3434684/42‑1 (MQ=255)
                                                      aTCACTGGGTTCCAATTCTACTAAGTCCTTTTCGTCAAAAACCTCAACTCCGTACTCATCAAATTc                >  1:1071746/1‑66 (MQ=255)
                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCa            <  1:980765/70‑1 (MQ=255)
                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt           >  1:560445/1‑71 (MQ=255)
                                                       tCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc          >  1:840402/1‑71 (MQ=255)
                                                         aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt        <  1:4570375/71‑1 (MQ=255)
                                                         aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt         <  1:665972/70‑1 (MQ=255)
                                                             ggTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                     <  1:1259961/38‑1 (MQ=255)
                                                               ttCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa  <  1:2229057/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: