Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,507,604 G→C 14.3% A289P (GCG→CCG)  cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,507,6040GC14.3% 106.7 / 11.8 42A289P (GCG→CCG) cycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base C (3/3);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTC  >  NC_002947/5507460‑5507749
                                                                                                                                                |                                                                                                                                                  
ccGGGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTTCTGGCGCTGATCGCGATCATGGCCGTGACCGCGTGGCGCGATGTGGTGCACGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCTGTCGTGCcggcgg                                                                                                                                                 <  8:82240/148‑1 (MQ=255)
 ctcgGGCGACCAAATCGATCCCGATCCGTATCATCGTGTTATATGTCCTGGCGGTGATCGCGATCTTGGCCGTGACCCCGTGGCGCGATGTGTTGCCCGGAAAAATCCCGTTCGTTGATTTGTTCGTGCTGGCAGGCTTGCcggcggcg                                                                                                                                               <  7:167937/147‑1 (MQ=255)
     ggCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                           <  7:95025/149‑1 (MQ=255)
     ggCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                           <  7:19713/149‑1 (MQ=255)
        gATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGca                                                                                                                                        >  7:128860/1‑149 (MQ=255)
         aTCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGcct                                                                                                                                       <  4:84738/149‑3 (MQ=255)
         aTCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGcct                                                                                                                                       <  4:2832/149‑3 (MQ=255)
             aCTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCAGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatcat                                                                                                                                    >  5:80888/1‑148 (MQ=255)
               tCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCaa                                                                                                                                 >  6:134466/1‑149 (MQ=255)
                             tATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGa                                                                                                                   >  6:211720/1‑149 (MQ=255)
                               tcatcGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                  <  5:134466/148‑1 (MQ=255)
                                              tgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtcg                                                                                                  <  6:80888/149‑1 (MQ=255)
                                                          tCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                      <  2:25205/149‑1 (MQ=255)
                                                           cgcgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgt                                                                                                    <  3:213116/134‑1 (MQ=255)
                                                           cgcgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgt                                                                                                    >  4:213116/1‑134 (MQ=255)
                                                                  aTGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGG‑AAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                  <  1:236503/112‑1 (MQ=255)
                                                                  aTGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGG‑AAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                  >  2:236503/1‑112 (MQ=255)
                                                                                    cgcgATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                      <  6:14720/123‑1 (MQ=255)
                                                                                    cgcgATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                      >  5:14720/1‑123 (MQ=255)
                                                                                             gtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcggcg                                                                                                           >  3:84843/1‑93 (MQ=255)
                                                                                             gtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcggcg                                                                                                           <  4:84843/93‑1 (MQ=255)
                                                                                             gtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTc                                                   >  1:107108/1‑149 (MQ=255)
                                                                                                cccGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGACAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCg                                                  >  1:69889/1‑147 (MQ=255)
                                                                                                    gAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCATGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTTGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGaa                                            <  7:48564/149‑1 (MQ=255)
                                                                                                         aGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGa                                       <  1:188949/149‑1 (MQ=255)
                                                                                                                  gttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGccc                               >  2:84933/1‑147 (MQ=255)
                                                                                                                      aGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGcc                                >  4:59981/1‑142 (MQ=255)
                                                                                                                      aGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGcc                                <  3:59981/142‑1 (MQ=255)
                                                                                                                       gttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGcc                                                                      <  8:69286/104‑1 (MQ=255)
                                                                                                                       gttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGcc                                                                      >  7:69286/1‑104 (MQ=255)
                                                                                                                        ttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGc                                                 <  8:234521/124‑1 (MQ=255)
                                                                                                                        ttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGc                                                 >  7:234521/1‑124 (MQ=255)
                                                                                                                          gttCGTGCAGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTc                      <  8:95521/148‑1 (MQ=255)
                                                                                                                            tCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTc                      >  6:109181/1‑146 (MQ=255)
                                                                                                                               tGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAg                >  2:3814/1‑149 (MQ=255)
                                                                                                                                   ggCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGt            >  1:256451/1‑149 (MQ=255)
                                                                                                                                    gCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTc           >  4:38172/1‑149 (MQ=255)
                                                                                                                                     gatgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGaa                 <  3:181129/139‑1 (MQ=255)
                                                                                                                                     gatgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGaa                 >  4:181129/4‑142 (MQ=255)
                                                                                                                                     cAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACCGCGCCCAAGGCTTTCGAGAAGCTGTc           <  1:3814/149‑1 (MQ=255)
                                                                                                                                       tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTCCGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCGTTCGAGCATCTGTCGAg        >  6:201805/2‑149 (MQ=255)
                                                                                                                                       tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTCGAg        >  5:254715/2‑149 (MQ=255)
                                                                                                                                            gCCGGCGGCGGCAAGCATCATCAACTTCGTGGGGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGt    <  8:128860/148‑1 (MQ=255)
                                                                                                                                             ccGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTc  >  8:100403/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGA‑CGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTC  >  NC_002947/5507460‑5507749

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: ‑