Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,816:1 +A 25.0% intergenic (‑66/+63) 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,8161.A25.0% 52.1 / 17.9 20intergenic (‑66/+63)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (6/9);  new base A (3/2);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

CTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAC  >  NC_002947/337676‑337948
                                                                                                                                             |                                                                                                                                    
ctgGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGTTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCa                                                                                                                               <  8:49170/148‑1 (MQ=255)
ctgGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCa                                                                                                                               <  6:161143/148‑1 (MQ=255)
    tGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGtt                                                                                                                          <  5:348490/149‑1 (MQ=255)
         gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGTTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgcgc                                                                                                                     <  4:212012/149‑1 (MQ=255)
           tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACACGTATAGCGCCCCATTCAAGGTTGgcgcaa                                                                                                                   >  5:288411/1‑147 (MQ=255)
              ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                <  7:200765/149‑1 (MQ=255)
              ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                <  6:246900/149‑1 (MQ=255)
                               tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                       >  5:176049/1‑125 (MQ=255)
                               tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                       <  6:176049/125‑1 (MQ=255)
                                 gAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                <  6:36184/146‑1 (MQ=255)
                                 gAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                >  5:36184/1‑146 (MQ=255)
                                                     ggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTgg                                                                                                          >  6:369833/1‑116 (MQ=17)
                                                     ggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTgg                                                                                                          <  5:369833/116‑1 (MQ=17)
                                                                                                                             ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAc                                                <  1:72091/102‑1 (MQ=14)
                                                                                                                             ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAc                                                >  2:72091/1‑102 (MQ=14)
                                                                                                                             ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc  >  2:116416/1‑148 (MQ=255)
                                                                                                                               acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc  >  2:105201/1‑146 (MQ=255)
                                                                                                                                 aGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAGGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc  >  2:120635/1‑144 (MQ=25)
                                                                                                                                 aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTc                                 >  1:175885/1‑113 (MQ=255)
                                                                                                                                 aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTc                                 <  2:175885/113‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                    
CTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAC  >  NC_002947/337676‑337948

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: