Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,730,937 C→G 10.9% R77R (CGC→CGG)  polB → DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,730,9370CG10.9% 164.7 / 12.0 55R77R (CGC→CGG) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base G (3/3);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTG  >  NC_002947/2730790‑2731081
                                                                                                                                                   |                                                                                                                                                
cGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGccc                                                                                                                                                >  2:605040/1‑150 (MQ=255)
       gggCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCATGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCACGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                         <  1:1570596/150‑1 (MQ=255)
          ccGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGc                                                                                                                                      <  1:24675/150‑1 (MQ=255)
            gcgACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGTCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCct                                                                                                                                    >  2:1102258/1‑150 (MQ=255)
            gcgACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGTCct                                                                                                                                    >  1:471698/1‑150 (MQ=255)
              gACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                      >  1:1159577/1‑144 (MQ=255)
              gACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                      <  2:1159577/146‑3 (MQ=255)
              gACTGTTGCGGCTGGCGGCGCAGGACCGGTTTGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCctct                                                                                                                                  <  1:528/150‑1 (MQ=255)
                  gTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATAc                                                                                                                              <  1:748938/150‑1 (MQ=255)
                    tGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGTCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATAccc                                                                                                                            <  1:1102258/150‑1 (MQ=255)
                             cggcgCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTAAAGCCACTGCGCCTGAAGGACTTCTAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATc                                                                                                                   <  1:446926/150‑1 (MQ=255)
                                 gcAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCa                                                                                                               >  2:531577/1‑150 (MQ=255)
                                           gTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCaa                                                                                                     >  1:146552/1‑150 (MQ=255)
                                           gTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCaa                                                                                                     >  1:166495/1‑150 (MQ=255)
                                             ggCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACt                                                                                                   <  2:1480374/150‑1 (MQ=255)
                                                tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTggg                                                                                                >  2:449420/1‑149 (MQ=255)
                                                tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                >  2:1399744/1‑150 (MQ=255)
                                                tttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGa                                                                                                >  1:755902/1‑150 (MQ=255)
                                                 ttCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGag                                                                                               <  1:449420/150‑1 (MQ=255)
                                                         tCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGTTCGCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCt                                                                                       >  1:514509/1‑150 (MQ=255)
                                                                  ggCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGc                                                                              <  2:713204/150‑1 (MQ=255)
                                                                   gCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGag                                                                                               <  1:1391655/132‑1 (MQ=255)
                                                                   gCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGag                                                                                               >  2:1391655/1‑132 (MQ=255)
                                                                   gCCGAGCACGCCCGGGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCCCTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGcc                                                                             >  1:292840/1‑150 (MQ=255)
                                                                    ccGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                            <  2:755902/150‑1 (MQ=255)
                                                                     cGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                         >  1:657419/1‑120 (MQ=255)
                                                                     cGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGAt                                                                                                         <  2:657419/120‑1 (MQ=255)
                                                                              ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                  >  2:1560206/1‑150 (MQ=255)
                                                                                        ctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGcc                                                        <  1:605040/150‑1 (MQ=255)
                                                                                        ctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                      >  2:142786/1‑70 (MQ=37)
                                                                                        ctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTCCGAGCATCGGGCGATGCTCgaa                                                                                                                                      <  1:142786/72‑3 (MQ=37)
                                                                                         tgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTACCGAAGTGTTCGAGGCCg                                                       >  2:913163/1‑150 (MQ=255)
                                                                                              aaCGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATc                                                  >  1:1604157/1‑150 (MQ=255)
                                                                                                cGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCg                                                <  2:146552/150‑1 (MQ=255)
                                                                                                      cGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGGGCGATGCTCgaa                                                                                                                                      <  1:817477/58‑3 (MQ=37)
                                                                                                      cGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGGGCGATGCTCgaa                                                                                                                                      >  2:817477/1‑56 (MQ=37)
                                                                                                            gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgtgt                                                                >  1:678765/1‑122 (MQ=255)
                                                                                                            gAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgtgt                                                                <  2:678765/122‑1 (MQ=255)
                                                                                                                 tCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGNCCGATATCCGGCCGCCGGAGCGTTAcc                               <  1:1560206/150‑1 (MQ=255)
                                                                                                                   aaGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                  <  2:718744/113‑1 (MQ=255)
                                                                                                                   aaGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                  >  1:718744/1‑113 (MQ=255)
                                                                                                                          tGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGa                                                                     <  2:126374/103‑1 (MQ=255)
                                                                                                                          tGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGa                                                                     >  1:126374/1‑103 (MQ=255)
                                                                                                                             gcCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCt                   >  1:424169/1‑150 (MQ=255)
                                                                                                                                   aGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTtcatca             <  1:1076479/150‑1 (MQ=255)
                                                                                                                                        ttCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCa        >  1:1113052/1‑150 (MQ=255)
                                                                                                                                           gAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGgcc                                             <  1:684734/110‑1 (MQ=255)
                                                                                                                                           gAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGgcc                                             >  2:684734/1‑110 (MQ=255)
                                                                                                                                           gAGCATCGCCCGATGCGCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCg     <  2:166495/150‑1 (MQ=255)
                                                                                                                                             gCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCGCGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCt                              <  2:1216620/123‑1 (MQ=255)
                                                                                                                                             gCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGgccgcc                                          <  2:305192/111‑1 (MQ=255)
                                                                                                                                             gCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGgccgcc                                          >  1:305192/1‑111 (MQ=255)
                                                                                                                                             gCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCt                              >  1:1216620/1‑123 (MQ=255)
                                                                                                                                              cATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTg  <  2:424169/150‑1 (MQ=255)
                                                                                                                                              cATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTg  >  2:1468548/1‑150 (MQ=255)
                                                                                                                                              cATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCCCCGGAGCGTTACCCCATGGAGCGCTTCATCACGGCACCGGtc  >  2:291438/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTG  >  NC_002947/2730790‑2731081

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