Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,456,312 C→A 9.0% E51* (GAG→TAG)  eco ← ecotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,456,3120CA 9.0% 155.7 / 2.4 51E51* (GAG→TAG) ecoecotin
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGG  >  NC_002947/3456174‑3456452
                                                                                                                                          |                                                                                                                                            
cgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTgg                                                                                                                                    >  4:70491/1‑149 (MQ=255)
   tGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGACGACTCGTCGGCCTGCGcc                                                                                                                                 >  7:37551/1‑149 (MQ=255)
       gTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCa                                                                                                                             <  7:253247/149‑1 (MQ=255)
       gTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCa                                                                                                                             <  2:207773/149‑1 (MQ=255)
           tCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTGCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTg                                                                                                                         >  6:12583/1‑149 (MQ=255)
           tCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGt                                                                                                                          <  5:316145/148‑1 (MQ=255)
             cAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGa                                                                                                                       >  3:204145/1‑149 (MQ=255)
               ggCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATa                                                                                                                     <  7:301698/149‑1 (MQ=255)
                      gtaATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGcc                                                                                                              >  6:37233/1‑149 (MQ=255)
                         aTTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCggg                                                                                                           <  4:44761/149‑1 (MQ=255)
                             tAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGaa                                                                                                       <  5:37233/149‑1 (MQ=255)
                               gCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAg                                                                                                      <  7:214329/148‑1 (MQ=255)
                                  ccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGTCGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCtt                                                                                                  >  1:48967/1‑149 (MQ=255)
                                     gCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCttttt                                                                                               <  3:70491/149‑1 (MQ=255)
                                      ccTTCCAGGGTGTGTGCTTCCTGGTTGCCGCCGAGGCTCTTACGATTGAAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTc                                                                                              <  4:188206/149‑1 (MQ=255)
                                            aGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTAGACCTGCAGCGTCTTTACCGCGAGGAGTTCCAGTTTAGAGGACGACTAGTCTGCCTGGGCCGGCAGGGGGATAACATGCCGGGTGAAGCCTTGTTCCgacgc                                                                                        >  4:9188/1‑145 (MQ=255)
                                                 gtgtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgggg                                                                                    >  5:304481/1‑148 (MQ=255)
                                                    tgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAgg                                                                                <  6:304481/149‑1 (MQ=255)
                                                       ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACGAGTCGGCCtataa                                                                                                                                  >  7:42177/1‑92 (MQ=255)
                                                       ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACGAGTCGGCCtataa                                                                                                                                  <  8:42177/96‑5 (MQ=255)
                                                        tCCAGGTTG‑CGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgcgc                                                                            >  6:20043/1‑148 (MQ=255)
                                                          ccgggtGCCGCCCCGGCGCGGACGATTGCAGGCGACCTGCAGCGTCTTGCCAGCGAGGCTTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATCCCCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                          <  3:9188/144‑1 (MQ=255)
                                                          cAGGTTGCCGCCCAGGCGCTGTCGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                          <  2:48967/149‑1 (MQ=255)
                                                          cAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGa                                                                          <  4:204145/149‑1 (MQ=255)
                                                                     ccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAgg                                                               >  6:81564/1‑149 (MQ=255)
                                                                     ccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAgg                                                               >  4:130134/1‑149 (MQ=255)
                                                                       aGGCGCTGACGATTGCAGTCGACCTACAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCt                                                             >  5:8212/1‑149 (MQ=255)
                                                                                cGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCAt                                                    <  3:248862/149‑1 (MQ=255)
                                                                                 gATTGCAGTCGACCTACAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATc                                                   <  6:8212/149‑1 (MQ=255)
                                                                                   ttGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCGATCg                                                  >  1:115140/1‑148 (MQ=255)
                                                                                        gTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTg                                            >  2:282714/1‑149 (MQ=255)
                                                                                             ccTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcg                                       >  6:317443/1‑149 (MQ=255)
                                                                                                gCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcggcg                                    >  7:42328/1‑149 (MQ=255)
                                                                                                 cAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAAGCTGGCGGCCATCGCAGGTgcggcggc                                     <  5:12583/147‑1 (MQ=255)
                                                                                                      tCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                              <  8:42328/149‑1 (MQ=255)
                                                                                                              gCGAGGATTTCCAGTTTATAGGTCGACTCGTCGGCCTGGGCCGGCAGCTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGc                      >  8:68890/1‑149 (MQ=255)
                                                                                                              gCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGc                      <  5:81564/149‑1 (MQ=255)
                                                                                                                gAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCa                    <  3:130134/149‑1 (MQ=255)
                                                                                                                   gATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGaa                 >  2:11562/1‑149 (MQ=255)
                                                                                                                         cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgcg            <  4:43937/148‑1 (MQ=255)
                                                                                                                         cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgcg            >  3:43937/1‑148 (MQ=255)
                                                                                                                         cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCgg           >  7:219436/1‑149 (MQ=255)
                                                                                                                          aGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGt          >  2:278403/1‑149 (MQ=255)
                                                                                                                                taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgg    <  1:282714/149‑1 (MQ=255)
                                                                                                                                taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg     <  8:219436/148‑1 (MQ=255)
                                                                                                                                taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg     <  1:130121/148‑1 (MQ=255)
                                                                                                                                 aGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg   >  8:212390/1‑149 (MQ=255)
                                                                                                                                 aGGCCGACGAGTCGGCCTGGGCCGGCAGGAGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg   >  1:26092/1‑149 (MQ=255)
                                                                                                                                   gCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgggg  >  7:33938/1‑148 (MQ=255)
                                                                                                                                       aCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCttt                                                                                                 >  7:101070/1‑49 (MQ=255)
                                                                                                                                       aCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCttt                                                                                                 <  8:101070/49‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
CGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGG  >  NC_002947/3456174‑3456452

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