Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,730,937 C→G 13.1% 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,9370CG13.1% 123.9 / 11.8 46R77R (CGC→CGGpolBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base G (3/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 = 5.90e-01

ACCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCA  >  NC_002947/2730788‑2731083
                                                                                                                                                     |                                                                                                                                                  
aCCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGc                                                                                                                                                    >  2:444888/1‑150 (MQ=255)
   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:593602/150‑1 (MQ=255)
                                                                      ccGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCg                                                                              <  2:74440/150‑1 (MQ=255)
                                                                        gAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGt                                                                            <  1:99754/150‑1 (MQ=255)
                                                                                ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                           >  2:733949/1‑79 (MQ=255)
                                                                                ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTc                                                                                                                                           <  1:733949/79‑1 (MQ=255)
                                                                                ccGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACgt                                                                    >  1:15774/1‑150 (MQ=255)
                                                                                       ctgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAg                                                             <  1:280982/150‑1 (MQ=255)
                                                                                        tgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                              <  1:1248602/148‑1 (MQ=255)
                                                                                        tgctgGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGa                                                              >  2:1248602/1‑148 (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)
                                                                                                                                                     |                                                                                                                                                  
ACCGATGCCGGGCCGCGACTGTTGCGGCTGGCGGCGCAGGAGTCGGTGGCTTTCATTCCTCAGGCGCAGGCCGAGCACGCCCGTGTGCTGCTGGCCAACGAAACCGGGACGAAGCTCAAGCCACTGCGCCTGAAGGACTTCGAGCATCGCCCGATGCTCGGCCTCTATACCCGCCAGCATCGCCAATTGATGCAACTGGAGCAGCGTCTGCGCACGGCCGGTATCGACGTGTTCGAGGCCGATATCCGGCCGCCGGAGCGTTACCTCATGGAGCGCTTCATCACGGCACCGGTGCA  >  NC_002947/2730788‑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: