Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,429,808 C→A 8.0% E24* (GAG→TAG) ‡ PP_4768 ← exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,429,8080CA 8.0% 154.0 / 3.0 50E24* (GAG→TAG) ‡PP_4768exonuclease
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base A (2/2);  total (26/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 = 6.06e-01

TTCTCCACCGGCGGCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGCC  >  NC_002947/5429663‑5429956
                                                                                                                                                 |                                                                                                                                                    
ttCTCCACCGGCGGCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTg                                                                                                                                                  <  2:141070/150‑1 (MQ=255)
    ccACCGGCGGCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCGAGGTGGcgg                                                                                                                                              <  2:724533/150‑3 (MQ=255)
     cACCGGCGGCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTCGCGGCTGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACg                                                                                                                                             >  1:379512/1‑150 (MQ=255)
        cggcggCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCa                                                                                                                                          <  1:467915/150‑1 (MQ=255)
              cgcACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGt                                                                                                                                    <  2:262929/150‑1 (MQ=255)
                       cagcaTGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCCCCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGa                                                                                                                           >  1:854542/1‑150 (MQ=255)
                       cagcaTGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGa                                                                                                                           >  1:150859/1‑150 (MQ=255)
                             ggCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGAAGGAACCGGCAGGCTAGTATTCAGCAGGCTTTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATAcc                                                                                                                     >  2:740791/1‑150 (MQ=255)
                               cGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCgg                                                                                                                   >  2:1082688/1‑150 (MQ=255)
                                ggtggtAATGCCGGTCAGCCCTGCGACGAAGGCTGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCggt                                                                                                                  <  2:307703/150‑1 (MQ=255)
                                 gtggtAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCggt                                                                                                                  <  2:75824/149‑1 (MQ=255)
                                 gtggtAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCggt                                                                                                                  >  1:75824/1‑149 (MQ=255)
                                     tAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACGGCCACCTCGGTGCAACGGCAGCCGGTACCGGGGGAAATACCGGTGGttt                                                                                                             >  2:528682/1‑150 (MQ=255)
                                              tCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGa                                                                                                    <  1:540742/150‑1 (MQ=255)
                                                aGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCGAGGTGGcggtg                                                                                                                                            <  2:653010/108‑5 (MQ=255)
                                                aGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCGAGGTGGcggtg                                                                                                                                            >  1:653010/1‑104 (MQ=255)
                                                 gCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGa                                                                                                 <  1:949369/150‑1 (MQ=255)
                                                    cTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCgggg                                                                                                                             <  2:903475/119‑1 (MQ=255)
                                                    cTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCgggg                                                                                                                             >  1:903475/1‑119 (MQ=255)
                                                     tGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGaa                                                                                                                          >  2:359353/1‑121 (MQ=255)
                                                     tGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGaa                                                                                                                          <  1:359353/121‑1 (MQ=255)
                                                        gacgaAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGAt                                                                                          >  2:159293/1‑150 (MQ=255)
                                                                ccggcaccggcaCGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCa                                                                                  <  2:1064557/150‑1 (MQ=255)
                                                                 cggcaccggcacGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCAt                                                                                 >  2:982831/1‑150 (MQ=255)
                                                                       cggcacGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAg                                                                           >  1:1078602/1‑150 (MQ=255)
                                                                               cAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTc                                                                                                      >  2:656610/1‑115 (MQ=255)
                                                                               cAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTc                                                                                                      <  1:656610/115‑1 (MQ=255)
                                                                               cAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAgt                                                                   <  1:159293/150‑1 (MQ=255)
                                                                                   aTTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGttt                                                               >  2:1049282/1‑150 (MQ=255)
                                                                                             ctctGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAgc                                                     >  2:830649/1‑150 (MQ=255)
                                                                                               ctGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAgcgc                                                   <  2:1078602/150‑1 (MQ=255)
                                                                                                  gTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGcc                                                >  2:416582/1‑150 (MQ=255)
                                                                                                       gcTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCCCGCCTTCGa                                           >  2:1052180/1‑150 (MQ=255)
                                                                                                        cTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGAt                                          <  1:830649/150‑1 (MQ=255)
                                                                                                              aaTGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATTGGCAc                                    >  2:641740/1‑150 (MQ=255)
                                                                                                              aaTGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCAc                                    <  1:740791/150‑1 (MQ=255)
                                                                                                                       aCCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGa                                           <  1:783123/134‑1 (MQ=255)
                                                                                                                       aCCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGt                           <  1:416582/150‑1 (MQ=255)
                                                                                                                       aCCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAAGACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGa                                           >  2:783123/1‑134 (MQ=255)
                                                                                                                             aaGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGgcagc                     >  1:803978/1‑150 (MQ=255)
                                                                                                                                  ttACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGGCAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGa                <  1:8861/150‑1 (MQ=255)
                                                                                                                                  ttACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACTTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGa                <  1:364797/150‑1 (MQ=255)
                                                                                                                                     ccaccGCCACCTCGGTGGCACGGCAGCCGGTACCGCGGGAAATACCGGAGGTTTCAAAGTCGATGACAGCGGTACGTTCCATGAAAAGACCTTATGGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCa             >  2:440692/1‑150 (MQ=255)
                                                                                                                                      caccGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGtt                                                                                                              <  2:1122537/52‑1 (MQ=255)
                                                                                                                                      caccGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGtt                                                                                                              >  1:1122537/1‑52 (MQ=255)
                                                                                                                                       cgtgCCACCGAGGGGGCACGGCAGGCGGTACCGGGGGCAATACCGGTGGTTTCAAAGTCGACGACAGCGATACGTTCCTGGTATAGACCTTAGTGTATGCGCAACAGCGCGCCGTCGATGGTGACGTAACGGGTGGTAGCACGGATctgg           >  2:212728/4‑147 (MQ=255)
                                                                                                                                       cccGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCGTTGAAAAGACCGTAGGGTTTTAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGt           <  1:1052180/149‑1 (MQ=255)
                                                                                                                                             aCCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGTGACAAc     >  1:715349/1‑150 (MQ=255)
                                                                                                                                             aCCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGTGACAAc     <  2:803978/150‑1 (MQ=255)
                                                                                                                                                tCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGcc  >  1:1148570/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCTCCACCGGCGGCGCACTGCGCAGCATGGCGGTGGTAATGCCGGTCAGCCCTGCGACGAAGGCCGGCACCGGCACGCCAGCATTCATCAGGCTCTGGTAGCGCTCGACAATGCGCCCACCCTCAAGCATTACCACCGCCACCTCGGTGGCACGGCAGCCGGTACCGGGGGAAATACCGGTGGTTTCAAAGTCGATGACAGCGATACGTTCCATGAAAAGACCTTAGTGTTTGAGCAACAGCGCGCCTTCGATGGGCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGCC  >  NC_002947/5429663‑5429956

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