Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,264,180 T→C 8.7% Y76C (TAT→TGT)  PP_1105 ← ATP‑dependent DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,264,1800TC 8.7% 126.8 / 3.0 46Y76C (TAT→TGT) PP_1105ATP‑dependent DNA ligase
Reads supporting (aligned to +/- strand):  ref base T (20/22);  new base C (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCGCTGGCGGACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGGCGCCATGTCGAA  >  NC_002947/1264038‑1264325
                                                                                                                                              |                                                                                                                                                 
gcgcTGGCGGACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAACGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCttt                                                                                                                                             >  7:228/1‑149 (MQ=255)
   cTGGCGGACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAACGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTttctt                                                                                                                                          <  8:228/149‑1 (MQ=255)
        ggACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGc                                                                                                                                                <  4:217829/138‑1 (MQ=255)
        ggACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGc                                                                                                                                                >  3:217829/1‑138 (MQ=255)
                tCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGTAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTc                                                                                                                             >  2:155809/1‑149 (MQ=255)
                  aGGCGTCAGCCCGCGCAGGGGTTGCAAGTGCTGTTCCACGCAGGTGGCCAGGCCTTCCTCGTTGCCATGGGCGTGCAGCGGCAGCAGCAGCGTAATGGTCTCGGCGAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTcc                                                                                                                            <  2:188263/148‑1 (MQ=255)
                      gTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAg                                                                                                                       >  6:234398/1‑149 (MQ=255)
                          gCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCACACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGgcca                                                                                                                   <  6:243905/149‑1 (MQ=255)
                           cccGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGgccag                                                                                                                  <  3:55615/149‑1 (MQ=255)
                                gcAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGTAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCcagccagc                                                                                                             <  1:155809/149‑1 (MQ=255)
                                         gCATGTGCTGTTCGACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGACAGGTTTCTTCGAAAAGCGAGTCCGGCAGGCCAGCCAGCGTCGgcgc                                                                                                     <  7:162127/148‑1 (MQ=255)
                                         gCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGTCCTGACAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGc                                                                                                     <  7:157598/148‑1 (MQ=255)
                                         gCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGACAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGc                                                                                                     >  8:157598/1‑148 (MQ=255)
                                         gCAAGTGCTGTTCCACCCAGGTGCCCACGCCTTCCTCGTTGCCATCGGCGTCCAGCGGCAGCAGCAGCCAAATGGTCTCGGCCAGGTCGCCAACGGCATGACAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCACGCCAGCGACGTAGcc                                                                                                    >  8:162127/1‑149 (MQ=255)
                                          cAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCa                                                                                                   <  8:83123/149‑1 (MQ=255)
                                                       aCCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAg                                                                                      >  1:291805/1‑149 (MQ=255)
                                                        cccAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCgg                                                                                                                          <  8:128791/112‑1 (MQ=255)
                                                        cccAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCgg                                                                                                                          >  7:128791/1‑112 (MQ=255)
                                                                     gCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCcagccagc                                                                                                             >  6:39913/1‑112 (MQ=255)
                                                                     gCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCcagccagc                                                                                                             <  5:39913/112‑1 (MQ=255)
                                                                        ttCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCt                                                                     >  4:101347/1‑149 (MQ=255)
                                                                          ccTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTgg                                                                   >  6:138932/1‑149 (MQ=255)
                                                                           cTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGc                                                                  >  4:99833/1‑149 (MQ=255)
                                                                              gTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTg                                                               >  4:174640/1‑149 (MQ=255)
                                                                               ttGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTgg                                                              <  3:99833/149‑1 (MQ=255)
                                                                               ttGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTgg                                                              <  5:234398/149‑1 (MQ=255)
                                                                                   cATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCCGCCAGCGTCGTCGCCAGTGCGCGGAGAAGGCGGGTGGGCCCCAGCTGGCGTGGgcgc                                                          >  8:4481/1‑149 (MQ=255)
                                                                                   cATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGgcgc                                                          <  3:101347/149‑1 (MQ=255)
                                                                                     tCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGcgccc                                                        >  5:126120/1‑149 (MQ=255)
                                                                                       ggCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGcgcccgc                                                      >  3:72493/1‑149 (MQ=255)
                                                                                        gCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGcgcccgcc                                                     >  1:46768/1‑149 (MQ=255)
                                                                                        gCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGcgcccgcc                                                     <  2:291805/149‑1 (MQ=255)
                                                                                             ccgCGGCAGCAGCAGCGAACTGGGCTCGGCCAGGGCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCa                                                <  7:4481/147‑1 (MQ=255)
                                                                                                  gcagcagcaGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGaaa                                           <  2:19207/149‑1 (MQ=255)
                                                                                                      cagcagcGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATaca                                       <  3:174640/149‑1 (MQ=255)
                                                                                                       agcagcGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCATCGTCGTCGCCAGTTCGCGTAGAAGGCTGGTGGGCACCCGCTGGCGTGGGCGCCCGCCGGCCAGGAAATacac                                      >  4:21165/1‑149 (MQ=255)
                                                                                                             gAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAACCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTCGCGTGGGCGCCCGCCGGCCAGGAAATACACCGcc                                  >  7:83415/1‑147 (MQ=255)
                                                                                                                 tGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGgca                            <  4:72493/149‑1 (MQ=255)
                                                                                                                              gTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTcggcgg               >  3:189822/1‑149 (MQ=255)
                                                                                                                                 gCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCggcgggcg            >  3:300630/1‑149 (MQ=255)
                                                                                                                                 gCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCggcgggc             <  4:189822/148‑1 (MQ=255)
                                                                                                                                    cACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGgcgc           >  7:228504/1‑147 (MQ=255)
                                                                                                                                    cACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGgcgc           <  8:228504/147‑1 (MQ=255)
                                                                                                                                             aTAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGGCGCCATGTCGaa  <  5:138932/147‑1 (MQ=255)
                                                                                                                                             aTAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGGCGCCATGTCGaa  <  6:126120/147‑1 (MQ=255)
                                                                                                                                             aTAGCTCTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGGCGCCATGTCGaa  <  1:16103/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCGCTGGCGGACTTGCTCAGGCGTCAGCCCGCGCAGGGGTAGCAAGTGCTGTTCCACCCAGGTGGCCAGGCCTTCCTCGTTGCCATCGGCGTGCAGCGGCAGCAGCAGCGAAATGGTCTCGGCCAGGTCGCCCACGGCCTGATAGCTTTCTTCGAAAAGCCAGTCCGGCAGGCCAGCCAGCGTCGTCGCCAGTTCGCGTAGAAGGCGGGTGGGCACCAGCTGGCGTGGGCGCCCGCCGGCCAGGAAATACACCGCCCAGGCAGCATCTTCGGCGGGCGCCATGTCGAA  >  NC_002947/1264038‑1264325

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