Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,730,938 C→G 10.9% P78A (CCG→GCG)  polB → DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,730,9380CG10.9% 137.3 / 7.5 46P78A (CCG→GCG) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base G (2/3);  total (21/25)
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.95e-01

GATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCA  >  NC_002947/2730791‑2731083
                                                                                                                                                   |                                                                                                                                                 
gATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCg                                                                                                                                                 >  1:901270/1‑150 (MQ=255)
    ccGGGCCGCGACTGTTGCGGCTGGCCAGGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGc                                                                                                                                             >  1:38447/1‑150 (MQ=255)
      gggCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                           <  1:429106/150‑1 (MQ=255)
         ccGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                        <  1:1133348/150‑3 (MQ=255)
                  ttGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCg                                                                                                                                          >  2:838776/1‑139 (MQ=255)
                  ttGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCg                                                                                                                                          <  1:838776/139‑1 (MQ=255)
                    gCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCg                                                                                                                             >  1:164903/1‑150 (MQ=255)
                             ggcgCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCg                                                                                                                    >  2:1051428/1‑150 (MQ=255)
                              gcgcAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGc                                                                                                                   >  1:593602/1‑150 (MQ=255)
                                gcAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGGGCGATGCTCgaa                                                                                                                                        >  1:1228859/1‑125 (MQ=255)
                                gcAGGAGTCGGTGGCTTTCATTCATCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGGGCGATGCTCgaa                                                                                                                                        <  2:1228859/127‑3 (MQ=255)
                                   ggAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAAtt                                                                                                              <  2:562482/150‑1 (MQ=255)
                                    gAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTg                                                                                                             >  2:1234042/1‑150 (MQ=255)
                                      gTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                           <  1:920955/150‑1 (MQ=255)
                                      gTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                           <  1:278520/150‑1 (MQ=255)
                                              cTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTgg                                                                                                   <  2:901270/150‑1 (MQ=255)
                                                       cTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGa                                                                                                            <  1:665206/132‑1 (MQ=255)
                                                       cTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGa                                                                                                            >  2:665206/1‑132 (MQ=255)
                                                                aGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACgg                                                                                 <  1:1427218/150‑1 (MQ=255)
                                                                 ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACgg                                                                                 <  1:803734/149‑1 (MQ=255)
                                                                 ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACgg                                                                                 >  2:803734/1‑149 (MQ=255)
                                                                   ccGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                              <  2:74440/150‑1 (MQ=255)
                                                                   ccGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                              <  2:593602/150‑1 (MQ=255)
                                                                     gAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGt                                                                            <  1:99754/150‑1 (MQ=255)
                                                                             ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                           <  1:733949/79‑1 (MQ=255)
                                                                             ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                           >  2:733949/1‑79 (MQ=255)
                                                                             ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                    >  1:15774/1‑150 (MQ=255)
                                                                                    ctgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAg                                                             <  1:280982/150‑1 (MQ=255)
                                                                                     tgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                              >  2:1248602/1‑148 (MQ=255)
                                                                                     tgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                              <  1:1248602/148‑1 (MQ=255)
                                                                                          gCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGat                                                       <  2:164903/150‑1 (MQ=255)
                                                                                               cGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCg                                                  <  1:1234042/150‑1 (MQ=255)
                                                                                                          ggAAGCTCCAGCCACTGCACCTTAAGGACTTCTAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGc                                       <  2:660822/149‑1 (MQ=255)
                                                                                                          cGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCACTTGAATGAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGGGTTCGAGGCCGATATCCGGCCGCCGGAGc                                       >  1:1315026/1‑150 (MQ=255)
                                                                                                                  aaGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTc                               >  2:1383805/1‑150 (MQ=255)
                                                                                                                      cACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTcctag                           >  2:1288785/1‑146 (MQ=255)
                                                                                                                            gcCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCt                     <  1:444888/150‑1 (MQ=255)
                                                                                                                             ccTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCTGATAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCtt                    <  2:38447/150‑1 (MQ=255)
                                                                                                                              cTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTtc                   <  1:1051428/150‑1 (MQ=255)
                                                                                                                                 aaGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTtcatc                >  2:157568/1‑150 (MQ=255)
                                                                                                                                   ggACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAc                                  <  1:1194637/130‑1 (MQ=255)
                                                                                                                                   ggACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTAc                                  >  2:1194637/1‑130 (MQ=255)
                                                                                                                                         cGAGCCTCGGCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGCGCGTTACCTCATGGGGCGCTTCATCACGGCAcc        >  2:248675/1‑150 (MQ=255)
                                                                                                                                         cGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAg                                        <  2:1145633/118‑1 (MQ=255)
                                                                                                                                         cGAGCATCGGGCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAg                                        >  1:1145633/1‑118 (MQ=255)
                                                                                                                                               tCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCa  >  1:1413736/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCA  >  NC_002947/2730791‑2731083

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