Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	2,103,915	(CCAG)7→6	coding (182‑185/216 nt)	ECB_01992 →	hypothetical protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	2103888	2103891	4	61 [0]	[0] 61	ECB_01992	hypothetical protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2103892 =	0 (0.000)	37 (0.990)	30/96	0.2	100%	coding (159/216 nt)	ECB_01992	hypothetical protein
	?	REL606	= 2103887	0 (0.000)					coding (154/216 nt)	ECB_01992	hypothetical protein

ATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTA  >  REL606/2103817‑2103966                                                                            |                                                                           aaacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcac                                                                            >  2:813683‑M1/76‑76 (MQ=255) aaacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcac                                                                            <  1:729480‑M1/1‑1 (MQ=255)  aacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaca                                                                           >  1:1164751‑M1/75‑76 (MQ=255)  aacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaca                                                                           >  2:203456‑M1/75‑76 (MQ=255)  aacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaca                                                                           >  2:1296504‑M1/75‑76 (MQ=255)  aacatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaca                                                                           >  1:397160‑M1/75‑76 (MQ=255)   acatggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacag                                                                          <  2:115337‑M1/3‑1 (MQ=255)     atggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagcn                                                                        >  1:695115‑M1/72‑75 (MQ=255)     atggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagcc                                                                        >  1:1179355‑M1/72‑76 (MQ=255)      tggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagcca                                                                       >  2:1203379‑M1/71‑76 (MQ=255)       ggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccag                                                                      <  1:1266758‑M1/7‑1 (MQ=255)       ggctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccag                                                                      <  2:1591265‑M1/7‑1 (MQ=255)         ctatacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagcc                                                                    <  2:202024‑M1/9‑1 (MQ=255)           atacaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccag                                                                  <  1:1459611‑M1/11‑1 (MQ=255)             acaatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagcc                                                                >  2:862011‑M1/64‑76 (MQ=255)              caatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagcca                                                               <  1:581468‑M1/14‑1 (MQ=255)               aatgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccag                                                              <  1:759290‑M1/15‑1 (MQ=255)                 tgaaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccagcc                                                            >  1:972475‑M1/60‑76 (MQ=255)                   aaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccagccag                                                          >  2:1147008‑M1/58‑76 (MQ=255)                   aaaactgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccagccag                                                          <  2:358137‑M1/19‑1 (MQ=255)                       ctgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccagccagccag                                                      >  2:960170‑M1/54‑76 (MQ=255)                       ctgcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcacagccagccagccagccagccag                                                      <  1:675541‑M1/23‑1 (MQ=255)                         gcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaCagccagccagccagccagccagcc                                                    <  2:441672‑M1/25‑1 (MQ=255)                         gcttcgttattgacttcgaatataaaggggtgaaatttgcattagatttcaCagccagccagccagccagccagcc                                                    >  2:1774373‑M1/52‑76 (MQ=255)                            tcgttattgacttcgaatataagggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAgg                                                 >  2:1179340‑M1/49‑76 (MQ=255)                             cgttattgacttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGgt                                                >  2:1643075‑M1/48‑76 (MQ=255)                                tattgacttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGtgtg                                             >  1:1246266‑M1/45‑76 (MQ=255)                                 attgacttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGa                                            <  1:541420‑M1/33‑1 (MQ=255)                                   tgacttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAAt                                          <  2:1735886‑M1/35‑1 (MQ=255)                                   tgacttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAAt                                          <  1:34090‑M1/35‑1 (MQ=255)                                     acttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAAtgt                                        >  2:1579497‑M1/40‑76 (MQ=255)                                     acttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAAtgt                                        >  1:578687‑M1/40‑76 (MQ=255)                                       ttcgaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTgg                                      >  1:1641896‑M1/38‑76 (MQ=255)                                          gaatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGAtt                                   <  1:1658750‑M1/42‑1 (MQ=255)                                           aatataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCAAGCCAGCCAGGTGTGAATGTGGAttt                                  <  1:1461664‑M1/43‑1 (MQ=255)                                            atataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGAtttt                                 >  2:848216‑M1/33‑76 (MQ=255)                                              ataaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGt                               <  2:891442‑M1/46‑1 (MQ=255)                                               taaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGta                              >  1:697332‑M1/30‑76 (MQ=255)                                                aaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtac                             >  1:1771928‑M1/29‑76 (MQ=255)                                                aaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtac                             >  1:1795300‑M1/29‑76 (MQ=255)                                                aaaggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtac                             <  1:615865‑M1/48‑1 (MQ=255)                                                  aggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCTAGCCAGGTGTGAATGTGGATTTTGtacta                           >  2:1343860‑M1/27‑76 (MQ=255)                                                  aggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtacta                           <  1:645584‑M1/50‑1 (MQ=255)                                                   ggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtactac                          >  1:979556‑M1/26‑76 (MQ=255)                                                   ggggtgaaatttgcattagatttcaCAGCCAGCCAGCCAACCAGCCAGCCAGGTGTGAATGTGGATTTTGtactac                          >  2:751104‑M1/26‑76 (MQ=255)                                                        gaaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAaat                     >  2:1774796‑M1/21‑76 (MQ=255)                                                         aaatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAaata                    >  1:198112‑M1/20‑76 (MQ=255)                                                          aatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAaataa                   <  1:1303033‑M1/58‑1 (MQ=255)                                                          aatttgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAaataa                   <  1:653163‑M1/58‑1 (MQ=255)                                                              tgcattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGc               <  1:1027026‑M1/62‑1 (MQ=255)                                                                cattagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGccc             <  2:1630160‑M1/64‑1 (MQ=255)                                                                  ttagatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCtt           >  1:882389‑M1/11‑76 (MQ=255)                                                                    agatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCttat         >  2:1652485‑M1/9‑76 (MQ=255)                                                                    agatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCttat         <  2:1580313‑M1/68‑1 (MQ=255)                                                                    agatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTATGTCGTACGAAATAACTGCCCttat         <  2:976428‑M1/68‑1 (MQ=255)                                                                     gatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCttatt        <  2:373242‑M1/69‑1 (MQ=255)                                                                     gatttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCttatt        <  2:17996‑M1/69‑1 (MQ=255)                                                                       tttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTtattat      >  1:501781‑M1/6‑76 (MQ=255)                                                                       tttcaCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTtattat      <  1:1424084‑M1/71‑1 (MQ=255)                                                                          caCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACt   >  1:1420135‑M1/3‑76 (MQ=255)                                                                           aCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACta  >  1:1771870‑M1/2‑76 (MQ=255)                                                                            |                                                                           ATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTA  >  REL606/2103817‑2103966

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