Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,263,488 C→T 7.5% A307T (GCC→ACC)  PP_1105 ← ATP‑dependent DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,263,4880CT 7.5% 152.2 / 2.5 53A307T (GCC→ACC) PP_1105ATP‑dependent DNA ligase
Reads supporting (aligned to +/- strand):  ref base C (20/29);  new base T (2/2);  total (22/31)
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.52e-01

CCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGTTCTTCGCCGCGTGA  >  NC_002947/1263341‑1263632
                                                                                                                                                   |                                                                                                                                                
ccAGCTGGCTTCGGCGTTCGTGCTGCTGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGc                                                                                                                                                  <  4:258778/148‑1 (MQ=255)
        cTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTg                                                                                                                                         <  4:202292/149‑1 (MQ=255)
                   gTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGtt                                                                                                                              <  3:255892/149‑1 (MQ=255)
                     gctgcGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATAGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTcc                                                                                                                            >  8:204217/1‑149 (MQ=255)
                      ctgcGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTcc                                                                                                                            <  3:57946/148‑1 (MQ=255)
                      ctgcGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTcc                                                                                                                            <  3:62550/148‑1 (MQ=255)
                           gTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCg                                                                                                                      <  6:341869/149‑1 (MQ=255)
                              ggCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGt                                                                                                                    <  1:18917/148‑1 (MQ=255)
                                     ccAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCgg                                                                                                            <  4:118740/149‑1 (MQ=255)
                                     ccAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCgg                                                                                                            <  6:61235/149‑1 (MQ=255)
                                           cTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCa                                                                                                                >  4:133047/1‑139 (MQ=255)
                                           cTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCa                                                                                                                <  3:133047/139‑1 (MQ=255)
                                           cTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTc                                                                                                       >  6:309421/1‑148 (MQ=255)
                                            tCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTc                                                                                                       >  7:230133/1‑147 (MQ=255)
                                                       cTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGAtct                                                                                          >  3:209364/1‑149 (MQ=255)
                                                          cagcagGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGAc                                                                                                  <  2:161436/138‑1 (MQ=255)
                                                          cagcagGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGAc                                                                                                  >  1:161436/1‑138 (MQ=255)
                                                          cagcagGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGAtctc                                                                                         <  6:191689/147‑1 (MQ=255)
                                                          cagcagGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCAAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCACTTCCGTCTCGGTTTCACCGCGTTTCCAGACCACGATCTCGc                                                                                       >  2:160747/1‑149 (MQ=255)
                                                             cagGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCAt                                                                                    <  7:213313/149‑1 (MQ=255)
                                                                gTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCAt                                                                                    <  6:298704/146‑1 (MQ=255)
                                                                gTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATc                                                                                   >  2:219640/1‑147 (MQ=255)
                                                                gTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGa                                                                                 <  1:140656/149‑1 (MQ=255)
                                                                gTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGACTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCAt                                                                                    >  5:298704/1‑146 (MQ=255)
                                                                    tAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGcac                                                                             <  5:309421/149‑1 (MQ=255)
                                                                       gCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGcaccac                                                                          >  2:282240/1‑149 (MQ=255)
                                                                       gCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGcaccac                                                                          <  1:219640/149‑1 (MQ=255)
                                                                               caccacCGGCTCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCgg                                                                  <  8:180256/149‑1 (MQ=255)
                                                                                accaccGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCggg                                                                 <  4:47817/149‑1 (MQ=255)
                                                                                  caccGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTc                                                                                                       >  2:7003/1‑109 (MQ=255)
                                                                                  caccGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTc                                                                                                       <  1:7003/109‑1 (MQ=255)
                                                                                                   gcagTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCa                                              <  7:204217/149‑1 (MQ=255)
                                                                                                    cagTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGTCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCAt                                                                     <  1:58502/125‑1 (MQ=255)
                                                                                                    cagTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCAt                                                                     >  2:58502/1‑125 (MQ=255)
                                                                                                    cagTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGAAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCAt                                             >  2:226218/1‑149 (MQ=255)
                                                                                                    cagTTTTTAGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCAt                                             >  5:73496/1‑149 (MQ=255)
                                                                                                     agTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATg                                            <  6:212098/149‑1 (MQ=255)
                                                                                                              cccAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCggg                                   >  7:245140/1‑149 (MQ=255)
                                                                                                              cccAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgg                                    <  1:160747/148‑1 (MQ=255)
                                                                                                                cAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgggg                                  >  3:161198/1‑148 (MQ=255)
                                                                                                                cAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgggg                                  >  1:165550/1‑148 (MQ=255)
                                                                                                                           gACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTa                      <  2:165550/149‑1 (MQ=255)
                                                                                                                            aCCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTaa                     >  7:25259/1‑149 (MQ=255)
                                                                                                                              cGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGGCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGGCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAAcc                   <  6:73496/149‑1 (MQ=255)
                                                                                                                              cGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGGTTGACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACTCCATCGGGCAACGTCTGTTCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAAcc                   <  4:209364/149‑1 (MQ=255)
                                                                                                                                aTGCGCTGTTGCAGCAAGGTGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgg                                    >  4:309696/1‑130 (MQ=255)
                                                                                                                                aTGCGCTGTTGCAGCAAGGTGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgg                                    >  4:304651/1‑130 (MQ=255)
                                                                                                                                aTGCGCTGTTGCAGCAAGGTGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgg                                    <  3:309696/130‑1 (MQ=255)
                                                                                                                                aTGCGCTGTTGCAGCAAGGTGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCgg                                    <  3:304651/130‑1 (MQ=255)
                                                                                                                                     cTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGttctt            >  8:100631/1‑149 (MQ=255)
                                                                                                                                        ttGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGttcttc           <  7:100631/147‑1 (MQ=255)
                                                                                                                                           cagcaAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGTTCTTcgccgc      >  3:82544/1‑149 (MQ=255)
                                                                                                                                             gcaAGTCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGTCCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGTTCTTCGCCGCGt    <  1:282240/149‑1 (MQ=255)
                                                                                                                                               aaGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGTTCTTCGCCGCGTGa  <  8:25259/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATGCGCTGTTGCAGCAAGGCGAAAGGCTGCACGGCCAGTTCCGTCTCGGTTTCACCGGGTTTCCAGACCAGGATCTCGCCATCGAGCACCACGCCATCGGGCAACGTCTGTGCCAGGCCATGCAGTTCGGGGAAGCGCTCGGTAACCAGTTCTTCGCCGCGTGA  >  NC_002947/1263341‑1263632

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