Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,007,369 A→T L154* (TTG→TAG)  ECB_02804 ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,007,3690AT97.3% 136.2 / ‑3.3 37L154* (TTG→TAG) ECB_02804conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (17/19);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TTACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCAACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGAAAACCC  >  REL606/3007298‑3007440
                                                                       |                                                                       
ttACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGctacaa                                                                     <  2:980646/76‑1 (MQ=255)
ttACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGctacaa                                                                     <  1:516107/76‑1 (MQ=255)
ttACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGctacaa                                                                     >  1:737570/1‑76 (MQ=255)
ttACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGctacaa                                                                     >  2:1229531/1‑76 (MQ=255)
ttACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGctacaa                                                                     >  2:97601/1‑76 (MQ=255)
 tACGCCCTGCAGCATCCGGAATGCCGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCtacaac                                                                    >  2:8865/1‑76 (MQ=255)
 tACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCtacaac                                                                    >  2:1336322/1‑76 (MQ=255)
 tACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCtacaac                                                                    >  1:264395/1‑76 (MQ=255)
     cccTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAActttc                                                                <  1:931849/76‑1 (MQ=255)
       cTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAActttcct                                                              >  1:281708/1‑76 (MQ=255)
        tGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAActttcctt                                                             <  1:1281658/76‑1 (MQ=255)
         gcagcaTCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAActttccttt                                                            <  2:170219/76‑1 (MQ=255)
           agcaTCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCt                                                          <  2:776778/76‑1 (MQ=255)
                ccGGAATGCAGCGCCATTCTGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCg                                                     >  2:205632/1‑76 (MQ=255)
                   gAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAg                                                  <  2:1073574/76‑1 (MQ=255)
                    aaTGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGc                                                 >  1:1264360/1‑76 (MQ=255)
                    aaTGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGATACAAATTTCATTTCTGCCAGGAga                                                 >  2:1525791/1‑75 (MQ=255)
                     aTGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCg                                                >  2:827107/1‑76 (MQ=255)
                        cAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATg                                             >  2:26247/1‑76 (MQ=255)
                         aGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGa                                            <  2:983292/76‑1 (MQ=255)
                            gcCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGc                                         <  1:428036/76‑1 (MQ=255)
                              cATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGccc                                       <  1:1486945/76‑1 (MQ=255)
                               aTTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCg                                      >  1:292719/1‑76 (MQ=255)
                                 tATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGc                                    >  1:1049985/1‑76 (MQ=255)
                                  aTGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCg                                   >  1:118093/1‑76 (MQ=255)
                                  aTGCGCATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCg                                   <  1:1085754/76‑1 (MQ=255)
                                     cgcATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAgc                                <  2:1220473/76‑1 (MQ=255)
                                      gcATCATCGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAgcg                               >  2:923755/1‑76 (MQ=255)
                                         tcatcGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAg                            <  2:1348366/76‑1 (MQ=255)
                                           atcGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGtc                          <  2:406303/76‑1 (MQ=255)
                                            tcGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGtct                         <  2:977839/76‑1 (MQ=255)
                                             cGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGtctc                        <  1:1507216/76‑1 (MQ=255)
                                             cGACTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGtctc                        >  1:1366214/1‑76 (MQ=255)
                                              naCTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCa                       >  1:1098359/2‑76 (MQ=255)
                                                cTGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCAtt                     >  1:466213/1‑76 (MQ=255)
                                                 tGAATGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTc                    <  1:380772/76‑1 (MQ=255)
                                                   aaTGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTg                  <  1:215064/76‑1 (MQ=255)
                                                   aaTGGAGAGAAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTg                  <  2:756716/76‑1 (MQ=255)
                                                          agaAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGccc           >  1:1421501/1‑76 (MQ=255)
                                                          agaAGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGccc           <  1:114859/76‑1 (MQ=255)
                                                          agaAGCACAATGCAACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGccc           >  2:475839/1‑76 (MQ=255)
                                                             aGCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGa        >  2:284277/1‑76 (MQ=255)
                                                              gCACAATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGaa       >  1:1520293/1‑76 (MQ=255)
                                                               cacaATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGaca      >  2:285963/1‑74 (MQ=255)
                                                               cacaATGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGaaa      >  2:407595/1‑76 (MQ=255)
                                                                  aaTGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGAAAAcc   <  1:511265/76‑1 (MQ=255)
                                                                   aTGCTACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGAAAAccc  >  1:612891/1‑76 (MQ=255)
                                                                       |                                                                       
TTACGCCCTGCAGCATCCGGAATGCAGCGCCATTATGCGCATCATCGACTGAATGGAGAGAAGCACAATGCAACAACTTTCCTTTCTGCCCGGAGCGATGACGCCCGGCGAGCGCAGTCTCATTCTGCGGGCCCTGAAAACCC  >  REL606/3007298‑3007440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: