Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,493,865 T→C 19.3% E22G (GAA→GGA)  ECB_03279 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,493,8650T→C19.3% 3.6 ‑3.2 0.04665ECB_03279hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/10):  ref base (13/23):  total (13/33)
Fisher's exact test strand distribution p-value = 4.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGCTCTCTTCTTCCCTAACGCTGTTATAGTACCCTTCCCGGTCGACACATCCCGCCAGTAACGCCGTGCCT  >  REL606/3493830‑3493900
                                   |                                   
ggCTCTCTTCTTCCCTAACGCTGTTATAGTACCCtt                                     <  1:1833720/36‑1 (MQ=255)
 gCTCTCTTCTTCCCTAACGCTGTTATAGTACCCTTc                                    <  1:3208319/36‑1 (MQ=255)
  ctctctTCTTCCCTAACGCTGTTATAGTACCCTTcc                                   <  1:94703/36‑1 (MQ=255)
   tctctTCTTCCCTAACGCTGTTATAGTACCCTTccc                                  <  1:3276341/36‑1 (MQ=255)
   tctctTCTTCCCTAACGCTGTTATAGTACCCTTccc                                  <  1:2886888/36‑1 (MQ=255)
   tctctTCTTCCCTAACGCTGTTATAGTACCCTTccc                                  >  1:5212053/1‑36 (MQ=255)
    ctctTCTTCCCTAACGCTGTTATAGTACCCTTCCCg                                 >  1:2443604/1‑36 (MQ=255)
      cttcttCCCTAACGCTGTTATAGTACCCTTCCCGGt                               <  1:4352438/36‑1 (MQ=255)
       ttcttcCCTAACGCTGTTATAGTACCCTTCCCGGTc                              <  1:7815186/36‑1 (MQ=255)
       ttcttcCCTAACGCTGTTATAGTACCCTTCCCGGTc                              >  1:2585245/1‑36 (MQ=255)
         cttcCCTAACGCTGTTATAGTACCCTTCCCGGTCGa                            <  1:7218120/36‑1 (MQ=255)
          ttcCCTAACGCTGTTATAGTACCCTTCCCGGTCGac                           >  1:2746405/1‑36 (MQ=255)
          ttcCCTAACGCTGTTATAGTACCCTTCCCGGTCGac                           <  1:6243898/36‑1 (MQ=255)
          ttcCCTAACGCTGTTATAGTACCCTTCCCGGTCGac                           <  1:6245842/36‑1 (MQ=255)
           tcCCTAACGCTGTTATAGTACCCTTCCCGGTCGaca                          >  1:338276/1‑36 (MQ=255)
           tcCCTAACGCTGTTATAGTACCCTTCCCGGTCGaca                          <  1:7852717/36‑1 (MQ=255)
           tcCCTAACGCTGTTATAGTACCCTTCCCGGTCGaca                          <  1:4289266/36‑1 (MQ=255)
            cccTAACGCTGTTATAGTACCCTTCCCGGTCGacac                         <  1:6888633/36‑1 (MQ=255)
                aaCGCTGTTACAGTCCCCATCCCGGTCGACACATcc                     <  1:4525543/36‑1 (MQ=25)
                 aCGCTGTTATAGTACCCTTCCCGGTCGACACATccc                    <  1:1855841/36‑1 (MQ=255)
                 aCGCTGTTATAGTACCCTTCCCGGTCGACACATccc                    <  1:5927763/36‑1 (MQ=255)
                 aCGCTGTTATAGTACCCTTCCCGGTCGACACATccc                    >  1:5732752/1‑36 (MQ=255)
                 aCGCTGTTATAGTACCCTCCCCGGTCGACACATccc                    <  1:4432327/36‑1 (MQ=255)
                   gCTGTTATAGTACCCTTCCCGGTCGACACATCCCGc                  <  1:6959671/36‑1 (MQ=255)
                    cTGTTATAGTACCCTTCCCGGTCGACACATCCCGcc                 <  1:4426630/36‑1 (MQ=255)
                    cTGTTATAGTACCCTTCCCGGTCGACACATCCCGcc                 >  1:2197212/1‑36 (MQ=255)
                     tGTTATAGTACCCTTCCCGGTCGACACATCCCGCCa                <  1:276346/36‑1 (MQ=255)
                       ttATAGTCCCCTTCCCGGTCGACACATCCCGCCAGt              <  1:20165/36‑1 (MQ=255)
                       ttATAGTACCCTTCCCGGTCGACACATCCCGCCAGt              <  1:6943678/36‑1 (MQ=255)
                        tataGTACCCTTCCCGGTCGACACATCCCGCCAGTa             <  1:6510724/36‑1 (MQ=255)
                         ataGTACCCTTCCCGGTCGACACATCCCGCCAGTaa            >  1:1051208/1‑36 (MQ=255)
                           aGTACCCTTCCCGGTCGACACATCCCGCCAGTAACg          >  1:4836057/1‑36 (MQ=255)
                           aGTACCCTTCCCGGTCGACACATCCCGCCAGTAACg          >  1:6315021/1‑36 (MQ=255)
                           aGTACCCTTACCGGTCGACACATCCCGCCTGTAACg          >  1:6946074/1‑36 (MQ=37)
                           aGTACCCTCCCCGGTCGACACACCCCGCCAGTAACg          <  1:3321532/36‑1 (MQ=38)
                            gTACCCTCCCCGGTCGACACATCCCGCCAGTAACGc         <  1:36728/36‑1 (MQ=255)
                             tccccTCCCCGGTCGACACATCCCGCCAGTAACGcc        <  1:545371/34‑1 (MQ=37)
                              ccccTCCCCTGTCGACACATCCCGCCAGTAACGCCg       <  1:3633966/35‑1 (MQ=25)
                              ccccTCCCCGGTCGACACATCCCGCCAGTAACGCCg       <  1:5549187/35‑1 (MQ=255)
                              ccccTCCCCGGTCGACACATCCCGCCAGTAACGCCg       <  1:1814477/35‑1 (MQ=255)
                              ccccTCCCCGGTCGACACATCCCGCCAGTAACGCCg       <  1:1344672/35‑1 (MQ=255)
                               cccTTCCCGGTCGACACATCCCGCCAGTAACGCCgc      >  1:416428/1‑35 (MQ=255)
                               cccTTCCCGGTCGACACATCCCGCCAGTAACGCCGt      >  1:6601311/1‑36 (MQ=255)
                               cccTCCCCGGTGAACACACCCCGCCAGTAACGCCGt      <  1:4990183/36‑1 (MQ=16)
                               cccTCCCCGGGCGACAAACCCCTCCAGTAACGCCGt      <  1:2630738/36‑1 (MQ=16)
                               cccCTCCCAGCCGCCACATCCCGCCAGTAACGCCGt      <  1:6799377/36‑1 (MQ=16)
                                 cTTCCCGTTCGACACATCCCGCCAGTAACGCCGTGc    <  1:6078062/36‑1 (MQ=255)
                                   tCCCGGTCGACACATCCCCCCAGTAACGCCGTGCCt  >  1:500721/1‑36 (MQ=255)
                                   |                                   
GGCTCTCTTCTTCCCTAACGCTGTTATAGTACCCTTCCCGGTCGACACATCCCGCCAGTAACGCCGTGCCT  >  REL606/3493830‑3493900

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG