Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,618,848 C→G 7.5% G233G (GGC→GGG)  codA → cytosine deaminase/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,618,8480CG 7.5% 162.4 / 3.3 53G233G (GGC→GGG) codAcytosine deaminase/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base G (2/2);  total (26/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 = 9.26e-01

TACCCGCGACCAGGGGGTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCTGCAAGGGCGCTTC  >  NC_002947/3618701‑3618995
                                                                                                                                                   |                                                                                                                                                   
tACCCGCGACCAGGGGGTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCa                                                                                                                                                    <  2:117211/149‑1 (MQ=255)
            gggggTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCggg                                                                                                                                        >  8:69814/1‑149 (MQ=255)
             ggggTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCgg                                                                                                                                         >  4:147285/1‑147 (MQ=255)
             ggggTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCgg                                                                                                                                         >  6:222433/1‑147 (MQ=255)
             ggggTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCgg                                                                                                                                         <  5:222433/147‑1 (MQ=255)
             ggggTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCgg                                                                                                                                         <  3:147285/147‑1 (MQ=255)
                gTATGGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATCTGCACTGGGAGGTGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCTTGGGTCAc                                                                                                                                    <  7:249554/149‑1 (MQ=255)
                      tcgGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGgccag                                                                                                                              <  5:276996/149‑1 (MQ=255)
                                        gACCTGGCCGAGCGCACCGGCTGCCTGGGGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATgg                                                                                                             <  2:18276/148‑1 (MQ=255)
                                           cTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTc                                                                                                         >  2:34188/1‑149 (MQ=255)
                                                    cgcACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGTGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACaa                                                                                                <  2:78663/149‑1 (MQ=255)
                                                       aCCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGc                                                                                             <  7:69814/149‑1 (MQ=255)
                                                        ccGGCGGCCTGGTGGATGTGCACTGCTACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCGCGTATGACAATGCg                                                                                            <  7:254766/149‑1 (MQ=255)
                                                               ccTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGc                                                                                      >  2:25384/1‑148 (MQ=255)
                                                                 tggtggATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTa                                                                                    >  1:312859/1‑148 (MQ=255)
                                                                 tggtggATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACAACGTGGCCATGGGCTCGTATGACAATGCGTATTGCTa                                                                                    >  6:195195/1‑148 (MQ=255)
                                                                    tggATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAg                                                                                <  3:292004/149‑1 (MQ=255)
                                                                     ggATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGc                                                                               >  1:235827/1‑149 (MQ=255)
                                                                             aCTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGc                                                                       <  8:355988/149‑1 (MQ=255)
                                                                             aCTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGc                                                                       >  7:355988/1‑149 (MQ=255)
                                                                                   acgaGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTcc                                                                         <  4:80755/141‑1 (MQ=255)
                                                                                   acgaGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTcc                                                                         >  3:80755/1‑141 (MQ=255)
                                                                                      agaCCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGGGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAg                                                              <  1:373554/149‑1 (MQ=255)
                                                                                        aCCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCt                                                            >  6:184788/1‑149 (MQ=255)
                                                                                                cccGCAATCACTTTTCCTCGAAGTTTTGGCCTAAGAAGCCCTGGTAAGCGGCATGTGGGCCAGTGTCACGGCCAGCCACACCGTGGCCATGGGCTCTTATTACAATGCGTATTTCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCa                                                    <  7:168006/149‑1 (MQ=255)
                                                                                                cccGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGc                                                                   <  2:223728/134‑1 (MQ=255)
                                                                                                cccGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGc                                                                   >  1:223728/1‑134 (MQ=255)
                                                                                                    cAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCa                                                 >  7:228688/1‑148 (MQ=255)
                                                                                                     aaTCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTa                                                                                          >  1:256430/1‑106 (MQ=255)
                                                                                                     aaTCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTa                                                                                          <  2:256430/106‑1 (MQ=255)
                                                                                                         aCGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtc                                            <  2:312859/148‑1 (MQ=255)
                                                                                                               ccTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtcgtgtcgt                                     >  3:369973/1‑149 (MQ=255)
                                                                                                               ccTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtcgtgtcgt                                     >  7:78743/1‑149 (MQ=255)
                                                                                                                cTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtcgtgtcgt                                     <  5:195195/148‑1 (MQ=255)
                                                                                                                 tCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACCATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGc                                   >  5:100527/1‑149 (MQ=255)
                                                                                                                 tCGAAGTGTTGGCCGAAGAAGCACGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCACGCTGTCGGGCATCAACTTCGTGTCGTGc                                   >  3:173942/1‑149 (MQ=255)
                                                                                                                            gCCGAAGAAGCCCGGGTACGCGGGATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCccgaccga                           >  5:231814/1‑146 (MQ=255)
                                                                                                                            gCCGAAGAAGCCCGGGTACGCGGGATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGagag                        >  5:237245/1‑149 (MQ=255)
                                                                                                                            gCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGCGCATCAACTTCGTGTCGTGCCCGACCgaga                         <  1:34188/148‑1 (MQ=255)
                                                                                                                              cGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGc                                                                                                                  <  7:323215/57‑1 (MQ=255)
                                                                                                                              cGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGc                                                                                                                  >  8:323215/1‑57 (MQ=255)
                                                                                                                              cGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCa                      >  2:44102/1‑149 (MQ=255)
                                                                                                                               gaagaaGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTa                                                                                                      >  1:277880/1‑68 (MQ=255)
                                                                                                                               gaagaaGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTa                                                                                                      <  2:277880/68‑1 (MQ=255)
                                                                                                                                aagaagCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCAtt                    <  2:235827/149‑1 (MQ=255)
                                                                                                                                     gCCCGGGTACGCGGGATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCt               <  6:231814/149‑1 (MQ=255)
                                                                                                                                     gCCCGGGTACGCGGGATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCt               <  6:237245/149‑1 (MQ=255)
                                                                                                                                     gCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCt               >  3:86238/1‑149 (MQ=255)
                                                                                                                                          ggTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCTGCAAg          <  1:44102/149‑1 (MQ=255)
                                                                                                                                           ttACGCGGCATGGGTGCCCGGGTCAAGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCTGCAAgg         <  4:173942/148‑1 (MQ=255)
                                                                                                                                           gTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtcgtgtc                                       >  6:208930/1‑119 (MQ=255)
                                                                                                                                           gTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTtcgtgtc                                       <  5:208930/119‑1 (MQ=255)
                                                                                                                                                  gCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCTGCAAGGGCGCTTc  >  2:22707/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACCCGCGACCAGGGGGTATCGTCGGTGAAGTTTCTGATGGACCTGGCCGAGCGCACCGGCTGCCTGGTGGATGTGCACTGCGACGAGACCGACGACCCGCAATCACGTTTCCTCGAAGTGTTGGCCGAAGAAGCCCGGGTACGCGGCATGGGGGCCCGGGTCACGGCCAGCCACACCGTGGCCATGGGCTCGTATGACAATGCGTATTGCTACAAGCTGTTCCGCCTGCTCAAGCTGTCGGGCATCAACTTCGTGTCGTGCCCGACCGAGAGCATTCACCTGCAAGGGCGCTTC  >  NC_002947/3618701‑3618995

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