Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,827 Δ1 bp 14.9% intergenic (‑77/+52) PP_0279 ← / ← PP_0280 hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947337,8270T.14.9% 120.2 / 18.7 47intergenic (‑77/+52)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base . (5/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACA  >  NC_002947/337685‑337967
                                                                                                                                               |                                                                                                                                             
gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgcgc                                                                                                                                         >  3:101387/1‑149 (MQ=255)
    aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTAAGGTGGGGGCTGTGGAGGTGGTTGGTAATGGCAAGGATGTGTTGTGGTTTTGAGGGCCTCTTCGCGGGCCCGCCCGCTCCCACAGGTATAGCGCAACATTCAAGGTTGGCGCTaca                                                                                                                                     >  5:312016/1‑149 (MQ=255)
         cTGGTCTAGGGCCTGGTAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGTCCTCTTCGCTGTCACTCACGCTCCCACAGGTATAGCGCCAAATTCAAGTTTGGAGCTACACCtgt                                                                                                                                <  6:312016/149‑1 (MQ=255)
                           gCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGCCACATTCAAGGTTGGCGCTAACCCTTTGGGAGCGGGCGTGCCcgc                                                                                                              >  8:284772/1‑149 (MQ=255)
                                       gTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagtcagt                                                                                                 >  1:43227/1‑149 (MQ=255)
                                        tGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                      <  2:251174/112‑1 (MQ=255)
                                        tGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACAATCAAGGTTGGCGCTac                                                                                                                                      >  1:251174/1‑112 (MQ=255)
                                                     ggtggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagtca                                                                                                   <  1:56281/133‑1 (MQ=255)
                                                     ggtggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagtca                                                                                                   >  2:56281/1‑133 (MQ=255)
                                                             ggCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggca                                                                            <  7:284772/148‑1 (MQ=255)
                                                              gCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTTGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTAAGTGAGGTCTCAGCGGCCTTggcagg                                                                          >  1:202563/1‑149 (MQ=255)
                                                                   gCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGGCACATTGAAGGGTGGCGCGACACCGGGGGGGGCGGGCGTGCCCGCGCAGCAGTCAGTGAGGtcaga                                                                                        >  4:72119/1‑127 (MQ=37)
                                                                   gCAAGGATGTGTGGGGGTTTTGAGGCCCTCTTCGCGGCCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCa                                                                                        <  3:72119/130‑1 (MQ=255)
                                                                     aaGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAc                                                                   >  5:118718/1‑149 (MQ=255)
                                                                      aGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCACCCACAGGTTTAGCGCCACCTTGAAGGTTGGGGCGACCCCAGGGGGAGCGGGCGTGCCCGCGAGGGAGGGAGTGAGGGGGCAGCGGGCGTGGCAGGCAGCCAcc                                                                  >  6:301649/1‑148 (MQ=19)
                                                                        gATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAcgcg                                                                <  2:43227/149‑1 (MQ=255)
                                                                            tgggggggtTTTGCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAATGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCaa                                                             <  1:62386/143‑1 (MQ=37)
                                                                                 tggtTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATTGCGCCACATTCAAGGTTGGCGGTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAgcg                                                       >  7:8283/1‑149 (MQ=255)
                                                                                       tgtgggCCTTTTCGCGGGCACGCCGGATCCCACTGTTTTAGGCCAAAATTCATGGTTGGCGCTAAACCTGTGGTTGGGGGGGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCca                                                  <  8:40350/145‑1 (MQ=255)
                                                                                           ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTTGCAGGCATCCCCGCGCACCCAGCGCTCCagcagc                                             >  2:243286/1‑149 (MQ=255)
                                                                                              ctctTCGCGGGCACGCCCCCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCAGTGGGAGCCGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCATGGCAGGCATCCACGCGCAACCAGCGCTCCaccagcat                                           >  6:118059/1‑146 (MQ=255)
                                                                                                 ttCGCGGGCACGCCCGCTCCCACAGGTATAGTGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                       <  5:135412/149‑1 (MQ=255)
                                                                                                    gcgGCCACGCCCGCACCCACAGGAATACCGCCACAATCAAGCTTGGCTCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGCTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCt                                    <  8:234815/149‑1 (MQ=38)
                                                                                                      gggCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  <  8:8283/149‑1 (MQ=255)
                                                                                                       ggCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                                  >  8:206347/1‑148 (MQ=18)
                                                                                                             cccGCTCCCACAGGTATTGCGCCACATTCAAGGTTGGCGCTACACGTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGc                           <  4:101387/149‑1 (MQ=255)
                                                                                                               cGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCAt                         <  7:206347/149‑1 (MQ=21)
                                                                                                                    ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCagcagc                                             <  1:6320/124‑1 (MQ=18)
                                                                                                                    ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCagcagc                                             >  2:6320/1‑124 (MQ=18)
                                                                                                                    ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  5:100059/1‑148 (MQ=25)
                                                                                                                    ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  7:289565/1‑148 (MQ=25)
                                                                                                                    ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGc                           >  5:57966/1‑142 (MQ=255)
                                                                                                                    ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGc                           <  6:57966/142‑1 (MQ=255)
                                                                                                                     cacaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  3:205972/1‑147 (MQ=255)
                                                                                                                      acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  4:116635/1‑146 (MQ=25)
                                                                                                                      acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     >  6:133529/1‑146 (MQ=255)
                                                                                                                          gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCa              >  5:104413/1‑149 (MQ=255)
                                                                                                                          gTATAGCGCCACATTAATGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGGCTCAGCGGCCATGGCAGGCATAAACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCa              >  3:103516/1‑149 (MQ=255)
                                                                                                                           tataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCtt                                   >  7:215454/1‑127 (MQ=255)
                                                                                                                           tataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCtt                                   <  8:215454/127‑1 (MQ=255)
                                                                                                                             taGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAcc            <  4:205972/148‑1 (MQ=255)
                                                                                                                                 gcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCt                                            <  3:48832/112‑1 (MQ=255)
                                                                                                                                 gcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCt                                            >  4:48832/1‑112 (MQ=255)
                                                                                                                                     acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa    <  5:133529/149‑1 (MQ=255)
                                                                                                                                       aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAaca  >  1:162625/1‑149 (MQ=255)
                                                                                                                                        tttgaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGTCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     <  6:237052/125‑1 (MQ=255)
                                                                                                                                        ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                     <  6:241028/129‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACA  >  NC_002947/337685‑337967

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