Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,001,813 A→G 7.2% Q585R (CAG→CGG)  PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,001,8130AG 7.2% 158.5 / 2.3 55Q585R (CAG→CGG) PP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

ATTGAATCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTG  >  NC_002947/2001673‑2001961
                                                                                                                                             |                                                                                                                                                    
aTTGAATCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGtt                                                                                                                                              >  6:285258/1‑149 (MQ=255)
aTTGAATCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGtt                                                                                                                                              >  7:306174/1‑149 (MQ=255)
    aaTCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGt                                                                                                                                          >  3:132268/1‑149 (MQ=255)
    aaTCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGt                                                                                                                                          >  2:284567/1‑149 (MQ=255)
        aGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCgtt                                                                                                                                      <  5:285258/149‑1 (MQ=255)
              gCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGATTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTgg                                                                                                                                           <  1:106755/138‑1 (MQ=255)
              gCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTgg                                                                                                                                           >  2:106755/1‑138 (MQ=255)
              gCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTg                                                                                                                                            >  5:67739/1‑137 (MQ=255)
              gCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTg                                                                                                                                            <  6:67739/137‑1 (MQ=255)
                        cAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGAt                                                                                                                      <  4:132268/149‑1 (MQ=255)
                              gCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAg                                                                                                                >  8:261040/1‑149 (MQ=255)
                               cacaAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAgt                                                                                                               <  6:104230/149‑1 (MQ=255)
                                acaAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAgtg                                                                                                              >  2:214050/1‑149 (MQ=255)
                                  aagaagACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAgtgt                                                                                                             >  3:316208/1‑148 (MQ=255)
                                  aagaagACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCG‑TTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCa                                                                                                           <  4:220538/149‑1 (MQ=255)
                                   agaagaCCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCa                                                                                                           <  1:284567/149‑1 (MQ=255)
                                       gaCCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCAttt                                                                                                        >  4:19808/1‑148 (MQ=255)
                                                  tcttctATGGCAGCGGTACTCAAGCACACAACAACGACATTATCGAACAGGCGCTGCCTGCGATTGA‑GCGAGAGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTTCCATTGTCATTTACAGCAAAATAc                                                                                            >  5:47876/1‑149 (MQ=255)
                                                    ttctATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACAGAAAAATACGc                                                                                          >  6:110600/1‑149 (MQ=255)
                                                      ctATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACg                                                                                           <  4:310827/146‑1 (MQ=255)
                                                               cGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTa                                                                               <  8:306174/149‑1 (MQ=255)
                                                                  tACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTAccc                                                                            <  4:316208/149‑1 (MQ=255)
                                                                    cTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAgtg                                                                                                              <  6:294230/113‑1 (MQ=255)
                                                                    cTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAgtg                                                                                                              >  5:294230/1‑113 (MQ=255)
                                                                     tCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGt                                                                         >  6:52694/1‑149 (MQ=255)
                                                                       aaGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGt                                                                                                                                <  8:135208/92‑1 (MQ=255)
                                                                       aaGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGt                                                                                                                                >  7:135208/1‑92 (MQ=255)
                                                                       aaGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTg                                                                       >  6:241105/1‑149 (MQ=255)
                                                                        aGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGc                                                                      <  7:261040/149‑1 (MQ=255)
                                                                              caactCCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTAccc                                                                            <  2:357630/137‑1 (MQ=255)
                                                                               aacaCCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCAt                                                               >  1:110166/1‑149 (MQ=255)
                                                                                  acCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCGGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGa                                                            >  2:117210/1‑149 (MQ=255)
                                                                                  acCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGa                                                            >  1:15074/1‑149 (MQ=255)
                                                                                    cGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACa                                                          <  5:52694/149‑1 (MQ=255)
                                                                                    cGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACa                                                          <  2:110166/149‑1 (MQ=255)
                                                                                      aCTTTATCGAACAGTCGCTGCCTGCGTTTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGACGTTTTTGGTTACCTGCGATTGCCAGTGTAATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAgt                                                        <  2:50368/149‑1 (MQ=255)
                                                                                        tttATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAgtgt                                                      >  4:189342/1‑149 (MQ=255)
                                                                                           aTCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGc                                                    >  8:117115/1‑148 (MQ=255)
                                                                                            tCGAACAGGCGCTTCCTGCGATTGA‑GCGTGTGCTTAACGATACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCaa                                                  <  3:19808/149‑1 (MQ=255)
                                                                                            tCGAACAGGCGCTGCCTGCGATTGAGGCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGc                                                    >  1:188057/1‑148 (MQ=255)
                                                                                              gAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCAT‑GACAGTGTGCAAGGt                                               >  8:61400/1‑149 (MQ=255)
                                                                                                  aGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTAt                                            <  5:235577/149‑1 (MQ=255)
                                                                                                    gcgcTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCGGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTg                                          <  1:117210/149‑1 (MQ=255)
                                                                                                    gcgcTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGAtt                                                                                                                     >  8:43993/1‑74 (MQ=255)
                                                                                                    gcgcTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGAtt                                                                                                                     <  7:43993/74‑1 (MQ=255)
                                                                                                           ctgcGATTGA‑GCGTGTGCTTAACGAGACCCCTCGGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCgct                                    >  7:329841/1‑148 (MQ=255)
                                                                                                               gATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCAT‑GACAGTGTGCAAGGTTATTGGTCGCTGCTTGa                              <  7:61400/149‑1 (MQ=255)
                                                                                                                aTTGA‑GCGTGTGGTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTTTTGGTTACCTGAGATTGCCAGTTTCATTTAGCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGa                              <  6:47876/149‑1 (MQ=255)
                                                                                                                           tgCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGat                    <  7:117115/149‑1 (MQ=255)
                                                                                                                             cTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGata                   <  1:214050/148‑1 (MQ=255)
                                                                                                                                 aCGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCaa               <  5:241105/148‑1 (MQ=255)
                                                                                                                                   gagaCCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAatat            >  8:332517/1‑149 (MQ=255)
                                                                                                                                       ccccTCGGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAACATTGCg        <  8:329841/149‑1 (MQ=255)
                                                                                                                                       ccccTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCg        <  7:332517/149‑1 (MQ=255)
                                                                                                                                       ccccTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCg        <  3:189342/149‑1 (MQ=255)
                                                                                                                                             aGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGgtg     >  6:277115/1‑146 (MQ=255)
                                                                                                                                             aGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGgtg     <  5:277115/146‑1 (MQ=255)
                                                                                                                                             aGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTg  <  2:15074/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
ATTGAATCAGATTCGCACGAACGACAAATCGCACAAGAAGACCATCGATATCTTCTATGGCAGCGGTACTCAAGCACACAACACCGACTTTATCGAACAGGCGCTGCCTGCGATTGA‑GCGTGTGCTTAACGAGACCCCTCAGGCTCGTTTGGTCGTTGTTGGTTACCTGCGATTGCCAGTGTCATTTACCGCAAAATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTG  >  NC_002947/2001673‑2001961

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