Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,305,098	Δ13 bp	coding (1488‑1500/1986 nt)	yjcS ←	alkyl sulfatase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4305098	4305110	13	68 [0]	[0] 68	yjcS	alkyl sulfatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4305097	0 (0.000)	66 (0.920)	43/274	0.4	100%	coding (1501/1986 nt)	yjcS	alkyl sulfatase
	?	NC_000913	4305111 =	0 (0.000)					coding (1487/1986 nt)	yjcS	alkyl sulfatase

CGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGC  >  NC_000913/4304961‑4305236                                                                                                                                         |                                                                                                                                            cgAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcg                                                                                                                                          >  1:380893‑M1/1‑137 (MQ=255)      ggcgtccgGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctg                                                                                                                                     <  1:390392‑M1/137‑9 (MQ=255)       gtgtCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctgg                                                                                                                                    <  2:1146788‑M1/140‑10 (MQ=255)           ccGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcac                                                                                                                                <  1:855447‑M1/140‑14 (MQ=255)                 aaCCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggt                                                                                                                          >  2:510185‑M1/1‑121 (MQ=255)                      gtggtgCCGTGGCTGAACTTATGCACCCATTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcg                                                                                                                     >  1:397768‑M1/1‑116 (MQ=255)                             cGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatca                                                                                                              <  1:674069‑M1/140‑32 (MQ=255)                             cGTGGCTGAACTTATGCACCCATTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatca                                                                                                              <  2:397768‑M1/140‑32 (MQ=255)                              gTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcac                                                                                                             >  2:809059‑M1/1‑108 (MQ=255)                               tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacc                                                                                                            <  1:272908‑M1/140‑34 (MQ=255)                               tGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacc                                                                                                            <  1:633229‑M1/140‑34 (MQ=255)                                 gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctg                                                                                                          >  1:207918‑M1/1‑105 (MQ=255)                                 gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatca                                                                                                              >  1:882415‑M1/1‑105 (MQ=255)                                 gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatca                                                                                                              <  2:882415‑M1/136‑32 (MQ=255)                                      aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttca                                                                                                     >  1:483945‑M1/1‑100 (MQ=255)                                               ccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctg                                                                                            >  1:1102198‑M1/1‑91 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  2:436835‑M1/1‑88 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  2:214844‑M1/1‑88 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  2:1124045‑M1/1‑88 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  2:767772‑M1/1‑88 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  1:280909‑M1/1‑88 (MQ=255)                                                  cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccg                                                                                         >  1:1127375‑M1/1‑88 (MQ=255)                                                    tCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgac                                                                                       <  2:38362‑M1/140‑55 (MQ=255)                                                     cgcgCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgacc                                                                                      >  2:345810‑M1/1‑85 (MQ=255)                                                     cgcgCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgacc                                                                                      >  1:929232‑M1/1‑85 (MQ=255)                                                      gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgctacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgacca                                                                                     <  1:37096‑M1/140‑57 (MQ=255)                                                      gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgacca                                                                                     <  1:15663‑M1/140‑57 (MQ=255)                                                       cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccag                                                                                    <  2:1139995‑M1/140‑58 (MQ=255)                                                       cgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccag                                                                                    <  2:964248‑M1/140‑58 (MQ=255)                                                            ctctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggta                                                                               <  2:207918‑M1/140‑63 (MQ=255)                                                                    cgcCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcacctt                                                                       >  1:933797‑M1/1‑70 (MQ=255)                                                                                 gAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttc                                                          <  1:828753‑M1/140‑84 (MQ=255)                                                                                 gAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttc                                                          <  2:483945‑M1/140‑84 (MQ=255)                                                                                  aaaCCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttct                                                         <  1:1124045‑M1/140‑85 (MQ=255)                                                                                   aaCCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttctt                                                        <  2:409318‑M1/140‑86 (MQ=255)                                                                                    aCCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttg                                                       >  2:618431‑M1/1‑54 (MQ=255)                                                                                     ccGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgc                                                      <  2:846447‑M1/140‑88 (MQ=255)                                                                                      cgcgCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcg                                                     <  1:932146‑M1/140‑89 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgcc                                                   >  2:279479‑M1/1‑50 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgcc                                                   >  1:745835‑M1/1‑50 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgcc                                                   >  1:411413‑M1/1‑50 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgcc                                                   >  2:642353‑M1/1‑50 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgcc                                                   <  2:411413‑M1/140‑91 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcg                                                     >  2:538422‑M1/1‑50 (MQ=255)                                                                                        cgcCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcg                                                     <  1:538422‑M1/138‑89 (MQ=255)                                                                                                ggCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgat                                           >  1:309347‑M1/1‑42 (MQ=255)                                                                                                   ggACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgac                                        <  1:809059‑M1/140‑102 (MQ=255)                                                                                                   ggACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgac                                        <  2:396231‑M1/140‑102 (MQ=255)                                                                                                    gACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgaca                                       <  1:226925‑M1/140‑103 (MQ=255)                                                                                                         ggCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggc                                  >  1:156627‑M1/1‑33 (MQ=255)                                                                                                         ggCCTGATAGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggc                                  >  2:852484‑M1/1‑33 (MQ=255)                                                                                                                 aGCCCAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccggttgatgacacgggcagaaccgcc                          <  1:279479‑M1/140‑116 (MQ=255)                                                                                                                   cccAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgccc                        >  1:1106067‑M1/1‑23 (MQ=255)                                                                                                                     cAGCTGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccag                      >  2:740517‑M1/1‑21 (MQ=255)                                                                                                                         tGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcc                  >  1:305472‑M1/1‑17 (MQ=255)                                                                                                                         tGTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcc                  >  1:914052‑M1/1‑17 (MQ=255)                                                                                                                          gTTCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcct                 <  2:380893‑M1/140‑125 (MQ=255)                                                                                                                           ttCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctg                >  1:304048‑M1/1‑15 (MQ=255)                                                                                                                            tCAAAGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgc               >  2:341632‑M1/1‑14 (MQ=255)                                                                                                                               aaGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacg            <  2:743628‑M1/140‑130 (MQ=255)                                                                                                                                aGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgt           <  1:510185‑M1/140‑131 (MQ=255)                                                                                                                                aGTTATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgt           <  2:304048‑M1/140‑131 (MQ=255)                                                                                                                                   tATTCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaac        >  1:720135‑M1/1‑7 (MQ=255)                                                                                                                                     ttCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgcctaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaacgt      >  2:639865‑M1/1‑5 (MQ=255)                                                                                                                                     ttCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaacgt      >  2:209184‑M1/1‑5 (MQ=255)                                                                                                                                      tCGCgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaacgtt     <  1:767772‑M1/140‑137 (MQ=255)                                                                                                                                         cgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaacgtttac  <  1:642353‑M1/140‑140 (MQ=255)                                                                                                                                         cgcgacctggtcacccgggttggcggcaatcacctgtttcagcagttctgccgaccagcggtaatcaccttgcttgttcgcttcttgcgccaggttgatgacacgggcagaaccgcccagcgcctgcacgtaacgtttac  <  2:1102198‑M1/140‑140 (MQ=255)                                                                                                                                         |                                                                                                                                            CGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGC  >  NC_000913/4304961‑4305236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑