Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,620,028 G→A 10.0% R69R (CGC→CGT)  opdH ← tricarboxylate‑specific outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,620,0280GA10.0% 106.4 / 2.3 40R69R (CGC→CGT) opdHtricarboxylate‑specific outer membrane porin
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base A (2/2);  total (21/19)
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.62e-01

GGGTGGGCAGCAGGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTGGCGATG  >  NC_002947/1619883‑1620172
                                                                                                                                                  |                                                                                                                                                
gggtgggCAGCAGGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCgcgc                                                                                                                                               >  6:121573/1‑149 (MQ=255)
 ggtgggCAGCAGGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCt                                                                                                                                              >  4:241288/1‑149 (MQ=255)
   tgggCAGCAGGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTg                                                                                                                                            >  8:6690/1‑149 (MQ=255)
      gcagcaGGCCGGTGTCGGTACGGTCGGGGCTTGAGTTGAGCTTGATGCCGAGG‑ATGCCCAAGGCGTCCAGTCCAAAGCCCACTGTTCAATCGGTGTAGGCAGACTCGAAATTGAGTATGAAGCCATGGGCCCATTCATCGCGCTTggtt                                                                                                                                         <  8:252658/149‑3 (MQ=255)
        agcagGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGCGATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCACGCTTGGAtt                                                                                                                                        >  1:202749/1‑149 (MQ=255)
        agcagGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCACGCTTGGATtg                                                                                                                                       >  1:157325/1‑149 (MQ=255)
                      ggTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGc                                                                                                                         >  8:240201/1‑149 (MQ=255)
                       gTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATtg                                                                                                                                       <  7:197293/134‑1 (MQ=255)
                       gTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATtg                                                                                                                                       >  8:197293/1‑134 (MQ=255)
                                         tCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTTCTGCGAGCTGGTGCCGTCGCGGAAGTAGCgggt                                                                                                      >  4:142948/1‑147 (MQ=255)
                                            aGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGGTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTg                                                                                                     >  7:100370/1‑147 (MQ=255)
                                            aGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGaa                                                                                                   <  7:6690/149‑1 (MQ=255)
                                                   tGCCGAGC‑CTGCCCAACGAGTCCAGGCCAAAGCCCACGGTGCGATCGGGGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCGATTCATCGCGCGTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGGGGCGGTTGAAGTACATg                                                                                            >  6:249045/1‑149 (MQ=255)
                                                   tGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATg                                                                                            <  7:240201/149‑1 (MQ=255)
                                                         gC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCgg                                                                                      <  8:100370/149‑1 (MQ=255)
                                                               cccAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGttt                                                                                  <  5:121573/148‑1 (MQ=255)
                                                                 cAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCg                                                                                >  7:155935/1‑148 (MQ=255)
                                                                   aCGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAg                                                                             <  7:8181/149‑1 (MQ=255)
                                                                   aCGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGAATGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTAACGAAg                                                                             >  2:227618/1‑149 (MQ=255)
                                                                       gTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGt                                                                                    >  2:8394/1‑138 (MQ=255)
                                                                       gTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGt                                                                                    <  1:8394/138‑1 (MQ=255)
                                                                          cAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGTGGTTGAAGTACATGGTGCGGGTTTCGAAGGTGGcc                                                                       <  1:34198/148‑1 (MQ=255)
                                                                                      cACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCACGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCTGTGCTGTCTTCga                                                          <  2:202749/149‑1 (MQ=255)
                                                                                      cACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCACGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCTGTGCTGTCTTCga                                                          <  2:157325/149‑1 (MQ=255)
                                                                                              cATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAggcg                                                  <  3:241288/149‑1 (MQ=255)
                                                                                                         ccAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCAt                                       >  8:128199/1‑149 (MQ=255)
                                                                                                               tCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGc                                 >  2:22575/1‑149 (MQ=255)
                                                                                                                 gAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGa                               >  2:163120/1‑149 (MQ=255)
                                                                                                                  aaaTTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGa                               >  8:14710/1‑148 (MQ=255)
                                                                                                                       gAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGa                               >  5:6108/1‑143 (MQ=255)
                                                                                                                       gAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGa                               <  6:6108/143‑1 (MQ=255)
                                                                                                                        aGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCgg                         <  7:128199/148‑1 (MQ=255)
                                                                                                                          gATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGc                        >  5:207079/1‑147 (MQ=255)
                                                                                                                           aTGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAg                     >  7:63543/1‑149 (MQ=255)
                                                                                                                           aTGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAg                     >  3:2319/1‑149 (MQ=255)
                                                                                                                              aaGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAgcgc                  <  1:22575/149‑1 (MQ=255)
                                                                                                                                         ccATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTgg       <  7:14710/149‑1 (MQ=255)
                                                                                                                                         ccATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTgg       >  4:239148/1‑149 (MQ=255)
                                                                                                                                         ccATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTgg       <  1:163120/149‑1 (MQ=255)
                                                                                                                                          cATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTGGc      <  3:93029/149‑1 (MQ=255)
                                                                                                                                               aTCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTGGCGATg  <  8:155935/148‑1 (MQ=255)
                                                                                                                                                  gcgcTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTgg                                            >  2:183923/1‑103 (MQ=255)
                                                                                                                                                  gcgcTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTgg                                            <  1:183923/103‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGGTGGGCAGCAGGCCGCTGTCGGTACGGTCGGGGCTGGAGTCGAGCTTGATGCCGAGC‑ATGCCCAACGCGTCCAGGCCAAAGCCCACGGTGCCATCGGTGTAGCCAGACTCGAAATTGAGGATGAAGCCCTGGGCCCATTCATCGCGCTTGGATTGCTGCGCGCTGGTGCCGTCGCGGAAGTCGCGGTTGAAGTACATGTTGCGGGTTTCGAAGGTGGCCGTGCTGTCTTCGAAGAAGGCGGCCTGGCTCATCGGGGCGACACCGGCAAGCGCGAGGGCGCTGGCGATG  >  NC_002947/1619883‑1620172

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