Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,078,315 G→T 19.9% L316M (CTG→ATG)  PP_5327 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,078,3150GT19.9% 54.8 / 2.6 20L316M (CTG→ATG) PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (10/6);  new base T (2/2);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCAGGAACACCGGCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACA  >  NC_002947/6078188‑6078440
                                                                                                                               |                                                                                                                             
gCAGGAACACCGGCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGtcatca                                                                                                                       >  1:36453/1‑136 (MQ=255)
gCAGGAACACCGGCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGtcatca                                                                                                                       <  1:124249/136‑1 (MQ=255)
gCAGGAACACCGGCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGtcatca                                                                                                                       <  1:475204/136‑1 (MQ=255)
            gCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGa                                                                                                           <  2:36453/136‑1 (MQ=255)
             ccAGGTTGTTAACCGCAATGCGGATACGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGaa                                                                                                          >  2:392830/1‑136 (MQ=255)
               aGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGa                                                                                                        >  2:387206/1‑136 (MQ=255)
                                 cGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGc                                                                                      <  1:298838/136‑1 (MQ=255)
                                          cGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCATGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGgtgt                                                                             >  2:464588/1‑136 (MQ=255)
                                          cGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCATGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGgtgt                                                                             >  2:19312/1‑136 (MQ=255)
                                              tGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGc                                                                         >  2:517137/1‑136 (MQ=255)
                                                aCCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCtt                                                                       >  1:212514/1‑136 (MQ=255)
                                                     gCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCAcgc                                                                  <  2:459321/136‑1 (MQ=255)
                                                                 gTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCATGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTg                                                      <  1:464588/136‑1 (MQ=255)
                                                                 gTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCATGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTg                                                      <  1:19312/136‑1 (MQ=255)
                                                                   aTTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAg                                                    >  1:409143/1‑136 (MQ=255)
                                                                   aTTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAg                                                    >  1:40847/1‑136 (MQ=255)
                                                                       aCGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAAc                                                >  2:6796/1‑136 (MQ=255)
                                                                                                  gAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTg                                                      <  1:505671/103‑1 (MQ=255)
                                                                                                  gAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTg                                                      >  2:505671/1‑103 (MQ=255)
                                                                                                                     gCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACa  >  1:178510/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                             
GCAGGAACACCGGCCAGGTTGTTAACCGCAATGCGGATAAGCCGGGTGACCGGGCCCTGCCTGGCGTATTCACGCAGGTACACGGCAGCCAGCACGCCGAACGGGGTAACGATCACGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACA  >  NC_002947/6078188‑6078440

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