Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,128 T→G 13.8% E112A (GAG→GCG) ‡ eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,1280TG13.8% 118.2 / 14.9 51E112A (GAG→GCG) ‡ecoecotin
Reads supporting (aligned to +/- strand):  ref base T (26/18);  new base G (5/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCG  >  NC_002947/3455991‑3456275
                                                                                                                                         |                                                                                                                                                   
ctTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTg                                                                                                                                          <  2:169911/149‑1 (MQ=255)
 tttAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGc                                                                                                                                         >  2:84322/1‑149 (MQ=255)
      ccTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCggtg                                                                                                                                    <  3:337876/149‑4 (MQ=255)
        ttCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCgg                                                                                                                                  >  3:148628/1‑149 (MQ=255)
              acaTCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGc                                                                                                                            >  6:12375/1‑149 (MQ=255)
                 tCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCAt                                                                                                                         >  7:279247/1‑149 (MQ=255)
                 tCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCAt                                                                                                                         >  7:225494/1‑149 (MQ=255)
                                cacaCGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccgg                                                                                                          >  5:8584/1‑149 (MQ=255)
                                    cgcgGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccgggccg                                                                                                      <  3:79090/149‑1 (MQ=255)
                                         tAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGAc                                                                                                 <  6:8584/149‑1 (MQ=255)
                                               aCTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCgg                                                                                                                                      >  4:276780/1‑105 (MQ=38)
                                               aCTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCgg                                                                                                                                      <  3:276780/106‑2 (MQ=38)
                                                     aCATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCg                                                                                     <  8:279247/149‑1 (MQ=255)
                                                     aCATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCg                                                                                     <  8:225494/149‑1 (MQ=255)
                                                           ttGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCAATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGgtagt                                                                                >  7:3732/1‑148 (MQ=255)
                                                                   aTACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGc                                                                       >  8:109796/1‑149 (MQ=255)
                                                                     acacCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                     >  5:32208/1‑149 (MQ=255)
                                                                     acacAACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                     >  8:188860/1‑149 (MQ=255)
                                                                       accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                     >  6:219666/1‑147 (MQ=255)
                                                                       accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                     >  6:219663/1‑147 (MQ=255)
                                                                       accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCa                                                                   <  5:219663/149‑1 (MQ=255)
                                                                       accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCa                                                                   <  5:219666/149‑1 (MQ=255)
                                                                           ccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCt                                                               >  8:91002/1‑149 (MQ=255)
                                                                                 gCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAgg                                                         >  7:71724/1‑149 (MQ=255)
                                                                                        gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTg                                                                                                                  <  4:40380/85‑1 (MQ=255)
                                                                                        gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTg                                                                                                                  >  3:40380/1‑85 (MQ=255)
                                                                                        gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGgtgtgt                                                   <  7:188860/148‑1 (MQ=255)
                                                                                          tgttgtAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGc                                                 >  3:44659/1‑148 (MQ=255)
                                                                                                aGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAgg                                          >  8:186329/1‑149 (MQ=255)
                                                                                                  cgcAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGtt                                        <  5:12375/149‑1 (MQ=255)
                                                                                                   gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                       <  8:3732/149‑1 (MQ=255)
                                                                                                   gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                       <  7:109796/149‑1 (MQ=255)
                                                                                                   gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCAGACGATGTACACGCGGTAGTAATTGAAGCACCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                       >  2:272318/1‑149 (MQ=255)
                                                                                                           aaGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCCGCGTGCTGGCCGGGCCGCTGACGTTGTCCCGGCGGTAGGAAGTGTAGCCCCAGCCTTCCAGGGTGTGTGGTTCCAGGTTGCCGccccg                               >  8:216848/1‑147 (MQ=255)
                                                                                                             gCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccggg                                                                                                         >  8:116428/1‑73 (MQ=255)
                                                                                                               cATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgcgc                           >  5:307295/1‑149 (MQ=255)
                                                                                                                   gCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGAc                       <  6:307295/149‑1 (MQ=255)
                                                                                                                    ccGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGAc                       >  7:210342/1‑148 (MQ=255)
                                                                                                                    ccGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACg                      <  6:32208/149‑1 (MQ=255)
                                                                                                                     cGACCACCGGCAGGAAGGCCTCGGTCTTCTTGCCGTCCGGTCAGGCCATCAGCGTGCTGGCCGGGGCGCTGAGGTTGTCCAGGCGGTAGTTATTGTAGCCCCAGCCTTGCTGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGa                     <  1:272318/149‑1 (MQ=255)
                                                                                                                     cGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGAAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGa                     <  4:148628/149‑1 (MQ=255)
                                                                                                                      gACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGAt                    >  8:259249/1‑149 (MQ=255)
                                                                                                                        cggacGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTg                  >  3:42377/5‑149 (MQ=255)
                                                                                                                        cggacGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTg                  >  3:218950/5‑149 (MQ=255)
                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                             <  7:178248/135‑1 (MQ=255)
                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                             >  8:178248/1‑135 (MQ=255)
                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAg               >  5:317587/1‑149 (MQ=255)
                                                                                                                            cGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGt              >  2:221966/1‑149 (MQ=255)
                                                                                                                                  gAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCt        >  1:328314/1‑149 (MQ=255)
                                                                                                                                       ccTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGc   <  8:210342/149‑1 (MQ=255)
                                                                                                                                        cTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCg  >  3:78875/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCG  >  NC_002947/3455991‑3456275

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