Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,290 G→T 8.9% L20M (CTG→ATG) ‡ PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,2900GT 8.9% 137.8 / 4.2 45L20M (CTG→ATG) ‡PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base T (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

TGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAG  >  NC_002947/3057144‑3057438
                                                                                                                                                  |                                                                                                                                                      
tgACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAggg                                                                                                                                                      >  5:69431/1‑149 (MQ=255)
 gACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAg                                                                                                                                                        <  6:86868/146‑1 (MQ=255)
 gACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAg                                                                                                                                                        >  5:86868/1‑146 (MQ=255)
 gACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGgc                                                                                                                                                     >  2:185691/1‑149 (MQ=255)
         caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGgcgc                                                                                                                                                   >  4:238584/1‑143 (MQ=255)
         caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGgcgc                                                                                                                                                   <  3:238584/143‑1 (MQ=255)
         caAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGAca                                                                                                                                             <  7:127325/149‑1 (MQ=255)
            gCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGTCGGGCGCCCTGGCGCTGGCcagca                                                                                                                                          <  6:28377/149‑1 (MQ=255)
                     gcgcTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGt                                                                                                                                 <  1:185691/149‑1 (MQ=255)
                           gCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTaacg                                                                                                                             >  4:181176/1‑145 (MQ=255)
                              cACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAggg                                                                                                                       >  5:151932/1‑149 (MQ=255)
                                ccGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGcc                                                                                                                     >  3:207114/1‑149 (MQ=255)
                                 cGGCAGGGTTTCGGTCACTTCCAGGGTATGCAGGTCGATAAGACTGAGGCTGATGTCCTTATCGATGGAGACGCAGGCGGTCGAAGCGATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCACGATGTGCTCAAGCAGGGccc                                                                                                                    >  8:201413/1‑149 (MQ=255)
                                  ggCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCg                                                                                                                   <  4:133041/149‑1 (MQ=255)
                                   gCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAg                                                                                                                                                        >  1:103720/1‑112 (MQ=255)
                                   gCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAg                                                                                                                                                        <  2:103720/112‑1 (MQ=255)
                                     aGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTg                                                                                                                <  2:38040/149‑1 (MQ=255)
                                            cGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGgc                                                                                                                                                     >  6:232013/1‑106 (MQ=255)
                                            cGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGgc                                                                                                                                                     <  5:232013/106‑1 (MQ=255)
                                             ggTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAggg                                                                                                        >  4:231821/1‑149 (MQ=255)
                                             ggTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAgg                                                                                                         <  4:207114/148‑1 (MQ=255)
                                                 aCTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcg                                                                                                    >  6:232795/1‑149 (MQ=255)
                                                   ttCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                                  >  4:20604/1‑149 (MQ=255)
                                                     cccgCGTTTGCAGGGCGATCAGGCGGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGGGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGggcggc                                                                                                  <  4:255548/144‑1 (MQ=255)
                                                      cAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAag                                                                                                            >  6:235945/1‑136 (MQ=255)
                                                      cAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAag                                                                                                            <  5:235945/136‑1 (MQ=255)
                                                               cAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGAtt                                                                                     >  7:213935/1‑149 (MQ=255)
                                                                 ggTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTcc                                                                                   <  5:232795/149‑1 (MQ=255)
                                                                  gTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAgg                                                                                                                                                       >  7:220410/1‑82 (MQ=255)
                                                                  gTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCt                                                                                  <  2:184409/149‑1 (MQ=255)
                                                                     gATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                                <  3:181176/148‑1 (MQ=255)
                                                                     gATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                                <  6:151932/148‑1 (MQ=255)
                                                                       tCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                                <  8:180240/146‑1 (MQ=255)
                                                                       tCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTa                                                                                >  7:180240/1‑146 (MQ=255)
                                                                       tCAGGCTGAGGCCGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTAttg                                                                             <  3:172212/149‑1 (MQ=255)
                                                                                             ctcGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAgggg                                                        <  3:231821/148‑1 (MQ=255)
                                                                                                gTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGAGTCCTTATTGTTGGGGTTGTCGGGTCAGGGGTAGc                                                    >  6:231028/1‑149 (MQ=255)
                                                                                                         ccAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGTAAAGCAGTGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCAc                                          >  2:99774/1‑149 (MQ=255)
                                                                                                         ccAGGCGCTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGGATGGCGGAGGGCTTATTGTTGTGGTAGGCGGGGCAGGGGTAGCTACACGTCa                                           >  2:87232/1‑149 (MQ=255)
                                                                                                          cAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCAc                                          >  5:35732/1‑149 (MQ=255)
                                                                                                               ggggggaGCCATGGTGTGTCTGCTTGCCAGGGCCAGGGCGGCCGACATCACGGTGGGGGAAAGGCGGGCCCGGTGGAAGATGGGGCGGCGGATGGTCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCg                                      <  1:205401/143‑1 (MQ=255)
                                                                                                                      gCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTc                              <  8:213935/149‑1 (MQ=255)
                                                                                                                           ccccagcctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGATAGTTCCTATGAGGTGTGGGTGTGGGGTGAGGGGTAGCGACACGTCACATCGGGCGTGTGATAGc                         >  1:222141/8‑149 (MQ=12)
                                                                                                                              cagcctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGc                                                                                          >  1:23820/5‑81 (MQ=17)
                                                                                                                              cagcctgctgTCGGGCGCCCTGGCGCTGGCCAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGc                                                                                          <  2:23820/77‑1 (MQ=17)
                                                                                                                                 aCTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGCCGCTTGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAgg                    <  1:87232/149‑1 (MQ=255)
                                                                                                                                         ccAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAtt           <  5:139688/149‑1 (MQ=255)
                                                                                                                                          cAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCAtt           >  2:75592/1‑148 (MQ=255)
                                                                                                                                           aGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTc         >  2:3550/1‑149 (MQ=255)
                                                                                                                                              gcCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGtt      <  3:20604/149‑1 (MQ=255)
                                                                                                                                              gcCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGtt      >  7:57587/1‑149 (MQ=255)
                                                                                                                                                  gggCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAg  <  6:69431/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
TGACAGCAGCAAGCCACGCGGGCGCTGGCCCACCGGCAGGGTTTCGGTGACTTCCAGCGTTTGCAGGTCGATCAGGCTGAGGCTGTTGTCCTTCTCGTTGGAGACCCAGGCGGTGGCAGCCATGGTGTGACTGCTGGCCAGCGCCAGGGCGCCCGACAGCAGGCTGGGGTAAAGCAGGGCCCGGTGGAAGAGGGGGCGGCGCATGGCGGATTCCTTATTGTTGTGGTTGTCGGGTCAGGGGTAGCTACACGTCACTTCGGGCTTGTCATAGCCCAGGCTGTCCATTTCGTTGAAG  >  NC_002947/3057144‑3057438

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: