Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,173,836 A→C 9.5% I64S (ATC→AGC) ‡ PP_1926 ← phosphotyrosine‑protein phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,173,8360AC 9.5% 147.3 / 5.7 52I64S (ATC→AGC) ‡PP_1926phosphotyrosine‑protein phosphatase
Reads supporting (aligned to +/- strand):  ref base A (26/21);  new base C (2/3);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCGTAAAAATCA  >  NC_002947/2173694‑2173977
                                                                                                                                              |                                                                                                                                             
gATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGcc                                                                                                                                           >  2:159669/1‑147 (MQ=255)
    gCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCATCGCGGCCAAGca                                                                                                                                      <  3:142695/148‑2 (MQ=255)
      cAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTggg                                                                                                                                   >  2:35287/1‑149 (MQ=255)
       aGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGTTCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGc                                                                                                                                  <  2:184467/149‑1 (MQ=255)
       aGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCACCCCATGCTGCTTGGCCGCGATGCCAACATggg                                                                                                                                   >  8:214926/1‑148 (MQ=255)
       aGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTggg                                                                                                                                   >  4:284429/1‑148 (MQ=255)
       aGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGc                                                                                                                                  <  6:140765/149‑1 (MQ=255)
         ggCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCaa                                                                                                                                <  1:159669/149‑1 (MQ=255)
            cTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCATCGCGGCCAAGcagc                                                                                                                                    <  6:254660/142‑4 (MQ=255)
            cTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCATCGCGGCCAAGcagc                                                                                                                                    <  6:250521/142‑4 (MQ=255)
            cTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCATCGCGGCCAAGcagc                                                                                                                                    >  5:254660/1‑139 (MQ=255)
            cTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCATCGCGGCCAAGcagc                                                                                                                                    >  5:250521/1‑139 (MQ=255)
                    aaGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTcc                                                                                                                     <  3:284429/149‑1 (MQ=255)
                             aaGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGctct                                                                                                            <  7:214926/149‑1 (MQ=255)
                                 gTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGt                                                                                                        >  5:240974/1‑149 (MQ=255)
                                  gccAACTGTTGGCGAATCTCTACCAGCCGAGACTCGAAGTTTTTATCAGGAGCCTGATCATCCTCGATCGGCAGCGTGAACCCCTCCAACCCATTCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTc                                                                                                       <  4:278987/148‑1 (MQ=255)
                                    cAACTGGTGGCGAATCTCTCCCAGCCTAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGGAACTCCGTTTTCCTGAAGCTCTGCGTcc                                                                                                      <  2:90312/148‑1 (MQ=255)
                                          gtggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACAt                                                                                                 >  8:88103/1‑147 (MQ=255)
                                                aaTCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTaa                                                                                         <  4:144371/149‑1 (MQ=255)
                                                 aTCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTaa                                                                                         >  4:234531/1‑148 (MQ=255)
                                                 aTCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAAc                                                                                        <  7:61097/149‑1 (MQ=255)
                                                  tctctcCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAAc                                                                                        >  6:185729/1‑148 (MQ=255)
                                                      tcCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAgcg                                                                                   >  7:54058/1‑149 (MQ=255)
                                                        ccAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAgcgc                                                                                  >  8:40902/1‑148 (MQ=255)
                                                        ccAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCt                                                                                 <  3:234531/149‑1 (MQ=255)
                                                            ccGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTAc                                                                               <  8:155184/147‑1 (MQ=255)
                                                            ccGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTAc                                                                               >  7:155184/1‑147 (MQ=255)
                                                                   tCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGt                                                                                                                          >  6:38594/1‑97 (MQ=255)
                                                                   tCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGt                                                                                                                          <  5:38594/97‑1 (MQ=255)
                                                                        gTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTgat                                                                 >  1:210425/1‑147 (MQ=255)
                                                                        aggatGATCAGGCTCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCt                                                                 >  8:237321/5‑149 (MQ=255)
                                                                         ttttGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACg                                                                              >  2:216696/1‑135 (MQ=255)
                                                                         ttttGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCCTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACg                                                                              <  1:216696/135‑1 (MQ=255)
                                                                           ttGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCtt                                                                >  5:94106/1‑147 (MQ=255)
                                                                              aTCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGc                                                            <  7:88103/148‑1 (MQ=255)
                                                                              aTCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGc                                                            <  6:240974/148‑1 (MQ=255)
                                                                               tCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGc                                                            <  4:211940/147‑1 (MQ=255)
                                                                               tCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGc                                                            >  3:211940/1‑147 (MQ=255)
                                                                                               ccTCGATCGGCACCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACACCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGc                                           >  4:122849/1‑148 (MQ=255)
                                                                                                cTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTAc                                                                               <  3:16416/111‑1 (MQ=255)
                                                                                                cTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTAc                                                                               >  4:16416/1‑111 (MQ=255)
                                                                                                 tCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGc                                           >  5:277704/1‑146 (MQ=255)
                                                                                                  cGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCgaga                                       >  7:274821/1‑149 (MQ=255)
                                                                                                             tGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTg                            >  3:294373/1‑149 (MQ=255)
                                                                                                                        cAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGa                  >  7:106467/1‑148 (MQ=255)
                                                                                                                        cAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGa                  >  6:137317/1‑148 (MQ=255)
                                                                                                                            ccATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAgg             <  8:54058/149‑1 (MQ=255)
                                                                                                                             cATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGc            >  5:1279/1‑149 (MQ=255)
                                                                                                                              aTGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCg           >  2:277155/1‑149 (MQ=255)
                                                                                                                               tgctgcTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCGt          <  8:106467/149‑1 (MQ=255)
                                                                                                                                        gCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCGTAAAAATCa  <  1:35287/148‑1 (MQ=255)
                                                                                                                                        gCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCGTAAAAATCa  <  5:185729/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCCAGGCTTTCAGCGAGAGCGCTGGTCTGGGTGCCCGGACAAGGCGTAAAAATCA  >  NC_002947/2173694‑2173977

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