Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,661,297 Δ1 bp 10.5% coding (871/1701 nt) ligB → DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,661,2970G.10.5% 112.6 / 8.8 38coding (871/1701 nt)ligBDNA ligase B
Reads supporting (aligned to +/- strand):  ref base G (14/20);  new base . (2/2);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

GACTGGGTTTTCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGCCCGTTCTTCACGTGCAAC  >  NC_002947/5661151‑5661432
                                                                                                                                                  |                                                                                                                                        
gACTGGGTTTTCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAgc                                                                                                                                       >  1:39828/1‑150 (MQ=255)
     ggTTTTCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTgg                                                                                                                                 <  2:664718/150‑1 (MQ=255)
       ttttCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCa                                                                                                                               >  1:361677/1‑150 (MQ=255)
        tttCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAg                                                                                                                              <  1:140388/150‑1 (MQ=255)
                            cAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGAt                                                                                                          <  2:39828/150‑1 (MQ=255)
                            cAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCAAAGGCGCCATACTGGAt                                                                                                          <  1:509833/150‑1 (MQ=255)
                               cATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATcgc                                                                                                       >  1:163167/1‑150 (MQ=255)
                                aTCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATcgag                                                                                                      >  1:657112/1‑148 (MQ=255)
                                  cGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGc                                                                                                    <  2:361677/150‑1 (MQ=255)
                                   gCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCAGAGAAGGACAAGCCCTGCCGCGGGGGACGCGGGCGCCTGCCGAGCGCTGGCAGGCCAAGGCGCCATAATGGATGGCGGcc                                                                                                   >  1:379487/1‑150 (MQ=255)
                                    ccATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTgg                                                                                                                                 <  2:763011/119‑1 (MQ=255)
                                    ccATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCACTCGGCCACCTGCCGAGCGCTgg                                                                                                                                 >  1:763011/1‑119 (MQ=255)
                                         gACACGCTGGAGGAAGTAGCACACTGGCGCGAGCACTGGTCTTGATATGCGGGGGTGGTGGTGAGCGTCGTTGTGCTCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAg                                                                                             <  1:692442/150‑1 (MQ=255)
                                                tGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCt                                                                                      <  2:152460/150‑1 (MQ=255)
                                                tGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCt                                                                                      <  1:204914/150‑1 (MQ=255)
                                                  gaggaTGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTAt                                                                                    <  2:512095/150‑1 (MQ=255)
                                                          agcaCACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGc                                                                            >  1:263911/1‑150 (MQ=255)
                                                           gcaCACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGcc                                                                           >  2:76407/1‑150 (MQ=255)
                                                                   gcgcGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGa                                                                   >  2:644832/1‑150 (MQ=255)
                                                                    cgcgAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGaa                                                                  <  2:501742/150‑1 (MQ=255)
                                                                             tGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGAc                                                         >  2:736288/1‑150 (MQ=255)
                                                                               gTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGt                                                       >  1:45577/1‑150 (MQ=255)
                                                                                  tCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATcgcg                                                                                                      <  2:490541/100‑1 (MQ=255)
                                                                                  tCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATcgcg                                                                                                      >  1:490541/1‑100 (MQ=255)
                                                                                         gctctgcCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCAAAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCgg                                             <  2:379487/150‑1 (MQ=255)
                                                                                                 cTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCg                                     <  1:76407/150‑1 (MQ=255)
                                                                                                      ccACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTg                                <  2:680025/150‑1 (MQ=255)
                                                                                                        aCCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGc                              <  2:163167/150‑1 (MQ=255)
                                                                                                                    aTCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGAGCGGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGccc                  >  1:461621/1‑150 (MQ=255)
                                                                                                                           gcCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTCCGGGTCGGCCGCACTGGCAGGGTCACGCCCGttcttc           <  1:644832/150‑1 (MQ=255)
                                                                                                                           gcCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGAGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGCCCGttcttc           <  2:657112/150‑1 (MQ=255)
                                                                                                                               gggCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGcc                   <  2:283814/138‑1 (MQ=255)
                                                                                                                               gggCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGcc                   >  1:283814/1‑138 (MQ=255)
                                                                                                                                    gTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGCCCGTTCTTCACGTGCAAc  <  1:736288/150‑1 (MQ=255)
                                                                                                                                        ttggccTGCCAGCGCTCGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAggg                          <  1:521503/118‑1 (MQ=255)
                                                                                                                                        ttggccTGCCAGCGCTCGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAggg                          >  2:521503/5‑122 (MQ=255)
                                                                                                                                        ttggccTGCCAGCGCTCGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAggg                          >  2:521544/5‑122 (MQ=255)
                                                                                                                                        gtgtccTGCCAGCGCTCGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAggg                          <  1:521544/118‑1 (MQ=255)
                                                                                                                                                  gAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAg                                                                             >  2:792480/1‑61 (MQ=255)
                                                                                                                                                  gAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAg                                                                             <  1:792480/61‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GACTGGGTTTTCCGGACAGCCAGCTTTACAGCATCGCCATCGACACGCTGGAGGATGCAGCACACTGGCGCGAGCACTGGTATCGCTCTGCTCTGCCCTTCGCCACCGACGGTGTGATCCTGCGCCAGGGCAGTCGGCCACCTGCCGAGCGCTGGCAGGCCAAGGCGCCATACTGGATCGCGGCCTGGAAGTACCCCTATGCACAAGCCTTGGCCGAAGTACGCGACGTGCGCTTTCGGGTCGGCCGCACTGGCAGGGTCACGCCCGTTCTTCACGTGCAAC  >  NC_002947/5661151‑5661432

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