Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	147,727	Δ12 bp	coding (760‑771/903 nt)	yadD →	transposase_31 family protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	147727	147738	12	79 [0]	[0] 78	yadD	transposase_31 family protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 147726	0 (0.000)	69 (0.830)	48/186	0.9	100%	coding (759/903 nt)	yadD	transposase_31 family protein
	?	NC_000913	147739 =	0 (0.000)					coding (772/903 nt)	yadD	transposase_31 family protein

TACTGAACAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGC  >  NC_000913/147627‑147824                                                                                                    |                                                                                                   tACTGAACAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATT                                                                                                    >  1:297670‑M1/1‑100 (MQ=255)  aCTGAACAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTc                                                                                                   >  1:1533684‑M1/1‑99 (MQ=255)    tGAACAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcag                                                                                                 <  2:28195‑M1/100‑4 (MQ=255)      aaCAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagca                                                                                               <  2:622960‑M1/100‑6 (MQ=255)         aaGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcaggg                                                                                            <  2:1264680‑M1/100‑9 (MQ=255)              gATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagac                                                                                       <  2:158936‑M1/100‑14 (MQ=255)              gATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagac                                                                                       <  2:963179‑M1/100‑14 (MQ=255)                 ttGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacagg                                                                                    <  2:92933‑M1/100‑17 (MQ=255)                 ttGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacagg                                                                                    <  2:1645276‑M1/100‑17 (MQ=255)                   gTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaa                                                                                  >  1:789509‑M1/1‑82 (MQ=255)                   gTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaa                                                                                  <  1:510161‑M1/100‑19 (MQ=255)                       tACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaa                                                                              >  2:588459‑M1/1‑78 (MQ=255)                        aCGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaag                                                                             <  1:1575165‑M1/100‑24 (MQ=255)                             gtgtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcagg                                                                        <  2:778717‑M1/100‑29 (MQ=255)                             gtgtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcagg                                                                        >  2:1180643‑M1/1‑72 (MQ=255)                             gtgtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcagg                                                                        <  1:1809697‑M1/100‑29 (MQ=255)                               gtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaa                                                                      >  2:972529‑M1/1‑70 (MQ=255)                               gtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaa                                                                      >  1:1458652‑M1/1‑70 (MQ=255)                               gtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaa                                                                      >  1:1269077‑M1/1‑70 (MQ=255)                               gtTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcagga                                                                       >  2:1964378‑M1/1‑70 (MQ=255)                                 tGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaatt                                                                    >  2:1898028‑M1/1‑68 (MQ=255)                                 tGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaatt                                                                    >  1:348512‑M1/1‑68 (MQ=255)                                    gagaCAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgc                                                                 >  2:256424‑M1/1‑65 (MQ=255)                                    gagaCAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgc                                                                 >  2:1955507‑M1/1‑65 (MQ=255)                                    gagaCAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgc                                                                 >  1:506294‑M1/1‑65 (MQ=255)                                      gaCAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTGTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgccc                                                               >  2:362046‑M1/1‑63 (MQ=255)                                       aCAGGGAAACGGGAGGGGAGTCTATGAGGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgccca                                                              <  1:2208694‑M1/100‑39 (MQ=255)                                        cAGGGAGACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggacttcgcccag                                                             <  2:588687‑M1/100‑40 (MQ=255)                                        cAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccag                                                             >  1:1844052‑M1/1‑61 (MQ=255)                                        cAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccag                                                             >  1:1985686‑M1/1‑61 (MQ=255)                                         aGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagc                                                            <  2:1178360‑M1/100‑41 (MQ=255)                                          gggAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcg                                                           >  2:1731034‑M1/1‑59 (MQ=255)                                          gggAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcg                                                           <  2:1844040‑M1/100‑42 (MQ=255)                                          gggAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcg                                                           >  1:632213‑M1/1‑59 (MQ=255)                                                cGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttct                                                     <  1:1582624‑M1/100‑48 (MQ=255)                                                cGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttct                                                     >  2:2242004‑M1/1‑53 (MQ=255)                                                 gggaggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctg                                                    <  2:33052‑M1/100‑49 (MQ=255)                                                 gggaggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctg                                                    <  2:1146248‑M1/100‑49 (MQ=255)                                                 gggaggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTaagcaaggaagacaggaagtaagtcaggaattcgcccagcgtcttctg                                                    <  2:104617‑M1/100‑49 (MQ=255)                                                  ggaggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctga                                                   >  2:297175‑M1/1‑51 (MQ=255)                                                    aggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagt                                                 <  1:633370‑M1/100‑52 (MQ=255)                                                     ggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagta                                                >  1:1623569‑M1/1‑48 (MQ=255)                                                     ggggagTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagta                                                >  1:1938000‑M1/1‑48 (MQ=255)                                                          ggctaTGATGACGCGGGCGAAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaaggcaggaattcgcccagcgtcttctgagtaaagg                                            <  2:2153751‑M1/97‑57 (MQ=255)                                                          gTCTATGATGACGCTGGCGCCGTGGTTTGCAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggcattcgcccagcgtcttctgagtaaagga                                           >  2:1425495‑M1/1‑43 (MQ=255)                                                          gTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagaaaagga                                           >  2:1915633‑M1/1‑43 (MQ=255)                                                            cTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaat                                         >  2:710757‑M1/1‑41 (MQ=255)                                                               tgatggCGGTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtctgctgaggaaaggaatgtc                                      >  2:329392‑M1/1‑38 (MQ=255)                                                                gatgaCGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtct                                     >  1:173971‑M1/1‑37 (MQ=255)                                                                gatgaCGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtct                                     >  1:115030‑M1/1‑37 (MQ=255)                                                                gatgaCGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtct                                     >  1:2233365‑M1/1‑37 (MQ=255)                                                                          gcgcAGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaaggcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacg                           <  2:479955‑M1/100‑74 (MQ=255)                                                                              aGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgc                       <  1:256427‑M1/100‑78 (MQ=255)                                                                              aGTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgc                       <  1:1955521‑M1/100‑78 (MQ=255)                                                                               gTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgca                      <  2:716866‑M1/100‑79 (MQ=255)                                                                               gTGGTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgca                      >  1:1189755‑M1/1‑22 (MQ=255)                                                                                 ggTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcaga                    <  1:1898040‑M1/100‑81 (MQ=255)                                                                                  gTTTGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagag                   >  1:2027438‑M1/1‑19 (MQ=255)                                                                                   tttGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagaga                  >  1:1868004‑M1/1‑18 (MQ=255)                                                                                    ttGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagat                 <  2:1565431‑M1/100‑84 (MQ=255)                                                                                    ttGAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagat                 <  2:297667‑M1/100‑84 (MQ=255)                                                                                      gAAGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatg                <  2:1813105‑M1/99‑85 (MQ=255)                                                                                       aaGAGAAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggc              <  1:972537‑M1/100‑87 (MQ=255)                                                                                        agagAAAGGGATTcagcagggaagacaggaagtaagtaaggaattcgcccagcgtcttctgagtaaagaaatgtcccgggaagacgttgcagagatggca             >  2:709820‑M1/1‑13 (MQ=255)                                                                                         gagaAAGGGATTcagcagggaagacagaaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaa            >  1:1730134‑M1/1‑12 (MQ=255)                                                                                           gaAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaat          >  1:1100046‑M1/1‑10 (MQ=255)                                                                                           gaAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaat          <  1:329395‑M1/100‑91 (MQ=255)                                                                                           gaAAGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaat          >  1:1951016‑M1/1‑10 (MQ=255)                                                                                             aaGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaattt        >  2:570801‑M1/1‑8 (MQ=255)                                                                                             aaGGGATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaattt        >  2:807118‑M1/1‑8 (MQ=255)                                                                                               gggATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttac      >  1:1625681‑M1/1‑6 (MQ=255)                                                                                               gggATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttac      <  2:1269070‑M1/100‑95 (MQ=255)                                                                                                tctttcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttacc     <  1:1333172‑M1/97‑96 (MQ=255)                                                                                                ggATTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttacc     <  2:1697835‑M1/100‑96 (MQ=255)                                                                                                ggATTcagcagggaagaaaggaaataagtcaggaatacgcacagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaattta       <  2:1233661‑M1/98‑94 (MQ=255)                                                                                                  aTTcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttacct    >  2:560680‑M1/1‑3 (MQ=255)                                                                                                   ttcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttacctct  >  2:40307‑M1/1‑2 (MQ=255)                                                                                                   ttcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttacctct  >  2:1613157‑M1/1‑2 (MQ=255)                                                                                                   ttcagcagggaagacaggaagtaagtcaggaattcgcccagcgtcttctgagtaaaggaatgtctcgggaagacgttgcagagatggcaaatttaccaga  >  2:1333163‑M1/1‑2 (MQ=255)                                                                                                    |                                                                                                   TACTGAACAAGCGGATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGC  >  NC_000913/147627‑147824

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