Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,303 A→T 12.5% L170Q (CTG→CAG)  PP_4218 ← lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,768,3030AT12.5% 143.2 / 11.1 56L170Q (CTG→CAG) PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base T (2/5);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCACTGCGGCGC  >  NC_002947/4768161‑4768436
                                                                                                                                              |                                                                                                                                     
cgcTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTAGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGc                                                                                                                                    <  2:60406/146‑1 (MQ=255)
   ttCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGATaccac                                                                                                                              >  3:11914/1‑149 (MQ=255)
     cAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTaccac                                                                                                                              <  8:124094/147‑1 (MQ=255)
      aaGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCTGGCTaccattct                                                                                                                           <  6:150296/149‑5 (MQ=255)
      aaGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGa                                                                                                                           <  2:230976/149‑1 (MQ=255)
       aGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCTGGCTaccattctg                                                                                                                          <  8:293108/149‑6 (MQ=255)
       aGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCTGGCTaccattctg                                                                                                                          <  6:47486/149‑6 (MQ=255)
        gCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGaca                                                                                                                         >  8:124900/1‑149 (MQ=255)
                cggCTTCCATAGCGCTGCACCGAAGGCCGCCTGCAAGTGGCATCGGTCTCCGGGTAGATCATTACCTGCAGTCGCGGTGCCGGCAGCTCGCGATGCGCAGCCAACTTGTTGGCCAGAATGGTAGCCATGCTACCACCGACACTGTCGcc                                                                                                                 <  4:11914/147‑1 (MQ=255)
                ccGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCGACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTACCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGcc                                                                                                                 <  1:270266/149‑1 (MQ=255)
                 cGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCTGGCTaccattct                                                                                                                           <  8:252835/138‑5 (MQ=255)
                 cGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCTGGCTaccattct                                                                                                                           >  7:252835/1‑134 (MQ=255)
                         tAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGc                                                                                                        <  3:6826/149‑1 (MQ=255)
                          aGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGc                                                                                                                                    <  4:130825/120‑1 (MQ=255)
                          aGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGc                                                                                                                                    >  3:130825/1‑120 (MQ=255)
                           gcgcTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGc                                                                                                        >  5:187134/1‑147 (MQ=255)
                            cgcTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAg                                                                                                     >  6:26483/1‑149 (MQ=255)
                               tGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCg                                                                                                  >  4:179476/1‑149 (MQ=255)
                                   ccGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTTGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCccac                                                                                                              <  6:144088/133‑1 (MQ=255)
                                   ccGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCaccacc                                                                                                          >  3:139764/1‑137 (MQ=255)
                                   ccGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCaccacc                                                                                                          <  4:139764/137‑1 (MQ=255)
                                   ccGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTgc                                                                                              <  6:187134/149‑1 (MQ=255)
                                     gAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTgcgc                                                                                            <  5:26483/149‑1 (MQ=255)
                                      aaGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCt                                                                                                <  2:141447/144‑1 (MQ=255)
                                      aaGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCt                                                                                                >  1:141447/1‑144 (MQ=255)
                                       aGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCAt                                                                                          <  2:170750/149‑1 (MQ=255)
                                                 cAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAg                                                                                <  5:242902/149‑1 (MQ=255)
                                                   aCTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAgcg                                                                              <  3:179476/149‑1 (MQ=255)
                                                       gCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGcc                                                                            <  1:279124/147‑1 (MQ=255)
                                                         aTCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGGGCATCGATGCGCAgcga                                                                             <  4:156389/144‑2 (MQ=255)
                                                         aTCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGGGCATCGATGCGCAgcga                                                                             >  3:156389/1‑143 (MQ=255)
                                                                  cGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGa                                                               >  3:50910/1‑149 (MQ=255)
                                                                         aTCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACGACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGcccacc                                                          >  7:155249/1‑145 (MQ=255)
                                                                            aTTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAgcg                                                     <  8:74797/149‑1 (MQ=255)
                                                                             ttACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGaccaac                                                          >  5:258397/1‑143 (MQ=255)
                                                                             ttACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGaccaac                                                          <  6:258397/143‑1 (MQ=255)
                                                                                 cTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAgg                                                >  1:28532/1‑149 (MQ=255)
                                                                                   gCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCa                                              >  2:192219/1‑149 (MQ=255)
                                                                                               ccGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCg                                  >  2:56893/1‑149 (MQ=255)
                                                                                                    aGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCa                                          <  4:232579/136‑1 (MQ=255)
                                                                                                    aGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCa                                          >  3:232579/1‑136 (MQ=255)
                                                                                                     gCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTg                            <  4:50910/149‑1 (MQ=255)
                                                                                                            cTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCg                     <  4:12887/149‑1 (MQ=255)
                                                                                                              gcgcAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCgg                    >  2:281442/1‑148 (MQ=255)
                                                                                                               cgcAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCggg                   <  2:28532/148‑1 (MQ=255)
                                                                                                                 cAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGaaa                <  1:192219/149‑1 (MQ=255)
                                                                                                                      aaCTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCa           >  6:129181/1‑149 (MQ=255)
                                                                                                                         tGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCa                               >  2:155417/1‑126 (MQ=255)
                                                                                                                         tGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCa                               <  1:155417/126‑1 (MQ=255)
                                                                                                                           gTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCTGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCACTgcg      >  4:267683/1‑149 (MQ=255)
                                                                                                                              ggCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAgg                                        <  5:67119/112‑1 (MQ=255)
                                                                                                                              ggCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAgg                                        >  6:67119/1‑112 (MQ=255)
                                                                                                                              ggCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCACTgcggc    <  7:124900/148‑1 (MQ=255)
                                                                                                                                 cAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCACTGCGgcgc  >  2:281981/1‑147 (MQ=255)
                                                                                                                                   cggtGGTAGCCTGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTgg                           >  3:73434/4‑120 (MQ=255)
                                                                                                                                   cggtGGTAGCCTGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTgg                           <  4:73434/117‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                     
CGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACTGGCATCGGTCACCGGGTAGATCATTACCTGCAGGCGCGGTGCCGGCAGCTCGCGCTGCGCAGCCAACTGGTTGGCCAGAATGGTAGCCAGGCTACCACCGACACTGTCGCCCACCACCGCCAGTCGCTGCGCATCGATGCCCAGCGCCTCGGCCTGCTCGACCAACCAGCGCCAGGCATCCAGGGCATCGTCACTGGCTGTCGGGAAACGCCACTGCGGCGC  >  NC_002947/4768161‑4768436

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