Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,177,038 A→G 10.4% L316S (TTG→TCG)  PP_4553 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,177,0380AG10.4% 96.3 / 3.4 39L316S (TTG→TCG) PP_4553ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (2/2);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCATGTGCCAGGTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTGGCGAACCAGAGCACAGCCCT  >  NC_002947/5176897‑5177182
                                                                                                                                             |                                                                                                                                                
gCATGTGCCAGGTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGt                                                                                                                                            >  8:336825/1‑148 (MQ=255)
       ccAGGTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAg                                                                                                                                     <  7:177122/148‑1 (MQ=255)
           gTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCGACACGTCATACAGCCcga                                                                                                                                >  6:135005/1‑149 (MQ=255)
           gTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCGACACGTCATACAGCCcga                                                                                                                                >  6:140628/1‑149 (MQ=255)
           gTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCcga                                                                                                                                <  1:73128/149‑1 (MQ=255)
               tAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCGACACGTCATACAGCCCGACGAt                                                                                                                            <  5:140628/149‑1 (MQ=255)
               tAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCGACACGTCATACAGCCCGACGAt                                                                                                                            <  5:135005/149‑1 (MQ=255)
                aGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCACCTCCCCCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCACGGCAACACGGCATACAGCCCGACGAtt                                                                                                                           >  4:223898/1‑149 (MQ=255)
                    cAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTtctc                                                                                                                       <  6:168098/149‑1 (MQ=255)
                    cAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTtct                                                                                                                        <  6:364314/148‑1 (MQ=255)
                          cTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGt                                                                                                                 <  2:231554/149‑1 (MQ=255)
                               cccGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCa                                                                                                                   >  7:333811/1‑142 (MQ=255)
                               cccGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCa                                                                                                                   <  8:333811/142‑1 (MQ=255)
                                       aGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGc                                                                                                    >  6:181668/1‑149 (MQ=255)
                                              cTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTcgcc                                                                                             >  1:51755/1‑149 (MQ=255)
                                               ttCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTcgccg                                                                                            >  3:48799/1‑149 (MQ=255)
                                                     ttCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCgg                                                                                      >  1:223891/1‑149 (MQ=255)
                                                     ttCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCgg                                                                                      >  4:26621/1‑149 (MQ=255)
                                                                       cGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATcg                                                                    >  6:301190/1‑149 (MQ=255)
                                                                            cATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCgccgcc                                                               >  7:109126/1‑149 (MQ=255)
                                                                               cccGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTgg                                                            >  5:197807/1‑149 (MQ=255)
                                                                                   cTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACgg                                                        <  7:45357/149‑1 (MQ=255)
                                                                                      cGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATACACGGCAGCAATCGCCGCCTGGACGGta                                                      <  5:92531/148‑1 (MQ=255)
                                                                                             gcCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGcccc                                              >  8:331032/1‑149 (MQ=255)
                                                                                             gcCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATACGcccc                                              <  7:331032/149‑1 (MQ=255)
                                                                                                aCCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGaa                                           >  2:237853/1‑149 (MQ=255)
                                                                                                   gcgcgGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAAt                                          >  2:58632/1‑147 (MQ=255)
                                                                                                   gcgcgGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAAt                                          <  1:58632/147‑1 (MQ=255)
                                                                                                   gcgcgGTTAAGCGCCACCACCGCGGATGGCCAATCCCGCTGCCACTCGTCATACAGCCCGACGAGTTCTCCCCCGTCCGTCGCTTCAGCCGGCGCCGCTTCGGCAGGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGc                                        <  3:223898/149‑1 (MQ=255)
                                                                                                       ggTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACgg                                    <  7:336825/149‑1 (MQ=255)
                                                                                                                 caccccCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAggg                          <  1:237853/149‑1 (MQ=255)
                                                                                                                     accGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTg                      <  5:301190/149‑1 (MQ=255)
                                                                                                                      ccGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGc                              >  7:31847/1‑140 (MQ=255)
                                                                                                                      ccGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGc                              <  8:31847/140‑1 (MQ=255)
                                                                                                                      ccGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTgg                     <  4:48799/149‑1 (MQ=255)
                                                                                                                      ccGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTgg                     <  2:223891/149‑1 (MQ=255)
                                                                                                                               gCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTGGCGAACCaga            <  5:181668/149‑1 (MQ=255)
                                                                                                                                    tACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTGGCGAACCAGAGcaca       <  8:109126/149‑1 (MQ=255)
                                                                                                                                         cTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTGGCGAACCAGAGCACAGCCCt  <  5:285739/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCATGTGCCAGGTGATAGTCCAGCAACTCCCCCCGCGCCAGAATCCCTTCGACTTCCCGCAAGCCAGCCTCCGGCCCATCCCGCTTCGCCAGCGCCACCGCGCGGTTAAGCGCCACCACCGCTGATGGCCAATACCGCTGCAACACGTCATACAGCCCGACGATTTCTCCCCAGTCCGTCTCTTCAGCCGTCGCCGCTTCGGCATGCACGGCAGCAATCGCCGCCTGGACGGTATAGGCCCCGAATGCACGGCTCTGCAGGGCCTGGCGAACCAGAGCACAGCCCT  >  NC_002947/5176897‑5177182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑