Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,120 C→A 14.3% V115L (GTG→TTG)  eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,1200CA14.3% 133.3 / 12.4 49V115L (GTG→TTG) ecoecotin
Reads supporting (aligned to +/- strand):  ref base C (24/17);  new base A (5/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

GGTCCTGCCTACTCTACTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGT  >  NC_002947/3455975‑3456263
                                                                                                                                                 |                                                                                                                                               
gggccTGCCTACTCTACTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGaa                                                                                                                                              <  3:40419/146‑1 (MQ=255)
                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:225494/1‑149 (MQ=255)
                                 tCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCAt                                                                                                             >  7:279247/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:225494/149‑1 (MQ=255)
                                                                     aCATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCg                                                                         <  8:279247/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                                                                                                      >  3:40380/1‑85 (MQ=255)
                                                                                                        gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTg                                                                                                      <  4:40380/85‑1 (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                           <  7:109796/149‑1 (MQ=255)
                                                                                                                   gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                           <  8:3732/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:218950/5‑149 (MQ=255)
                                                                                                                                        cggacGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTg      >  3:42377/5‑149 (MQ=255)
                                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                 >  8:178248/1‑135 (MQ=255)
                                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                 <  7:178248/135‑1 (MQ=255)
                                                                                                                                           ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAg   >  5:317587/1‑149 (MQ=255)
                                                                                                                                            cGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGt  >  2:221966/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGTCCTGCCTACTCTACTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGT  >  NC_002947/3455975‑3456263

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