Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,500,117 T→C 7.8% R315R (AGA→AGG)  PP_3989 ← DNA‑cytosine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,500,1170TC 7.8% 141.5 / 2.5 51R315R (AGA→AGG) PP_3989DNA‑cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (2/2);  total (26/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.70e-01

TCGGTAGCGTCAACTGCATCTAATGAGACGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGCGACCA  >  NC_002947/4499969‑4500264
                                                                                                                                                    |                                                                                                                                                   
tcgGTAGCGTCAACTGCATCTAATGAGACGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAAt                                                                                                                                                     <  7:21386/149‑1 (MQ=255)
                   cTAATGAGACGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATc                                                                                                                                  >  8:16340/1‑149 (MQ=255)
                      aTGAGACGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCt                                                                                                                               <  6:10037/149‑1 (MQ=255)
                       tGAGACGTCCAAGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGACGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTg                                                                                                                              <  4:213278/149‑1 (MQ=255)
                       tGAGACCTCCAAGGAGCTCCCCCTTACAAAAGGAACTCCTTAAAACCTATCGCCAGTGCTGCTGCCACACGTGGAGGCACAGCGTTAGAAATTTGCTAAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTg                                                                                                                              <  8:303639/149‑1 (MQ=255)
                            cGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCt                                                                                                                         <  5:85536/149‑1 (MQ=255)
                                    tATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGcc                                                                                                                 >  6:34305/1‑149 (MQ=255)
                                         ccccTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGATGTCTGGGAATCCCTGCAGCTTCATTGCCTCGta                                                                                                            >  7:295811/1‑149 (MQ=255)
                                          cccTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGta                                                                                                            <  7:16340/148‑1 (MQ=255)
                                          cccTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGt                                                                                                             <  5:92292/147‑1 (MQ=255)
                                            cTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCAtt                                                                                                                    >  5:112342/1‑138 (MQ=255)
                                            cTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCAtt                                                                                                                    <  6:112342/138‑1 (MQ=255)
                                                aCAAAAGGAACTCTTTAACAAATTTGGCCAGTGCTGCTGCCACAGGTGGAGGAACAGAGTTATAAATTTGCTTAACTTCGGAGCTCAAAGAAACCATCAATCTATAGCTGTATGTGCATCCCTGCAGCTTCATTTCCTCGTATATACt                                                                                                      >  6:332663/1‑148 (MQ=255)
                                                  aaaaGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAg                                                                                                   >  2:52486/1‑149 (MQ=255)
                                                   aaaGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGCCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGc                                                                                                  >  3:156885/1‑149 (MQ=255)
                                                       gAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTc                                                                                              >  2:78770/1‑149 (MQ=255)
                                                        aaCTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCt                                                                                             <  5:34305/149‑1 (MQ=255)
                                                             tttAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTc                                                                                        <  1:78770/149‑1 (MQ=255)
                                                                aaaaCATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTg                                                                                                                              <  6:263045/108‑1 (MQ=255)
                                                                aaaaCATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTg                                                                                                                              >  5:263045/1‑108 (MQ=255)
                                                                  aaCATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCa                                                                                   <  4:351519/149‑1 (MQ=255)
                                                                     aTTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCa                                                                                   >  3:260675/1‑146 (MQ=255)
                                                                     aTTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCa                                                                                   <  4:260675/146‑1 (MQ=255)
                                                                          gCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCa                                                                                                                      >  2:288817/1‑106 (MQ=255)
                                                                          gCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCa                                                                                                                      <  1:288817/106‑1 (MQ=255)
                                                                           ccAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACa                                                                          >  7:79164/1‑149 (MQ=255)
                                                                                 ctgctgcCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAAtc                                                                    <  8:295811/149‑1 (MQ=255)
                                                                                  tgctgcCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTAAATACTAAGCCTTCTATGTCCATCAGGATGAACATGAAtct                                                                   >  1:184915/1‑149 (MQ=255)
                                                                                                gAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAg                                                     <  8:4391/149‑1 (MQ=255)
                                                                                                     acaGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCAc                                                >  2:175494/1‑149 (MQ=255)
                                                                                                          gTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTcc                                           <  1:155360/149‑1 (MQ=255)
                                                                                                           ttAGAAATTTGATTAACTTGGGCGATCAAATTCCCCATCAATCTATAGATGTCTGGGATTCCATGCAGCTTCATTGCCTCGTATATAATAAGCCTACTATGTCCATCAGGATCAACATGAATCTCCATGATACCATAGCCCACCgaacg                                          >  2:229215/1‑145 (MQ=255)
                                                                                                            tAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCgg                                         >  2:239764/1‑149 (MQ=255)
                                                                                                               aaaTTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGc                                        >  5:205782/1‑147 (MQ=255)
                                                                                                               aaaTTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAACCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGc                                        >  1:91595/1‑147 (MQ=255)
                                                                                                               aaaTTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAACCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGc                                        >  1:86993/1‑147 (MQ=255)
                                                                                                                 aTTTGCTTAACTTGGGCGCTCTAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTGGTATATACTAAGCCTTCTATGTCCATGAGGATGAACATGAATCTCCCTGTTACCATAGGCAAGCGTCCGGCTAgg                                    >  4:44807/1‑149 (MQ=255)
                                                                                                                       ttAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGtttttt                              >  5:225927/1‑149 (MQ=255)
                                                                                                                       ttAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGttttt                               >  7:333026/1‑148 (MQ=255)
                                                                                                                        tatcTTGGGCGCTCAAAGTCACCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGAATATACTAAGCCTTCTATGTCCATCAGTATGAACATGAATCTCCCTTTTACCATAGGCCACCGTCCGGCTAGGTTTTTTc                             <  1:229215/146‑1 (MQ=255)
                                                                                                                         aaCTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTc                             <  1:239764/148‑1 (MQ=255)
                                                                                                                                gcgcTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTaa                     <  6:225927/149‑1 (MQ=255)
                                                                                                                                    tCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCATTCTAACtt                  >  3:303488/1‑148 (MQ=255)
                                                                                                                                      aaaGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTaa                <  1:52486/148‑1 (MQ=255)
                                                                                                                                          tCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAAc           <  8:79164/149‑1 (MQ=255)
                                                                                                                                           ccccATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATCAACATGAAAATACCTGTTACCATAGGACACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAAc           >  3:250510/1‑148 (MQ=255)
                                                                                                                                             ccATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAg        >  2:84420/1‑149 (MQ=255)
                                                                                                                                              cATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGc       <  6:205782/149‑1 (MQ=255)
                                                                                                                                              cATCAACCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGc       <  2:86993/149‑1 (MQ=255)
                                                                                                                                              cATCAACCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGc       <  2:91595/149‑1 (MQ=255)
                                                                                                                                                   aTCTATAGCTGTCTGGGAATCCGTGAAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGCGacca  <  3:44807/149‑1 (MQ=255)
                                                                                                                                                   aTCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGCGacca  >  2:94612/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGTAGCGTCAACTGCATCTAATGAGACGTCCATGTATCTCCCCTTTACAAAAGGAACTCTTTAAAACATTTGGCCAGTGCTGCTGCCACAGGTGGAGGCACAGCGTTAGAAATTTGCTTAACTTGGGCGCTCAAAGTCCCCATCAATCTATAGCTGTCTGGGAATCCCTGCAGCTTCATTGCCTCGTATATACTAAGCCTTCTATGTCCATCAGGATGAACATGAATCTCCCTGTTACCATAGGCCACCGTCCGGCTAGGTTTTTTCCACTCTAACTTAATGAACGAGCGACCA  >  NC_002947/4499969‑4500264

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