Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 53,278 T→C 7.4% G295G (GGT→GGC opdT‑I → tyrosine‑specific outer membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294753,2780TC 7.4% 153.9 / 2.0 55G295G (GGT→GGCopdT‑Ityrosine‑specific outer membrane porin D
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (2/2);  total (28/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 = 9.88e-01

GAATCAGTACTACGCCGGCACAAGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGCGTTCCTGGA  >  NC_002947/53131‑53426
                                                                                                                                                   |                                                                                                                                                      
gAATCAGTACTACGCCGGCACAAGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGgt                                                                                                                                                        <  5:356070/148‑1 (MQ=255)
         ctacGCCGGCACAAGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTAc                                                                                                                                               <  2:164154/148‑1 (MQ=255)
             gCCGGCACAAGCGTGGTTTACACGATGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGc                                                                                                                                          >  2:87785/1‑149 (MQ=255)
                ggCACAAGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCca                                                                                                                                       >  4:43516/1‑149 (MQ=255)
                      aGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCt                                                                                                                                 >  4:277573/1‑149 (MQ=255)
                        cGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTgg                                                                                                                               <  1:43917/149‑1 (MQ=255)
                         gTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGc                                                                                                                              <  4:292660/149‑1 (MQ=255)
                                      tGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAg                                                                                                                 >  6:173005/1‑149 (MQ=255)
                                              gaTCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGcaacaa                                                                                                         >  7:208237/1‑149 (MQ=255)
                                                    gcgcTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACg                                                                                                       <  8:40401/145‑1 (MQ=255)
                                                    gcgcTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGc‑‑ga                                                                                                 >  3:280923/1‑149 (MQ=255)
                                                       cTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGc‑‑gacga                                                                                              <  3:43516/149‑1 (MQ=255)
                                                              gTGGCTTCAACTACTACAACGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTc                                                                                       <  1:112662/149‑1 (MQ=255)
                                                              gTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGACGCGCCTGGGTGAATTCGACAACAACATCTGGNGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTc                                                                                       <  1:116999/149‑1 (MQ=255)
                                                              gTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTc                                                                                       <  8:72699/149‑1 (MQ=255)
                                                              gTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTc                                                                                       <  5:128017/149‑1 (MQ=255)
                                                              gTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGCGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTc                                                                                       >  2:166129/1‑149 (MQ=255)
                                                                 gCTTCTACTAGTACAAAGCCTTAGTTGAAGGTAAGAAGCGACTGGGTGAATTCGACTACAACATCTGTAGTGCGAAGTTCGGGGTGCTTTACGTTTCCCACACCCTGGCCCTGTCGCATCAGCTCAACAACGGC‑‑GACGACTATTTCGag                                                                                    <  8:35379/149‑2 (MQ=255)
                                                                     cAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTAc                                                                                 <  1:87785/148‑1 (MQ=255)
                                                                            tacAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGt                                                                           >  4:103339/1‑147 (MQ=255)
                                                                            tacAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGAAGATTTCGACTACCTGCGTCa                                                                         >  5:221121/1‑149 (MQ=255)
                                                                            tacAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGCGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGt                                                                           >  7:94933/1‑147 (MQ=255)
                                                                              cAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGt                                                                       >  6:118267/1‑149 (MQ=255)
                                                                              cAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGCTGCCCACACCCTGGCCCTGTCGCATCATCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGt                                                                       >  5:10445/1‑149 (MQ=255)
                                                                               aaaGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTc                                                                      >  7:211024/1‑149 (MQ=255)
                                                                                  gCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGa                                                                   <  4:280923/149‑1 (MQ=255)
                                                                                    cGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACt                                                                 <  8:208237/149‑1 (MQ=255)
                                                                                     gTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTc                                                                >  8:91128/1‑149 (MQ=255)
                                                                                          tGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGAttt                                                           <  7:91128/149‑1 (MQ=255)
                                                                                          tGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGAttt                                                           <  5:173005/149‑1 (MQ=255)
                                                                                             aGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACg                                                                                                                                              <  6:350858/65‑1 (MQ=255)
                                                                                             aGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACg                                                                                                                                              >  5:350858/1‑65 (MQ=255)
                                                                                             aGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGCGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtc                                                         <  8:94933/148‑1 (MQ=255)
                                                                                              ggTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtcgt                                                       >  6:159674/1‑149 (MQ=255)
                                                                                               gTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtcgtc                                                      >  5:47817/1‑149 (MQ=255)
                                                                                                tAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTa                                                                                  >  6:144471/1‑120 (MQ=255)
                                                                                                 aagaagCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGAGTTAAGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtcgtcgc                                                    >  5:350729/1‑148 (MQ=255)
                                                                                                            gggTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCg                                                                     <  4:343422/121‑1 (MQ=255)
                                                                                                            gggTGAATTCGACAACAACATATGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGCCCCTGTCGCATAAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTAGTCGACAACtccatc                                          >  5:467/1‑148 (MQ=255)
                                                                                                             ggTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACtccatc                                          <  3:277573/147‑1 (MQ=255)
                                                                                                                 aaTTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTc                                                                                                                        <  6:123840/67‑1 (MQ=255)
                                                                                                                  aTTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGc‑‑gacgacg                                                                                            <  1:129544/92‑1 (MQ=255)
                                                                                                                  aTTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGc‑‑gacgacg                                                                                            >  2:129544/1‑92 (MQ=255)
                                                                                                                    tCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtc                                                         >  6:344231/1‑125 (MQ=255)
                                                                                                                    tCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTtc                                                         <  5:344231/125‑1 (MQ=255)
                                                                                                                    tCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGt                                      >  4:85066/1‑144 (MQ=255)
                                                                                                                    tCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGt                                      <  3:85066/144‑1 (MQ=255)
                                                                                                                       acaacaacaTCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGAc                              >  7:20554/1‑149 (MQ=255)
                                                                                                                       acaacaacaTCTGGAGTGCCAAGGTCGGTGTGCGTTACGGGGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGAc                              <  6:50269/149‑1 (MQ=255)
                                                                                                                         aacaacaTCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACtt                            >  2:94724/1‑149 (MQ=255)
                                                                                                                                 cTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGCCAGGCGACGATTTCGACTACATGCGTCAGTCGGCCTCGATTTTCGTCGCCAAGTCAAGCAAGTACAGCGAGTTCACCTc                      >  1:115983/1‑149 (MQ=255)
                                                                                                                                 cTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGc                    >  6:370055/1‑149 (MQ=255)
                                                                                                                                    gAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGc                    <  6:47817/146‑1 (MQ=255)
                                                                                                                                         ccAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAaga             <  1:94724/148‑1 (MQ=255)
                                                                                                                                          cAAGGTCGGCGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGc           <  1:166129/149‑1 (MQ=255)
                                                                                                                                             ggTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGAGt         >  4:49691/1‑148 (MQ=255)
                                                                                                                                                   tgtgCGTTACGGTGCCCTCACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCTACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGCGTTCCTGGa  <  6:10445/149‑1 (MQ=255)
                                                                                                                                                   tgtgCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGCGTTCCTGGa  <  1:61732/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GAATCAGTACTACGCCGGCACAAGCGTGGTTTACCCGCTGAGCGACGATCTGGCGCTGATCGGTGGCTTCAACTACTACAAAGCCGTAGATGAAGGTAAGAAGCGCCTGGGTGAATTCGACAACAACATCTGGAGTGCCAAGGTCGGTGTGCGTTACGGTGCCCACACCCTGGCCCTGTCGCATCAGCGCAACAACGGC‑‑GACGACGATTTCGACTACCTGCGTCAGTCGGACTCGATTTTCGTCGACAACTCCATCCAGTACAGCGACTTCAACTCGCCGAAAGAGCGTTCCTGGA  >  NC_002947/53131‑53426

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: