Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 795,321 (G)6→7 intergenic (+26/+16) ECB_00736 → / ← ECB_00737 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606795,3151.G90.6% 244.0 / 16.9 64intergenic (+20/+22)ECB_00736/ECB_00737hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (18/40);  total (20/44)
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.13e-01

TTTGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGTGATAA  >  REL606/795241‑795388
                                                                           |                                                                         
tttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAg                                                                          >  1:1017860/1‑76 (MQ=255)
tttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAg                                                                          <  1:910375/76‑1 (MQ=255)
tttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAg                                                                          <  2:640760/76‑1 (MQ=255)
 ttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgg                                                                         >  2:1129123/1‑76 (MQ=255)
 ttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgg                                                                         <  1:1853236/76‑1 (MQ=255)
 ttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgg                                                                         <  1:1853215/76‑1 (MQ=255)
 ttGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgg                                                                         >  2:756096/1‑76 (MQ=255)
  tGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAggg                                                                        >  1:239930/1‑76 (MQ=255)
     ccTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGCTAATACGTTGAAATTATGTGGTGTATTATTAgggggg                                                                     >  1:652526/1‑76 (MQ=255)
       ttAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggt                                                                   <  2:1242300/76‑3 (MQ=255)
       ttAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggt                                                                   <  2:975653/76‑3 (MQ=255)
       ttAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggt                                                                   <  1:825511/76‑3 (MQ=255)
        tAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggta                                                                  >  1:1193254/1‑73 (MQ=255)
        tAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggta                                                                  >  2:931780/1‑73 (MQ=255)
        tAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggta                                                                  <  2:220131/76‑4 (MQ=255)
        tAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAgggggggta                                                                  <  2:1611635/76‑4 (MQ=255)
           ttCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTAccc                                                                <  2:1578130/76‑1 (MQ=255)
           ttCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTAccc                                                                >  2:862830/1‑76 (MQ=255)
               aGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTAccccccn                                                            <  1:1187917/76‑2 (MQ=255)
               aGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCt                                                            <  2:1249894/76‑1 (MQ=255)
                gAATTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGTACCCCCCtaa                                                          <  1:555743/76‑1 (MQ=255)
                gAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCta                                                           <  1:539224/76‑1 (MQ=255)
                gAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCta                                                           <  1:558609/76‑1 (MQ=255)
                gAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCta                                                           <  2:1632431/76‑1 (MQ=255)
                gAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCta                                                           <  2:684004/76‑1 (MQ=255)
                gAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCta                                                           <  2:845617/76‑1 (MQ=255)
                  aaTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaat                                                         <  2:934795/76‑1 (MQ=255)
                  aaTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaat                                                         <  2:1296351/76‑1 (MQ=255)
                  aaTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaat                                                         <  2:1296345/76‑1 (MQ=255)
                  aaTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaat                                                         <  2:934803/76‑1 (MQ=255)
                  aaTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaat                                                         <  2:1134236/76‑1 (MQ=255)
                   aTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaata                                                        <  2:1330563/76‑1 (MQ=255)
                   aTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaata                                                        <  2:703550/76‑1 (MQ=255)
                   aTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaata                                                        <  2:723562/76‑1 (MQ=255)
                   aTAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaata                                                        >  1:895270/1‑76 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       <  2:1022788/76‑1 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       <  2:946756/76‑1 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       <  1:1469723/76‑1 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       <  2:143813/76‑1 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       >  2:303879/1‑76 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCtaataa                                                       <  2:1746090/76‑1 (MQ=255)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGAATTATTAGGGGGGGTACCCCCCtaataa                                                       <  1:1469716/76‑1 (MQ=39)
                    tAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGAATTATTAGGGGGGTACCCCCCTAATaaa                                                      >  1:1113607/1‑76 (MQ=255)
                     aaCCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATaaa                                                      <  2:321118/76‑1 (MQ=255)
                     aaCCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATaaa                                                      >  2:680447/1‑76 (MQ=255)
                     aaCCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATaaa                                                      >  2:299316/1‑76 (MQ=255)
                     aaCCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATaaa                                                      >  2:1252879/1‑76 (MQ=255)
                        cTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCa                                                   >  2:866046/1‑76 (MQ=255)
                        cTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCa                                                   >  2:1512012/1‑76 (MQ=255)
                         tGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAg                                                  <  1:744683/76‑1 (MQ=255)
                         tGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAg                                                  <  2:1196208/76‑1 (MQ=255)
                         tGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAg                                                  <  2:220944/76‑1 (MQ=255)
                          ggCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAgg                                                 <  2:372768/76‑1 (MQ=255)
                             cAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATg                                              <  2:938651/76‑1 (MQ=255)
                             cAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATg                                              <  2:1622891/76‑1 (MQ=255)
                             cAGAAAGGCCAGAGATAAAACGTGGAAATTATGTGGTTTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATg                                              <  2:1573615/76‑1 (MQ=38)
                              aGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGt                                             <  2:1580656/76‑1 (MQ=255)
                                 aaGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAat                                          <  1:292861/76‑1 (MQ=255)
                                 aaGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAat                                          <  1:1380832/76‑1 (MQ=255)
                                 aaGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAat                                          >  1:1101750/1‑76 (MQ=255)
                                   ggCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAattat                                       <  1:339847/76‑1 (MQ=255)
                                      cAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCGCCTAATAAATCAGGATGTAattaaac                                     <  2:80337/76‑4 (MQ=255)
                                      cAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTAttt                                     >  2:502635/1‑76 (MQ=255)
                                       agagATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTAtttt                                    >  1:111949/1‑76 (MQ=255)
                                       agagATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTAtttt                                    <  2:279028/76‑1 (MQ=255)
                                       agagATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTAtttt                                    <  2:178774/76‑1 (MQ=255)
                                           ataataCGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTg                                <  1:1447875/76‑1 (MQ=255)
                                             aataCGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGAt                              >  1:1173972/1‑76 (MQ=255)
                                              ataCGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGAtt                             >  2:1290753/1‑76 (MQ=255)
                                              ataCGTTGAAATTATGTGGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGAttt                            <  1:616727/76‑1 (MQ=255)
                                                 cGTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTa                          <  2:383412/76‑1 (MQ=255)
                                                  gTTGAAATTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTAt                         <  2:77683/76‑1 (MQ=255)
                                                      aaaTTATGTGGTGTATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGtt                     <  2:1870731/76‑1 (MQ=255)
                                                         ttATGTGGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGttcttc                 >  1:1623344/1‑76 (MQ=255)
                                                            tgtgGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTcc              <  1:1003933/76‑1 (MQ=255)
                                                                gtgtATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGa           <  1:123477/76‑1 (MQ=255)
                                                                  gtATTATTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGaaa         >  2:452025/1‑76 (MQ=255)
                                                                   tattatTAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAg        >  2:1330253/1‑76 (MQ=255)
                                                                    attattAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGt       >  2:986558/1‑76 (MQ=255)
                                                                    attattAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGt       >  1:1589796/1‑76 (MQ=255)
                                                                    attattAGGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGt       >  2:45417/1‑76 (MQ=255)
                                                                      tattagGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGTga     >  2:381708/1‑76 (MQ=255)
                                                                        ttaGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGTgataa  <  1:478296/76‑1 (MQ=255)
                                                                           |                                                                         
TTTGTCCTTAATTCCAGAAATAACCTGGCCAGAAAGGCCAGAGATAATACGTTGAAATTATGTGGTGTATTATTAGGGGGGTACCCCCCTAATAAATCAGGATGTAATTATTTTTCTGATTTTATCGTTCTTCTCCAGAAAGTGATAA  >  REL606/795241‑795388

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