Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,811,541 G→T 8.2% Q181H (CAG→CAT)  dsbG → protein disulfide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,811,5410GT 8.2% 151.0 / 4.0 49Q181H (CAG→CAT) dsbGprotein disulfide isomerase
Reads supporting (aligned to +/- strand):  ref base G (17/28);  new base T (2/2);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CAACATGTTCTGGGAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAG  >  NC_002947/4811397‑4811685
                                                                                                                                                |                                                                                                                                                
cAACATGTTCTGGGAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTgcagcag                                                                                                                                               >  3:64789/1‑148 (MQ=255)
  aCATGTTCTGGGAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCAt                                                                                                                                            >  1:6637/1‑149 (MQ=255)
  aCATGTTCTGGGAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCAt                                                                                                                                            >  7:195367/1‑149 (MQ=255)
            ggAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCg                                                                                                                                  <  4:64789/149‑1 (MQ=255)
                   gcgcGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGCCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGc                                                                                                                           >  8:172457/1‑149 (MQ=255)
                     gcgGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCa                                                                                                                         <  3:131644/149‑1 (MQ=255)
                        gCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCa                                                                                                                      <  5:56948/149‑1 (MQ=255)
                                tGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAg                                                                                                              <  6:152651/149‑1 (MQ=255)
                                  gAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAgg                                                                                                             <  3:253997/148‑1 (MQ=255)
                                    tacGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCAc                                                                                                          <  2:83378/147‑1 (MQ=255)
                                       gggCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTc                                                                                                        <  4:162718/148‑1 (MQ=255)
                                          cAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTc                                                                                                        <  5:102358/145‑1 (MQ=255)
                                          cAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTc                                                                                                        >  6:102358/1‑145 (MQ=255)
                                          cAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGc                                                                                                    <  3:244617/149‑1 (MQ=255)
                                              gTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGt                                                                                                <  5:40161/149‑1 (MQ=255)
                                                gCAGTTGCGCCATATCAT‑GTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGc                                                                                              >  2:266647/1‑148 (MQ=255)
                                                       cgcCATATCATGGTCTGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAgg                                                                                        <  5:139507/148‑1 (MQ=255)
                                                        gcCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCt                                                                                      <  2:9265/149‑1 (MQ=255)
                                                         ccATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCtg                                                                                     >  5:175909/1‑149 (MQ=255)
                                                           atatCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCtgtg                                                                                   <  8:121445/149‑1 (MQ=255)
                                                                       gtatcatcCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTg                                                                       <  5:239709/147‑1 (MQ=255)
                                                                        catcatCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTgg                                                                      >  3:139367/1‑149 (MQ=255)
                                                                         atcatcCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTggc                                                                     <  4:139367/149‑1 (MQ=255)
                                                                            atcCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTggcggc                                                                  <  8:195367/149‑1 (MQ=255)
                                                                            atcCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTggcgg                                                                   <  6:175909/148‑1 (MQ=255)
                                                                                gcgAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGaa                                                               >  4:108147/1‑148 (MQ=255)
                                                                                gcgAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGaa                                                               >  7:60846/1‑148 (MQ=255)
                                                                                   aGGAAAGCCCGGGCATGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATg                                                           <  3:174530/149‑1 (MQ=255)
                                                                                          cccGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGt                                                                                                >  4:171441/1‑105 (MQ=255)
                                                                                          cccGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGt                                                                                                <  3:171441/105‑1 (MQ=255)
                                                                                          cccGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGa                                                    >  1:213110/1‑149 (MQ=255)
                                                                                                cAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCAGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAag                                              <  3:210443/149‑1 (MQ=255)
                                                                                                     cggcggCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCATCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGAt                                                   >  4:69298/1‑139 (MQ=255)
                                                                                                     cggcggCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCATCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGAt                                                   >  4:73493/1‑139 (MQ=255)
                                                                                                     cggcggCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCATCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGAt                                                   <  3:73493/139‑1 (MQ=255)
                                                                                                     cggcggCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCATCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGAt                                                   <  3:69298/139‑1 (MQ=255)
                                                                                                            gctgctGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGTCGTTGATGGAAGAGATGGGGTTGc                                  <  1:266647/149‑1 (MQ=255)
                                                                                                              tgctgGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAg                                >  7:138584/1‑149 (MQ=255)
                                                                                                                ctgGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGc                              <  2:6637/149‑1 (MQ=255)
                                                                                                                   gcggcgAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGa                            <  7:172457/148‑1 (MQ=255)
                                                                                                                   gcggcgAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGAc                           >  3:105792/1‑149 (MQ=255)
                                                                                                                     ggcgAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTc                         >  3:119661/1‑149 (MQ=255)
                                                                                                                        gAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGTCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAgg                               <  3:179822/140‑1 (MQ=255)
                                                                                                                        gAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAgg                               >  4:179822/1‑140 (MQ=255)
                                                                                                                                     aaGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACTCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCa         <  5:311651/149‑1 (MQ=255)
                                                                                                                                       ggCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAg        >  7:75771/1‑148 (MQ=255)
                                                                                                                                        gCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCAGGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGa      <  8:72534/149‑1 (MQ=255)
                                                                                                                                           cTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCt             >  4:68119/1‑139 (MQ=255)
                                                                                                                                           cTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCt             <  3:68119/139‑1 (MQ=255)
                                                                                                                                           cTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGa   <  2:213110/149‑1 (MQ=255)
                                                                                                                                             gcagcagcaTGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAg  >  5:271561/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAACATGTTCTGGGAACAGGCGCGGCCCTGGGTGGAGTCGGGCAAGGTGCAGTTGCGCCATATCATGGTCGGCATCATCCGCGAGGACAGCCCGGGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAG  >  NC_002947/4811397‑4811685

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