Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,141,329 Δ1 bp intergenic (‑84/+110) ECB_02035 ← / ← ECB_02036 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,141,3290A.96.3% 246.9 / 3.7 54intergenic (‑84/+110)ECB_02035/ECB_02036hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (19/33);  total (20/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TCTTTGCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGATCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGATTATTT  >  REL606/2141254‑2141396
                                                                           |                                                                   
tctTTGCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGa                                                                     >  2:10112/1‑76 (MQ=255)
    tGCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAAc                                                                >  1:1500255/1‑76 (MQ=255)
     gCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACg                                                               >  1:1783168/1‑76 (MQ=255)
     gCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACg                                                               <  1:1098212/76‑1 (MQ=255)
      cTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGATCAACg                                                               >  2:1650530/1‑76 (MQ=255)
        aaGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCt                                                            >  2:972471/1‑76 (MQ=255)
        aaGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCt                                                            >  2:197825/1‑76 (MQ=255)
        aaGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCt                                                            >  2:215454/1‑76 (MQ=255)
         aGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTc                                                           >  2:263224/1‑76 (MQ=255)
              cGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAg                                                      <  1:1376000/76‑1 (MQ=255)
               gCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGt                                                     <  1:1192622/76‑1 (MQ=255)
                 tGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTAt                                                   <  1:453818/76‑1 (MQ=255)
                  gTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATc                                                  >  1:1003255/1‑76 (MQ=255)
                  gTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATc                                                  >  2:1761366/1‑76 (MQ=255)
                   ttATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTACGTATCt                                                 >  2:572193/1‑76 (MQ=255)
                   ttATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCt                                                 <  1:1432634/76‑1 (MQ=255)
                    tATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtg                                                <  2:294674/76‑1 (MQ=255)
                     aTGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtgt                                               <  1:1107122/76‑1 (MQ=255)
                      tGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtgtt                                              <  1:1297039/76‑1 (MQ=255)
                      tGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtgtt                                              <  2:780011/76‑1 (MQ=255)
                        aGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtgttgt                                            >  1:482302/1‑76 (MQ=255)
                        aGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCtgttgt                                            >  2:993774/1‑76 (MQ=255)
                         gTTCCGATGCCAGAGCGGCAGCGGGTTTATCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTtgtg                                           <  2:1261937/76‑1 (MQ=255)
                         gTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTtgtg                                           <  2:275056/76‑1 (MQ=255)
                         gTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTtgtg                                           <  2:1538529/76‑1 (MQ=255)
                         gTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTtgtg                                           <  2:1261926/76‑1 (MQ=255)
                          ttCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGATCAACGCCTCGTAAGTATCTGTtgtg                                           <  1:10112/76‑1 (MQ=255)
                           tCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGaa                                         <  1:812299/76‑1 (MQ=255)
                           tCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGaa                                         <  1:65839/76‑1 (MQ=255)
                           tCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGaa                                         >  1:722502/1‑76 (MQ=255)
                            ccGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGaaa                                        <  1:1184320/76‑1 (MQ=255)
                            ccGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGaaa                                        <  1:1184326/76‑1 (MQ=255)
                              gATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAAtt                                      <  2:1442560/76‑1 (MQ=255)
                                   cAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAttt                                 <  1:496433/76‑1 (MQ=255)
                                   cAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAttt                                 <  2:1032629/76‑1 (MQ=255)
                                   cAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAttt                                 <  2:755841/76‑1 (MQ=255)
                                    agagCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAtttt                                <  1:1450833/76‑1 (MQ=255)
                                    agagCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAtttt                                <  2:283270/76‑1 (MQ=255)
                                    agagCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGAtttt                                <  2:1671282/76‑1 (MQ=255)
                                      agCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCa                              <  2:314630/76‑1 (MQ=255)
                                        cGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATc                            >  2:1072221/1‑76 (MQ=255)
                                               gggTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCtt                     >  1:554677/1‑76 (MQ=255)
                                               gggTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCtt                     <  1:509468/76‑1 (MQ=255)
                                                ggTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTa                    <  2:556302/76‑1 (MQ=255)
                                                 gTTTGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTAt                   >  1:841244/1‑76 (MQ=255)
                                                    tGTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCt                >  1:703315/1‑76 (MQ=255)
                                                     gTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCtt               <  2:1094176/76‑1 (MQ=255)
                                                     gTCTTCTTTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCtt               >  1:1060515/1‑76 (MQ=255)
                                                       cttcttTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCg             >  1:803125/1‑76 (MQ=255)
                                                       cttcttTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCg             <  1:337700/76‑1 (MQ=255)
                                                       cttcttTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCg             <  1:337693/76‑1 (MQ=255)
                                                         tcttTGTTGGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTg           >  1:948742/1‑76 (MQ=255)
                                                              gttgGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGatt      >  1:625819/1‑76 (MQ=255)
                                                               ttgGCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGatta     >  2:1555994/1‑76 (MQ=255)
                                                                  gCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGAttattt  <  2:1431788/76‑1 (MQ=255)
                                                                  gCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGAttattt  <  2:97845/76‑1 (MQ=255)
                                                                  gCACGGTTGTCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGAttattt  <  1:981301/76‑1 (MQ=255)
                                                                           |                                                                   
TCTTTGCTAAGGATCGCTGTTATGAGTTCCGATACCAGAGCGGCAGCGGGTTTGTCTTCTTTGTTGGCACGGTTGATCAACGCCTCGTAAGTATCTGTTGTGAAATTGATTTTCATCGTGTCTTATCCTTCGTGAGATTATTT  >  REL606/2141254‑2141396

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