Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 922,569 C→T 15.4% L154L (CTG→CTA PP_0803 ← protein secretion ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947922,5690CT15.4% 62.5 / 4.2 26L154L (CTG→CTAPP_0803protein secretion ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (8/14);  new base T (2/2);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

CGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGC  >  NC_002947/922444‑922703
                                                                                                                             |                                                                                                                                      
ctcTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcggc                                                                                                                              <  1:345803/134‑1 (MQ=255)
cGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcggc                                                                                                                              <  2:59258/136‑1 (MQ=255)
cGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcggc                                                                                                                              >  2:305829/1‑136 (MQ=255)
cGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcggc                                                                                                                              >  2:173325/1‑136 (MQ=255)
cGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcggc                                                                                                                              <  1:165646/136‑1 (MQ=255)
                  aTCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAAcg                                                                                                            >  1:282798/1‑136 (MQ=255)
                     aCTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCgg                                                                                                         <  1:472094/136‑1 (MQ=255)
                          aaCCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCg                                                                                                    <  1:306369/136‑1 (MQ=255)
                          aaCCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCg                                                                                                    <  1:253667/136‑1 (MQ=255)
                               ggCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATc                                                                                               <  1:305829/136‑1 (MQ=255)
                               ggCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATc                                                                                               >  1:426043/1‑136 (MQ=255)
                               ggCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATc                                                                                               >  1:196667/1‑136 (MQ=255)
                               ggCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATc                                                                                               <  1:173325/136‑1 (MQ=255)
                               ggCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCACCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATc                                                                                               <  2:392382/136‑1 (MQ=255)
                                                     cAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCgcgc                                                                         <  2:282798/136‑1 (MQ=255)
                                                            tcGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCTAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGctct                                                                  >  2:416096/1‑136 (MQ=255)
                                                            tcGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCTAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGctct                                                                  >  2:386794/1‑136 (MQ=255)
                                                                                               cGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcgg                                                                                                                               >  2:72945/1‑40 (MQ=255)
                                                                                               cGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCgcgg                                                                                                                               <  1:72945/40‑1 (MQ=255)
                                                                                                  gTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTctgcctgc                            <  2:426043/136‑1 (MQ=255)
                                                                                                      tGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTAc                                                                                                                >  1:498930/1‑48 (MQ=255)
                                                                                                      tGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTAc                                                                                                                <  2:498930/48‑1 (MQ=255)
                                                                                                       gAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCg                       <  2:287192/136‑1 (MQ=255)
                                                                                                           aCCCGCCAGGGTCTGTTCTAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGcc                   <  1:416096/136‑1 (MQ=255)
                                                                                                           aCCCGCCAGGGTCTGTTCTAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGcc                   <  1:386794/136‑1 (MQ=255)
                                                                                                                ccAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCctggct              >  1:31927/1‑136 (MQ=255)
                                                                                                                            ccAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAgc  >  2:48926/1‑136 (MQ=255)
                                                                                                                             |                                                                                                                                      
CGCTGGCGGTTCAGGCGGATCACTTCAACCCGGCTGGAGGCACCCATCTTTGCCAGGTTCTCGGTGATCTTCAGCTCGCTGCGCACCAGCTGCAACGAGTCCTGAATACCCGCCAGGGTCTGTTCCAGGCCGCGGCGGCGGGTCTGGTACAACGCGGTTTCGGCATCGATCAGGTCAGGAGAGTCGCGCAGGCTCTCCGGGAAGGTCAGTGGCTTGCCGGTCACCTCTGCCTGCAGGCGGGCCTGGCTGGCCCTGGCAGC  >  NC_002947/922444‑922703

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: