Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,050,010 C→A 9.1% L15M (CTG→ATG) ‡ oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,050,0100CA 9.1% 123.5 / 3.8 44L15M (CTG→ATG) ‡oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base A (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGC  >  NC_002947/4049882‑4050152
                                                                                                                                |                                                                                                                                              
cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGTGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGccc                                                                                                                              <  5:184474/147‑1 (MQ=255)
cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCCGCTGGCCAGGGCCAgcagcagggtgg                                                                                                                              <  3:71324/147‑7 (MQ=38)
  acgcgGGCAAGCCCGCTCCTACAAGGCCCCTTGAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                          <  7:280605/148‑1 (MQ=255)
          aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac                                                                                                                   >  1:149698/1‑148 (MQ=255)
              ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct                                                                                                              >  7:76395/1‑149 (MQ=255)
              ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCCCCACTGCCCTTGCCAGCGACGAtct                                                                                                              >  2:251748/1‑149 (MQ=255)
              ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct                                                                                                              >  1:284867/1‑149 (MQ=255)
               cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct                                                                                                              >  8:251769/1‑148 (MQ=255)
               cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                                                                                                             <  8:76395/149‑1 (MQ=255)
                   ccTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa                                                                                                         >  4:296783/1‑149 (MQ=255)
                        aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                          <  6:87852/127‑1 (MQ=255)
                        aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                          >  5:87852/1‑127 (MQ=255)
                        aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCAGCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCCCTGCCCTTGCCAGCGCCGATCTCAGGACCCgc                                                                                                    >  4:247806/1‑149 (MQ=255)
                           gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag                                                                                                  >  5:9219/1‑148 (MQ=255)
                                  tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACt                                                                                                                                  >  8:14318/1‑109 (MQ=255)
                                  tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACt                                                                                                                                  <  7:14318/109‑1 (MQ=255)
                                       gCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGt                                                                                     <  7:229156/149‑1 (MQ=255)
                                         gggCTTGCCCTCGAATGCCCCGCAAACTCATCCCGATATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGGTGCTGGCCCTGGCCAGCACCACTACCCTTTCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGTTGa                                                                                   <  8:224540/149‑1 (MQ=255)
                                          ggCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGAc                                                                                  <  3:167950/149‑1 (MQ=255)
                                                cccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                                                                                                                 <  6:73177/112‑1 (MQ=255)
                                                cccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                                                                                                                 >  5:73177/1‑112 (MQ=255)
                                                 ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGCCTGGggcg                                                                           >  2:340442/1‑149 (MQ=255)
                                                 ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACACCCCCCCTGCTGCTGGCCCTCGCAAGAACCACTGCCCTTGCCACCGACGATCTCATGAACCGCAGCACACTGAACGCTGACTGGggcg                                                                           >  8:205044/1‑149 (MQ=255)
                                                   gcgAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcggc                                                                         <  2:149698/149‑1 (MQ=255)
                                                           cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCa                                                                 >  3:137966/1‑149 (MQ=255)
                                                                     cAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaa                                                       <  3:296783/149‑1 (MQ=255)
                                                                      aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTgcg                                                                    >  5:169445/1‑135 (MQ=255)
                                                                      aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTgcg                                                                    <  6:169445/135‑1 (MQ=255)
                                                                           gagaTTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGCCCGGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                 >  5:364011/1‑149 (MQ=255)
                                                                           gagaTTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                 >  5:310019/1‑149 (MQ=255)
                                                                                    tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAg                                        <  3:335284/149‑1 (MQ=255)
                                                                                               gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                             <  7:251769/149‑1 (MQ=255)
                                                                                                   ttCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGCTCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGCCGGCGTCAAGCTCCCCGGCGACTa                          >  2:277434/1‑148 (MQ=255)
                                                                                                     cGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca                                                                                                   >  1:355147/1‑73 (MQ=255)
                                                                                                     cGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca                                                                                                   <  2:355147/73‑1 (MQ=255)
                                                                                                       cATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCg                     <  4:137966/149‑1 (MQ=255)
                                                                                                         tCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                   <  2:284867/149‑1 (MQ=255)
                                                                                                             ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGaaa               >  6:331736/1‑149 (MQ=255)
                                                                                                             ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGaaa               <  1:251748/149‑1 (MQ=255)
                                                                                                                   agtggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGcc          >  6:71381/6‑148 (MQ=255)
                                                                                                                    ctgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATAAACTCGAAGAAGACGGCGTCAACATCACCGGCGAATACAGAGGTGAAACCGcat         >  6:180527/1‑146 (MQ=255)
                                                                                                                     tggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTa        <  1:101471/145‑1 (MQ=255)
                                                                                                                     tggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTa        >  2:101471/4‑148 (MQ=255)
                                                                                                                      gctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa      >  7:162451/1‑149 (MQ=255)
                                                                                                                          ctggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACg   <  6:9219/148‑1 (MQ=255)
                                                                                                                          ctggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACAGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGc  >  4:75130/1‑149 (MQ=255)
                                                                                                                             gccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                             >  7:274842/1‑119 (MQ=255)
                                                                                                                             gccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                             <  8:274842/119‑1 (MQ=255)
                                                                                                                                ctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACaa     <  3:197831/140‑1 (MQ=255)
                                                                                                                                ctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACaa     >  4:197831/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                              
CTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGC  >  NC_002947/4049882‑4050152

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: ‑