Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,069,243 T→C 8.0% Q21R (CAA→CGA)  PP_4465 ← porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,069,2430TC 8.0% 131.5 / 2.1 50Q21R (CAA→CGA) PP_4465porin
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base C (2/2);  total (23/27)
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.00e+00

TAGCCCGACTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCACGAACGAACATCATCCATCCA  >  NC_002947/5069096‑5069386
                                                                                                                                                   |                                                                                                                                               
tAGCCCGACTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCttgtt                                                                                                                                                 <  6:42069/148‑1 (MQ=255)
  gCCCGACTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGCACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTgc                                                                                                                                               <  7:44272/148‑1 (MQ=255)
      gACTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGcc                                                                                                                                            >  4:316250/1‑147 (MQ=255)
       aCTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAg                                                                                                                                         <  3:316250/149‑1 (MQ=255)
       aCTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAg                                                                                                                                         <  7:295658/149‑1 (MQ=255)
                                 aatggcTGCCTCTGAGCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGGGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTATCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTa                                                                                                                <  1:120320/143‑1 (MQ=255)
                                    gTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGATAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTCGCGCCAAGGTGc                                                                                                                                     <  1:98405/124‑1 (MQ=255)
                                    gTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTCGCGCCAAGGTGc                                                                                                                                     <  1:102407/124‑1 (MQ=255)
                                    gTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTCGCGCCAAGGTGc                                                                                                                                     >  2:102407/1‑124 (MQ=255)
                                    gTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTCGCGCCAAGGTGc                                                                                                                                     >  2:98405/1‑124 (MQ=255)
                                     tCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCtt                                                                                                           <  4:168663/149‑1 (MQ=255)
                                      cTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCttt                                                                                                          >  6:281395/1‑149 (MQ=255)
                                            cTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACtt                                                                                                    >  7:138674/1‑149 (MQ=255)
                                                            gTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCttttt                                                                                                        <  4:48933/129‑1 (MQ=255)
                                                            gTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCttttt                                                                                                        >  3:48933/1‑129 (MQ=255)
                                                              gCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTc                                                                                  >  4:265138/1‑149 (MQ=255)
                                                                 ccGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGc                                                                                       >  3:139127/1‑141 (MQ=255)
                                                                 ccGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGc                                                                                       <  4:139127/141‑1 (MQ=255)
                                                                  cGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCa                                                                                <  3:267735/147‑1 (MQ=255)
                                                                           ccATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGt                                                                                                <  8:98062/122‑1 (MQ=255)
                                                                           ccATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGt                                                                                                >  7:98062/1‑122 (MQ=255)
                                                                                 tagtaCAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCa                                                                                                                            <  3:320814/88‑1 (MQ=255)
                                                                                 tagtaCAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCa                                                                                                                            >  4:320814/1‑88 (MQ=255)
                                                                                       aaGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTAttgt                                                         <  5:281395/149‑1 (MQ=255)
                                                                                         ggTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTAttgttg                                                       <  4:321277/149‑1 (MQ=255)
                                                                                              gagaACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTg                                                   >  4:13847/1‑148 (MQ=255)
                                                                                              gagaACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGt                                                  <  3:265138/149‑1 (MQ=255)
                                                                                               agaACTGATGGATGCCTTACTGTCATCAAGGAAACCCACAGCGTAGCCTTGTTGCGCCAAGGTTCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTc                                                 <  1:22798/149‑1 (MQ=255)
                                                                                                         gatgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTAttgttg                                                       <  4:195458/133‑1 (MQ=255)
                                                                                                         gatgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTAttgttg                                                       >  3:195458/1‑133 (MQ=255)
                                                                                                         gatgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTttatttat                                        <  3:13847/148‑1 (MQ=255)
                                                                                                         gatgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTtata                                       >  2:299138/1‑149 (MQ=255)
                                                                                                           tgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCAt                                                                               <  4:318955/107‑1 (MQ=255)
                                                                                                           tgCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCAt                                                                               >  3:318955/1‑107 (MQ=255)
                                                                                                              cTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCtttattta                                         <  8:221364/142‑1 (MQ=255)
                                                                                                              cTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCtttattta                                         >  7:221364/1‑142 (MQ=255)
                                                                                                               ttACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGAtt                                  >  8:304046/1‑148 (MQ=255)
                                                                                                                   tGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACt                              >  2:278003/1‑148 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTc                                                 <  2:44402/124‑1 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTc                                                 >  1:44402/1‑124 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATg                         >  5:212480/1‑148 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATg                         >  2:295122/1‑148 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGc                        <  1:295122/149‑1 (MQ=255)
                                                                                                                        tcaAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATGGTCTTTATTTATATGATTCACTTAATg                         >  5:207383/1‑148 (MQ=255)
                                                                                                                          aaGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCac                      <  1:278003/149‑1 (MQ=255)
                                                                                                                           aGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCac                      >  4:163939/1‑148 (MQ=255)
                                                                                                                                 cccACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCAcgaacgaa               <  3:163939/149‑1 (MQ=255)
                                                                                                                                  ccACCGCGTAGCCTTGTTGCGCCAATGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCACgaacgaac              <  6:207383/149‑1 (MQ=255)
                                                                                                                                  ccACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCACgaacgaac              <  6:212480/149‑1 (MQ=255)
                                                                                                                                    aCCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCta                                                                          >  4:44004/1‑87 (MQ=255)
                                                                                                                                    aCCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCta                                                                          <  3:44004/87‑1 (MQ=255)
                                                                                                                                               ttgttgCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCACGAACGAACATCAtccatcca  >  6:95413/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TAGCCCGACTGGAAGTCGAGCCTTAAACCGGTCAAGGTCTGCCGCTGATCCGCCCCACCCGTGCTCCGTTGATCCCCATCGTAGTACAAGGTGCGAGAACTGATGGATGCCTTACTGTCATCAAGGAAACCCACCGCGTAGCCTTGTTGCGCCAAGGTGCCGCATGCCAAAGCAGCCAGTATTCTTTTTACTTTTGTCACTACTGCCACTCCATTCCTATAAATATTACTTATTGTTGATTGTCTTTATTTATATGATTCACTTAATGCACGAACGAACATCATCCATCCA  >  NC_002947/5069096‑5069386

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