Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,139,185 C→T 10.8% Q442* (CAG→TAG)  PP_2754 → OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,139,1850CT10.8% 92.8 / 2.9 37Q442* (CAG→TAG) PP_2754OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CCACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAAT  >  NC_002947/3139040‑3139327
                                                                                                                                                 |                                                                                                                                              
ccACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTGGGCCTGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGt                                                                                                                                             <  2:73992/149‑1 (MQ=255)
               ccTGTGGCTGAAGTACGTGCTTCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCa                                                                                                                              <  3:16377/149‑1 (MQ=255)
               ccTGTGGCTGAAGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGACGCCGAGTCCAATCTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGatat                                                                                                                                >  1:129035/1‑147 (MQ=255)
                     gCTGAAGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCct                                                                                                                          >  7:91886/1‑147 (MQ=255)
                     gCTGAAGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCct                                                                                                                          <  8:91886/147‑1 (MQ=255)
                         aaGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGa                                                                                                                      <  3:97062/147‑1 (MQ=255)
                                             gccggccAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTcgc                                                                                                  >  7:85588/1‑147 (MQ=255)
                                             gccggccAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCga                                                                                                <  7:65821/149‑1 (MQ=255)
                                                ggccAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCAGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCgaga                                                                                              >  6:155017/1‑148 (MQ=255)
                                                    aaGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATc                                                                                                                               <  4:143025/111‑1 (MQ=255)
                                                    aaGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATc                                                                                                                               >  3:143025/1‑111 (MQ=255)
                                                      ggACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCga                                                                                                <  8:105532/140‑1 (MQ=255)
                                                      ggACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCga                                                                                                >  7:105532/1‑140 (MQ=255)
                                                               gTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCTAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAacac                                                                              >  4:239501/1‑149 (MQ=255)
                                                               gTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCTAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAacac                                                                              >  6:109862/1‑149 (MQ=255)
                                                                   aaCCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGAt                                                                          >  8:213373/1‑149 (MQ=255)
                                                                     ccTGATGCAGGCCACCCACCGCGTCGGTGGCGATCTCGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGAt                                                                          <  8:248467/147‑1 (MQ=255)
                                                                                  acccaccGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCaaa                                                            >  4:41290/1‑148 (MQ=255)
                                                                                  acccaccGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCaaa                                                            >  2:172678/1‑148 (MQ=255)
                                                                                     caccGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCTAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACa                                                        <  3:239501/149‑1 (MQ=255)
                                                                                     caccGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCTAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACa                                                        <  5:109862/149‑1 (MQ=255)
                                                                                        cgagTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGc                                                     <  2:129035/146‑1 (MQ=255)
                                                                                                  gCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGt                                           <  7:213373/149‑1 (MQ=255)
                                                                                                              ccGAGTCCAATGTCGACGAACTCCGCCAGATCGTCCAGTAACCGCTCGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCCAGATACCGCCATGCAACACCGAAGGCGGCCTGCCAAAAACAGCCGAGCAGCGGTCGGCAGGTACt                                >  5:166867/1‑148 (MQ=255)
                                                                                                                  gTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGc                           <  8:85588/149‑1 (MQ=255)
                                                                                                                   tCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGc                           <  7:272219/148‑1 (MQ=255)
                                                                                                                   tCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGc                           >  8:272219/1‑148 (MQ=255)
                                                                                                                   tCCAATCTCGACGAACTGCGCCTGATAGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCTAGGAGCGGTCGGCAGGTTCTGAGGc                           >  5:241846/1‑148 (MQ=255)
                                                                                                                      aaTGTGGAGGAACTGCGCCTGATCGTCCAGTAGCCGCTTGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCTTGCAACACCGATGGCGGCCTGGCAAAAACAGGGGAGGAGCGGTCGGCAGGTTCTGAGGCAAgt                       <  6:166867/149‑1 (MQ=255)
                                                                                                                            gacgaACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGc                 >  4:276348/1‑149 (MQ=255)
                                                                                                                                 aCTGCGGCTGATCGTCCAGTACCCTCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGc            <  6:241846/149‑1 (MQ=255)
                                                                                                                                        cTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAg                             >  7:267327/1‑125 (MQ=255)
                                                                                                                                        cTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAg                             <  8:267327/125‑1 (MQ=255)
                                                                                                                                         tGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAgtg                      <  7:244002/131‑1 (MQ=255)
                                                                                                                                         tGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGAGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAgtg                      >  8:244002/1‑131 (MQ=255)
                                                                                                                                           aTCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAAt  >  8:6617/1‑149 (MQ=255)
                                                                                                                                           aTCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAAt  >  4:137876/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCACTGGGAGCGCAACCTGTGGCTGAAGTACGTGCTGCAGAGTGGGCCGGCCAAGGACCTCTCGTTCAACCTGATGCAGGCCACCCACCGCGTCGGTGGCGATCACGTCGCCGAGTCCAATGTCGACGAACTGCGCCTGATCGTCCAGTACCCGCTGGATATCAGCCTCTGAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAAT  >  NC_002947/3139040‑3139327

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: