Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,657,604 C→T 8.4% A216T (GCC→ACC)  PP_3222 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,657,6040CT 8.4% 133.0 / 2.5 48A216T (GCC→ACC) PP_3222ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (2/2);  total (25/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 = 8.67e-01

CGTTGCGCGAGAAGATGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGCCGGAAACCAG  >  NC_002947/3657463‑3657752
                                                                                                                                             |                                                                                                                                                    
agTTGCGCGAGAAGATGGTTTCCACGAATACAGCCGCGGACACGGTTTTGCGGACTTGCAGGCCGACGAGCGTAATTAGCGGGAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCGGGCTAAGGCCGTTGGAGAAGGCg                                                                                                                                               <  4:182168/148‑1 (MQ=255)
    gcgcGAGAAGATGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCggtg                                                                                                                                            <  2:76054/148‑1 (MQ=255)
      gcgAGAAGATGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCggtggt                                                                                                                                          >  8:204223/1‑148 (MQ=255)
        gagaAGATGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGTGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCGTTGGCGAAGGCGGTGGTGAt                                                                                                                                       >  8:272971/1‑149 (MQ=255)
              aTGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTAc                                                                                                                                   <  4:332481/147‑1 (MQ=255)
              aTGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTAc                                                                                                                                   >  3:332481/1‑147 (MQ=255)
                  tttCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGGCCGGCTAAGGCCCTTGGCGAAGGCggtg                                                                                                                                            <  1:203715/134‑1 (MQ=255)
                  tttCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCggtg                                                                                                                                            >  2:203715/1‑134 (MQ=255)
                     ccaccaATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTg                                                                                                                           <  7:202557/148‑1 (MQ=255)
                      caccaATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGcc                                                                                                                         >  5:332681/1‑149 (MQ=255)
                       accaATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCca                                                                                                                        >  4:316539/1‑149 (MQ=255)
                            tACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGAt                                                                                                                                       <  2:42433/129‑1 (MQ=255)
                            tACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGAt                                                                                                                                       >  1:42433/1‑129 (MQ=255)
                            tACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGc                                                                                                                   <  4:75562/149‑1 (MQ=255)
                               aGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCt                                                                                                                >  8:123104/1‑149 (MQ=255)
                                 ccccGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTg                                                                                                               >  8:16076/1‑148 (MQ=255)
                                 ccccGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCAAACGCCCTg                                                                                                               >  2:116467/1‑148 (MQ=255)
                                   ccGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGTGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGa                                                                                                              >  4:199951/1‑147 (MQ=255)
                                       acacGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCcg                                                                                                        <  7:120587/149‑1 (MQ=255)
                                        cacGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCcgc                                                                                                       <  7:16076/149‑1 (MQ=255)
                                            gTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTg                                                                                                   <  3:316539/149‑1 (MQ=255)
                                            gTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTg                                                                                                   <  1:182791/149‑1 (MQ=255)
                                                     aCTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCTGCTAAGGCCCTTGGCGAAGGCGGTCGTGATGTACGGCTCCTCCCACACGCCCTGAAAGCCGCGATGCATCACCTg                                                                                          >  6:283514/1‑149 (MQ=255)
                                                     aCTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTg                                                                                          >  4:103896/1‑149 (MQ=255)
                                                       ttGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCt                                                                                                                <  6:36473/125‑1 (MQ=255)
                                                       ttGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCt                                                                                                                >  5:36473/1‑125 (MQ=255)
                                                            ggCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCGCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTa                                                                                    >  1:341026/1‑148 (MQ=255)
                                                                   cAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGTGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAg                                                                            >  5:362616/1‑149 (MQ=255)
                                                                    aGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGTGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAgc                                                                           <  3:199951/149‑1 (MQ=255)
                                                                        tCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTg                                                                       <  7:204223/149‑1 (MQ=255)
                                                                               cGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGTAGCACCTGGGCGTACACCGCAGCGCTGGGGATCg                                                                >  3:170602/1‑149 (MQ=255)
                                                                                                         ggCGCGCCTCAACCTGGGCCCGGCTCAGGCCCTTGGCGAAGGCGGGGGTGATGTACGGCGCCTGCCACACGCCCTGACAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAgg                                      <  8:320804/149‑1 (MQ=255)
                                                                                                            gTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTg                                   <  2:341026/149‑1 (MQ=255)
                                                                                                                cTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCgg                                                >  7:175908/1‑132 (MQ=255)
                                                                                                                cTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCgg                                                <  8:175908/132‑1 (MQ=255)
                                                                                                                cTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGc                                 >  2:73399/1‑147 (MQ=255)
                                                                                                                cTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGc                                 >  8:282004/1‑147 (MQ=255)
                                                                                                                 tGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAg                              >  3:147884/1‑149 (MQ=255)
                                                                                                                    aCCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGcc                           <  3:103896/149‑1 (MQ=255)
                                                                                                                       tGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTg                        >  2:212865/1‑149 (MQ=255)
                                                                                                                        gggCCCGGCTAAGGCCCTTGGTGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTgg                       <  6:362616/149‑1 (MQ=255)
                                                                                                                         ggCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTgg                       <  1:212865/148‑1 (MQ=255)
                                                                                                                               gCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCg                >  3:5659/1‑149 (MQ=255)
                                                                                                                                cTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACTTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGt               <  7:123104/149‑1 (MQ=255)
                                                                                                                                   aGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGc            >  5:309349/1‑149 (MQ=255)
                                                                                                                                    ggCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGcc           <  7:282004/149‑1 (MQ=255)
                                                                                                                                       ccTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCATGGCTGCCCGTGGCCgg         >  7:126801/1‑148 (MQ=255)
                                                                                                                                           ggCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGTAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGTGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGCCGGAAAc     <  4:170602/148‑1 (MQ=255)
                                                                                                                                             cGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGCCGGAAACCAg  <  3:249050/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CGTTGCGCGAGAAGATGGTTTCCACCAATACAGCCCCGGACACGGTATTGCCGACTTGCAGGCCGACCAGCGTCAATAGCGGCAGCGCCGCGTTGCGCAGGCCATGGCGTGCCTGAACCTGGGCCCGGCTAAGGCCCTTGGCGAAGGCGGTGGTGATGTACGGCTCCTGCCACACGCCCTGAAAGCCGCGCTGCAGCACCTGGGCGTACACCGCAGCGCTGGGGATCGCCAGGGTGATGGCCGGCAGCACCAGGCTGGCCAGGCCTTGGCTGCCCGTGGCCGGAAACCAG  >  NC_002947/3657463‑3657752

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