Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,385,129 C→T 7.1% E111K (GAG→AAG)  PP_2986 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,385,1290CT 7.1% 158.7 / 2.7 56E111K (GAG→AAG) PP_2986oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (30/22);  new base T (2/2);  total (32/24)
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.84e-01

GTCAGCGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGCCATG  >  NC_002947/3384980‑3385277
                                                                                                                                                     |                                                                                                                                                    
gTCAGCGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACCTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTc                                                                                                                                                      <  1:1223026/150‑1 (MQ=255)
   aGCGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTTGGCCTGCCCAAg                                                                                                                                                   <  2:1013408/150‑1 (MQ=255)
      gAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAg                                                                                                                                                <  1:862220/150‑1 (MQ=255)
              tgctgcCGTAGTAGTCATCGGCGTCGGCGTCCGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGACCCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCcgcg                                                                                                                                        <  1:1040704/150‑1 (MQ=255)
                ctgcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGTAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGgcgtc                                                                                                                                           <  1:44966/145‑3 (MQ=255)
                ctgcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGTAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGcgcc                                                                                                                                           >  2:44966/1‑145 (MQ=255)
                 tgcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCaa                                                                                                                                                    >  1:82362/1‑135 (MQ=255)
                 tgcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCaa                                                                                                                                                    <  2:82362/135‑1 (MQ=255)
                        gtagtCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGacac                                                                                                                              <  2:186917/150‑1 (MQ=255)
                              aTCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTAt                                                                                                                        <  2:528709/150‑1 (MQ=255)
                                  gcgtcggcgtcgGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTAGTATGCCGGCTCGCTAGGCATGTCCAAGCAGGGGCCGCGTACCCGACACGTTTATCTGt                                                                                                                    >  2:345599/1‑150 (MQ=255)
                                    gtcggcgtcggcCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTc                                                                                                                  <  2:433535/150‑1 (MQ=255)
                                                ccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcct                                                                                                                                            >  2:325235/1‑110 (MQ=255)
                                                ccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcct                                                                                                                                            <  1:965736/112‑3 (MQ=255)
                                                ccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcct                                                                                                                                            <  1:325235/112‑3 (MQ=255)
                                                ccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcct                                                                                                                                            >  2:965736/1‑110 (MQ=255)
                                                 cTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAg                                                                                                     <  1:345599/150‑1 (MQ=255)
                                                  tCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCCATGCGTGAGGTAGGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGt                                                                                                    >  2:437601/1‑150 (MQ=255)
                                                      gCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGa                                                                                                                >  2:834839/1‑134 (MQ=255)
                                                      gCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGa                                                                                                                <  1:834839/134‑1 (MQ=255)
                                                      gCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCt                                                                                                >  2:29693/1‑150 (MQ=255)
                                                      gCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCt                                                                                                >  2:753806/1‑150 (MQ=255)
                                                        aGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTgc                                                                                              >  1:346150/1‑150 (MQ=255)
                                                             tcctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGccc                                                                                         >  1:1144723/1‑150 (MQ=255)
                                                              cctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCt                                                                                        >  1:1024163/1‑150 (MQ=255)
                                                              cctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGACCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCAATTGGCTTCCGCCCt                                                                                        >  2:714349/1‑150 (MQ=255)
                                                               ctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCtt                                                                                       >  1:1076426/1‑150 (MQ=255)
                                                                tcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTg                                                                                      >  2:1098208/1‑150 (MQ=255)
                                                                  gATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTa                                                                                    <  2:1076426/150‑1 (MQ=255)
                                                                  gATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTa                                                                                    <  2:1144723/150‑1 (MQ=255)
                                                                          tGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                            <  1:305489/150‑1 (MQ=255)
                                                                          tGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                            >  2:644119/1‑150 (MQ=255)
                                                                           gCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgcg                                                                           >  2:390536/1‑150 (MQ=255)
                                                                                            gCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGa                                                                                                                >  1:4201/1‑96 (MQ=255)
                                                                                            gCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGa                                                                                                                <  2:4201/96‑1 (MQ=255)
                                                                                                   aaTGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAacac                                                   >  1:570486/1‑150 (MQ=255)
                                                                                                    aTGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACAcc                                                  <  1:29693/150‑1 (MQ=255)
                                                                                                       gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTa                                               >  2:690574/1‑150 (MQ=255)
                                                                                                       gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTa                                               <  1:753806/150‑1 (MQ=255)
                                                                                                                         cgtcgtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAt                             >  1:1147032/1‑150 (MQ=255)
                                                                                                                          gtcgtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAtg                            >  1:857296/1‑150 (MQ=255)
                                                                                                                           tcgtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAtgc                           >  2:514403/1‑150 (MQ=255)
                                                                                                                                   cGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGt                                                                                                                    >  1:32131/1‑53 (MQ=255)
                                                                                                                                   cGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGt                                                                                                                    <  2:32131/53‑1 (MQ=255)
                                                                                                                                    ggCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTAt                                              >  2:1071941/1‑122 (MQ=255)
                                                                                                                                    ggCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTAt                                              <  1:1071941/122‑1 (MQ=255)
                                                                                                                                    ggCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGc                  >  1:194648/1‑150 (MQ=255)
                                                                                                                                      ccccGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGt                <  2:194648/150‑1 (MQ=255)
                                                                                                                                      ccccGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGt                <  1:714349/150‑1 (MQ=255)
                                                                                                                                      ccccGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGt                >  2:854166/1‑150 (MQ=255)
                                                                                                                                         cGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGa             <  1:1098208/150‑1 (MQ=255)
                                                                                                                                          gTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGAc            >  1:229950/1‑150 (MQ=255)
                                                                                                                                            aGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCt          >  1:1237645/1‑150 (MQ=255)
                                                                                                                                                cTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGAAGCTGTGCAAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGc      >  2:675998/1‑150 (MQ=255)
                                                                                                                                                    tCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTc               <  1:1257640/137‑1 (MQ=255)
                                                                                                                                                    tCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTc               >  2:1257640/1‑137 (MQ=255)
                                                                                                                                                    tCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGCCATg  >  2:40921/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTCAGCGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGCCATG  >  NC_002947/3384980‑3385277

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