Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,102,264 G→T 14.4% Q187K (CAA→AAA)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,102,2640G→T14.4% 3.0 ‑4.5 0.054559ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  new base (9/1):  ref base (18/26):  total (27/27)
Fisher's exact test strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTCA  >  REL606/4102231‑4102298
                                 |                                  
aCTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTTCg                                  >  1:2438082/1‑36 (MQ=255)
 cTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGa                                 <  1:4609823/36‑1 (MQ=255)
 cTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGa                                 <  1:560310/36‑1 (MQ=255)
  tGGCGCGGTTGCGCCTGTACCGGTGCGGCTTTCggc                                >  1:3941422/1‑34 (MQ=37)
  tGGCGCGGTTGCGCCTGTACCGGTGCGGCTTTCggc                                >  1:3248722/1‑34 (MQ=37)
  tGGCGCGGTTGCGCCTGTACCGGTGCGGCTTTCGAc                                >  1:1655322/1‑36 (MQ=255)
   ggCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACg                               <  1:476003/36‑1 (MQ=255)
   ggCGCGGTTGCGACTGTAACGGAGCGGCTTTCgcgg                               >  1:8111832/1‑33 (MQ=16)
   ggCGCAGTTGCGCCTGTGCCGGTGCGGCTTGCGACg                               <  1:6626288/36‑1 (MQ=37)
    gcgcGGTTGCGCCTGTACCGGTGCGGCTTGTGGCGt                              >  1:527976/1‑36 (MQ=37)
    gcgcGGTTGCGCCTGTACCGGTGCGGCTTGCGACGt                              <  1:5089019/36‑1 (MQ=255)
    gcgcGGTTGCGCCTGTACCGGTGCGGCTTGCGACGt                              <  1:1951923/36‑1 (MQ=255)
      gcgGTTGCGCCTGTACCGGTGCGGCTTTCGACGTag                            >  1:6350279/1‑35 (MQ=255)
      gccgTTGCGCCTGTACCGTTGCGGCTTGCGACGTAc                            <  1:4699909/33‑1 (MQ=37)
       cgGTTGCGCCTGTACCGGTGCGGCTTGCGACGTAcg                           <  1:7808990/36‑1 (MQ=255)
       cgGTTGCGCCTGTACCGGTGCGGCTTGCGACGTAcg                           >  1:4287175/1‑36 (MQ=255)
       cgGTTGCGCCTGTACCGGTGCGGCTTGCGACGTAcg                           >  1:4129897/1‑36 (MQ=255)
        ggTTGCGCCTGTACCGGTGCGGCTTGCGACGTAcgc                          >  1:199060/1‑36 (MQ=255)
            gcgcCTGTACCGGTGCGGCTTTCGACGTACGCGTct                      >  1:6425863/1‑36 (MQ=255)
            gcgcCTGTACCGGTGCGGCTTGCGACGTACGCGTct                      >  1:7271563/1‑36 (MQ=255)
            gcgcCTGTACCGGCGAGGCTTTCGCCGTACGCGTct                      >  1:7526896/1‑36 (MQ=16)
             cgcCTGTACCGGTGCGGCTTGCGACGTACGCGTctg                     <  1:403006/36‑1 (MQ=255)
              gcCTGTACCGGTGCGGCTTGCGACGTACGCGTctgc                    <  1:6151810/36‑1 (MQ=255)
               ccTGTACCGGTGCGGCTTGCGACGTACGCGTctgct                   >  1:5512364/1‑36 (MQ=255)
               ccTGTACCGGTGCGGCTTGCGACGTACGCGTctgct                   <  1:7999128/36‑1 (MQ=255)
                cTGTACCGGTGCGGCTTGCGACGTACGCGTctgctg                  <  1:4043965/36‑1 (MQ=255)
                 tGTACCGGTGCGGCTTTCGACGTACGCGTctgctgc                 >  1:6269904/1‑36 (MQ=255)
                 tGTACCGGTGCGGCTTGCGACGTACGCGTctgctgc                 <  1:6157405/36‑1 (MQ=255)
                 tGTACCGGTGCGGCTTGCCACGTACGCGTctgcttt                 >  1:5070748/1‑34 (MQ=38)
                  gTACCGGTTCGGCTTGCGCCGTACGCGTctgctgct                <  1:3464738/36‑1 (MQ=37)
                  gTACCGGTGCGGCTTGCGACGTACGCGTctgctgct                <  1:2011373/36‑1 (MQ=255)
                   tACCGGTGCGGCTTGCGACGTACGCGTctgctgctg               <  1:2251981/36‑1 (MQ=255)
                    aCCGGTGCGGCTTGGGACGGAAGGGTCtactgctgc              >  1:5277447/1‑36 (MQ=16)
                    aCCGGTGCGGCTTGCGACGTACGCGTCtgctgctgc              >  1:1762519/1‑36 (MQ=255)
                    aCCGGTGCGGCTTGCGACGTACGCGTAtgcttctgc              >  1:6740411/1‑36 (MQ=255)
                     ccGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGcc             >  1:490239/1‑36 (MQ=255)
                       ggTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAg           <  1:4031651/36‑1 (MQ=255)
                        gTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAga          >  1:6581060/1‑35 (MQ=255)
                        gTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGc          <  1:3281022/36‑1 (MQ=255)
                        gTGCGGCTTGCGACGTACGCGTCTGCTGATGCCAGc          >  1:2304693/1‑36 (MQ=255)
                          gCGGCTTTCGACGTACGCGTCTGCTGCTGCCAGCtt        <  1:6737537/36‑1 (MQ=255)
                          gCGGCTTGCGACGTACGCGTCTGCTGCTGCCCGCtt        >  1:1967892/1‑36 (MQ=255)
                          gCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGCtt        <  1:2099285/36‑1 (MQ=255)
                          gCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGCtt        >  1:7091611/1‑36 (MQ=255)
                          gCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGCtt        >  1:4232701/1‑36 (MQ=255)
                             gCTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGt     <  1:8076794/36‑1 (MQ=255)
                              cTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGtt    <  1:6568230/36‑1 (MQ=255)
                              cTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGtt    <  1:25996/36‑1 (MQ=255)
                              cTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGtt    <  1:761544/36‑1 (MQ=255)
                               ttGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTc   <  1:7104685/36‑1 (MQ=255)
                               ttGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTc   <  1:4685941/36‑1 (MQ=255)
                                tGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTCa  >  1:5841256/1‑36 (MQ=255)
                                tGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTCa  <  1:1911064/36‑1 (MQ=255)
                                tGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTCa  >  1:8286161/1‑36 (MQ=255)
                                 |                                  
ACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAGCTTTGTTCA  >  REL606/4102231‑4102298

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG