Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,125 G→A 14.2% A113V (GCC→GTC)  eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,1250GA14.2% 119.5 / 9.2 49A113V (GCC→GTC) ecoecotin
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base A (5/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCT  >  NC_002947/3455991‑3456269
                                                                                                                                      |                                                                                                                                                
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)
                                                                                                                                      |                                                                                                                                                
CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCT  >  NC_002947/3455991‑3456269

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