Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,415,686 C→A 20.3% intergenic (‑287/‑155) kil ← / → sieB inhibitor of ftsZ, killing protein/phage superinfection exclusion protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,415,6860C→A20.3% 3.6 ‑4.6 0.033intergenic (‑287/‑155)kil/sieBinhibitor of ftsZ, killing protein/phage superinfection exclusion protein
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (14/11):  total (14/19)
Fisher's exact test strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AACCGGAGGCGCACTCCCACCATTTAAATTTAACAGACAAGACCGACTCTTTATGGATATCGGAAATGC  >  REL606/1415653‑1415721
                                 |                                   
aaCCGGAGGCGCACTCCCACCATTTAAATTTAACAg                                   >  1:5182422/1‑36 (MQ=255)
 aCCGGAGGCGCACTCCCACCAGTTAAATTTAACAGa                                  >  1:2012827/1‑36 (MQ=255)
  ccGGAGGCGCACTCCCACCATTTAAATTTAACAGaa                                 >  1:4619992/1‑35 (MQ=255)
  ccGGAGGCGCACTCCCACCATTTAAATTTAACAGAc                                 <  1:603810/36‑1 (MQ=255)
  ccGGAGGCGCACTCCCACCATTTAAATTTAACAGAc                                 <  1:4309717/36‑1 (MQ=255)
  ccGGAGGCGCACTCCAACCATTTAAATTTAACAGAc                                 <  1:4597422/36‑1 (MQ=255)
   cGGAGGCGCACTCCCACCCTTTAAATTTAACAGACa                                >  1:5855378/1‑36 (MQ=255)
   cGGAGGCGCACTCCCACCATTTAAATTTAACAGACa                                <  1:7290682/36‑1 (MQ=255)
   cGGAGGCGCACTCCCACCATTTAAATTTAACAGACa                                >  1:2424867/1‑36 (MQ=255)
   cGGAGGCGCACTCCCACCATTTAAATTTAACAGACa                                >  1:1142933/1‑36 (MQ=255)
    ggAGGCGCACTCCCACCATTTAAATTTAACAGACaa                               >  1:3667441/1‑36 (MQ=255)
     gAGGCGCACTCCCACCATTTAAATTTAACAGACAAg                              <  1:3952955/36‑1 (MQ=255)
       ggCGCACTCCCACAATTTAAATTTAACAGACAAGAc                            <  1:3559377/36‑1 (MQ=255)
        gcgcACTCCCACCATTTAAATTTAACAGACAAGAcc                           <  1:7179623/36‑1 (MQ=255)
        gcgcACTCCAACCATTTAAATTTAACAGACAAGAcc                           <  1:4941126/36‑1 (MQ=255)
          gcACTCCCACCATTTAAATTTAACAGACAAGACCGa                         <  1:2492044/36‑1 (MQ=255)
          gcACTCCCACCATTTAAATTTAACAGACAAGACCGa                         >  1:7008166/1‑36 (MQ=255)
           cACTCCCACCATTTAAATTTAACAGACAAGACCGAc                        >  1:3604724/1‑36 (MQ=255)
           cACTCCCACAATTTAAATTTAACAGACAAGACCGAc                        <  1:7179559/36‑1 (MQ=255)
               cccACAATTTAAATTTAAAAGAAAAGACCGACTCtt                    <  1:7737179/36‑1 (MQ=25)
                cgtccaTTCAAATTGCATACACAAGACCGACTCttt                   <  1:524807/33‑1 (MQ=255)
                 caccaTTTAAATTTAACAGACAAGACCGACTCTTTa                  >  1:1127324/1‑36 (MQ=255)
                      ttaaaattaaaAAGAAAAGACCGACTCTTTATGGat             <  1:4678471/33‑1 (MQ=16)
                       ttaaatttaaCAGACAAGACCGACTCTTTATGGata            >  1:2308318/1‑36 (MQ=255)
                       ttaaattaaaAAGAAAAGACCGACTCTTTATGGata            <  1:6059627/36‑1 (MQ=25)
                       taaaatttaaCAGAAAAGACCGACTCTTTATGGata            <  1:815297/34‑1 (MQ=38)
                        taaattaaaAAGAAAAGACCGACTCTTTATGGatat           <  1:2732129/36‑1 (MQ=25)
                        aaaatttaaAAGAAAAGACCGACTCTTTATGGatat           <  1:1149434/35‑1 (MQ=37)
                        aaaattaaaAAGAAAAGACCGACTCTTTATGGatat           <  1:4841795/35‑1 (MQ=21)
                          aattaaaAAGAAAAGACCGACTCTTTATGGATATCg         <  1:3441673/36‑1 (MQ=25)
                             ttaaCAGACAAGACCGACTCTTTATGGATATCGGaa      >  1:7615692/1‑36 (MQ=255)
                                aCAGACAAGACCGACTCTTTATGGATATCGGAAATg   >  1:4068328/1‑36 (MQ=255)
                                aCAGACAAGACCGACTCTTTATGGATATCAGAAAt  >  1:5956181/1‑35 (MQ=255)
                                 cAGACAAGACCGACTCTTTATGGATATCGGAAACgc  >  1:6271120/1‑36 (MQ=255)
                                 |                                   
AACCGGAGGCGCACTCCCACCATTTAAATTTAACAGACAAGACCGACTCTTTATGGATATCGGAAATGC  >  REL606/1415653‑1415721

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG