Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,836 C→T 17.6% intergenic (‑86/+43) 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,8360CT17.6% 68.3 / 9.0 35intergenic (‑86/+43)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base T (2/4);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

AGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAA  >  NC_002947/337689‑337983
                                                                                                                                                    |                                                                                                                                                   
aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCACTaca                                                                                                                                                    >  3:171641/1‑149 (MQ=255)
 ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGGCACATTCAAGGTTGGCGCTacac                                                                                                                                                   >  5:137658/1‑149 (MQ=255)
   aaCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCt                                                                                                                                                 <  5:28576/149‑1 (MQ=25)
     cTGGTCGAGGGCCTGGAAGATGGTGGTTAGTTGGGGGCTGTGGATGGGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACTCCCGCTCCGACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCtgt                                                                                                                                               <  5:151952/149‑1 (MQ=18)
               gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGGGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGc                                                                                                                                     <  1:235259/149‑1 (MQ=255)
                   ggAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGc                                                                                                                                 <  3:200570/149‑1 (MQ=255)
                                        ggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGACAAATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGAGTGCCCGCGAAGCAGTCAGTGAgg                                                                                                             >  7:249882/1‑148 (MQ=255)
                                           tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                                 >  8:253958/1‑109 (MQ=255)
                                           tgtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                                 <  7:253958/109‑1 (MQ=255)
                                            gtgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGGCCGCGAAGCAGTCAGTGAGGtct                                                                                                          >  7:153859/1‑147 (MQ=255)
                                                     gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                                            >  7:66759/1‑136 (MQ=255)
                                                     gtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTAAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGt                                                                                                            <  8:66759/136‑1 (MQ=255)
                                                         ggcatggcatGTATGTGTGGTGTTTTTGTGGGCCTCTTCGCGGGCAGGCCCGCTCCCACAGGTATAGCGCCACATTCAAGTTTGCCTCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                                           <  4:171641/149‑1 (MQ=37)
                                                           cATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCCCGCCCCCTCCCCCACGTTTAACGGCCCAGGCAACGTTGGCGCGAAACCTGTGGGCGCCGGCGTGCCCCGGCAAGAGTCAGGGAAGGCTCAGCGGCCGTgggctgc                                                                                         >  2:142652/1‑144 (MQ=255)
                                                               gCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGAGCCACATTCAAGGTTGGCGCTTCACCTGTGGGAGCGGGCCTGCCCGCGAAGCAGACAGGGCCGTCTCAGCGGCCTTGGCAGGCCTca                                                                                     >  3:39783/1‑148 (MQ=37)
                                                                cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               >  6:45228/1‑90 (MQ=255)
                                                                cAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               <  5:45228/90‑1 (MQ=255)
                                                                                                   ggaaCGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                                 <  8:249882/146‑1 (MQ=255)
                                                                                                     cACGCCCGCTCCCTCAGGTATTGCGCCACATTCATGGTGGGCGCTTCACGTGTGGGAGCGGGGGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                               <  4:39783/149‑1 (MQ=255)
                                                                                                        gCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGCGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                                            >  6:217000/1‑149 (MQ=255)
                                                                                                                ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  3:234911/1‑148 (MQ=255)
                                                                                                                ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  8:35833/1‑148 (MQ=255)
                                                                                                                 cacaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   <  5:218358/149‑1 (MQ=25)
                                                                                                                  acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  5:208137/1‑146 (MQ=25)
                                                                                                                  acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                               >  5:154393/1‑136 (MQ=255)
                                                                                                                  acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCAc                                               <  6:154393/136‑1 (MQ=255)
                                                                                                                   caGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  3:91219/1‑145 (MQ=25)
                                                                                                                   caGGTATAGCGCCCCATTCAAGGTGGCCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGGCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                                 <  6:137658/149‑1 (MQ=255)
                                                                                                                           gcgcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg                         >  3:88626/1‑149 (MQ=255)
                                                                                                                              ccACATTCACGGTGGGCGATACACGTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAg                      <  1:142652/149‑1 (MQ=255)
                                                                                                                                 acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa                    <  7:35833/149‑1 (MQ=255)
                                                                                                                                   aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAacac                 >  3:34456/1‑149 (MQ=255)
                                                                                                                                       aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGc               <  4:88626/148‑1 (MQ=255)
                                                                                                                                          tgtGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgc           <  4:91219/147‑1 (MQ=255)
                                                                                                                                           gtGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcg          <  5:217000/148‑1 (MQ=255)
                                                                                                                                              gcgcTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcggc        <  4:61694/148‑1 (MQ=255)
                                                                                                                                                gcTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcggcg       >  8:20276/1‑147 (MQ=255)
                                                                                                                                                    cacCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaagaa  <  4:34456/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAA  >  NC_002947/337689‑337983

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: