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

GATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCT  >  NC_000913/147640‑147836                                                                                                    |                                                                                                  actgaacaagcggatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattc                                                                                                   >  1:1533684‑M2/100‑100 (MQ=255)   tgaacaagcggatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattcag                                                                                                 <  2:28195‑M2/3‑1 (MQ=255)     aacaagcggatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattcagca                                                                                               <  2:622960‑M2/5‑1 (MQ=255)        aagcggatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAggg                                                                                            <  2:1264680‑M2/8‑1 (MQ=255)             gatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGAc                                                                                       <  2:963179‑M2/13‑1 (MQ=255)             gatttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGAc                                                                                       <  2:158936‑M2/13‑1 (MQ=255)                ttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAgg                                                                                    <  2:1645276‑M2/16‑1 (MQ=255)                ttgttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAgg                                                                                    <  2:92933‑M2/16‑1 (MQ=255)                  gttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGaa                                                                                  >  1:789509‑M2/83‑100 (MQ=255)                  gttttacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGaa                                                                                  <  1:510161‑M2/18‑1 (MQ=255)                      tacggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGaagtaa                                                                              >  2:588459‑M2/79‑100 (MQ=255)                       acggtgtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGaagtaag                                                                             <  1:1575165‑M2/23‑1 (MQ=255)                            gtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAgg                                                                        >  2:1180643‑M2/73‑100 (MQ=255)                            gtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAgg                                                                        <  2:778717‑M2/28‑1 (MQ=255)                            gtgttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAgg                                                                        <  1:1809697‑M2/28‑1 (MQ=255)                              gttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGaa                                                                      >  2:972529‑M2/71‑100 (MQ=255)                              gttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGaa                                                                      >  1:1458652‑M2/71‑100 (MQ=255)                              gttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGaa                                                                      >  1:1269077‑M2/71‑100 (MQ=255)                              gttgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGa                                                                       >  2:1964378‑M2/71‑99 (MQ=255)                                tgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAAtt                                                                    >  2:1898028‑M2/69‑100 (MQ=255)                                tgagagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAAtt                                                                    >  1:348512‑M2/69‑100 (MQ=255)                                   gagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGc                                                                 >  2:256424‑M2/66‑100 (MQ=255)                                   gagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGc                                                                 >  1:506294‑M2/66‑100 (MQ=255)                                   gagacagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGc                                                                 >  2:1955507‑M2/66‑100 (MQ=255)                                     gacagggaaacgggaggggagtctatgatgacgctggcgcagtggtgtgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGccc                                                               >  2:362046‑M2/64‑100 (MQ=255)                                      acagggaaacgggaggggagtctatgaggacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCa                                                              <  1:2208694‑M2/38‑1 (MQ=255)                                       cagggagacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGACTTCGCCCAg                                                             <  2:588687‑M2/39‑1 (MQ=255)                                       cagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAg                                                             >  1:1985686‑M2/62‑100 (MQ=255)                                       cagggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAg                                                             >  1:1844052‑M2/62‑100 (MQ=255)                                        agggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGc                                                            <  2:1178360‑M2/40‑1 (MQ=255)                                         gggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCg                                                           >  1:632213‑M2/60‑100 (MQ=255)                                         gggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCg                                                           <  2:1844040‑M2/41‑1 (MQ=255)                                         gggaaacgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCg                                                           >  2:1731034‑M2/60‑100 (MQ=255)                                               cgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGtcttct                                                     <  1:1582624‑M2/47‑1 (MQ=255)                                               cgggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGtcttct                                                     >  2:2242004‑M2/54‑100 (MQ=255)                                                gggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTg                                                    <  2:1146248‑M2/48‑1 (MQ=255)                                                gggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTg                                                    <  2:33052‑M2/48‑1 (MQ=255)                                                gggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattAAGCAAGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTg                                                    <  2:104617‑M2/48‑1 (MQ=255)                                                 ggaggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGa                                                   >  2:297175‑M2/52‑100 (MQ=255)                                                   aggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGt                                                 <  1:633370‑M2/51‑1 (MQ=255)                                                    ggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTa                                                >  1:1623569‑M2/49‑100 (MQ=255)                                                    ggggagtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTa                                                >  1:1938000‑M2/49‑100 (MQ=255)                                                         gtctatgatgacgctggcgccgtggtttgcagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGCATTCGCCCAGCGTCTTCTGAGTAAAGGa                                           >  2:1425495‑M2/44‑100 (MQ=255)                                                         gtctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGAAAAGGa                                           >  2:1915633‑M2/44‑100 (MQ=255)                                                         ggctatgatgacgcgggcgaagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGGCAGGAATTCGCCCAGCGTCTTCTGAGTAAAgg                                            <  2:2153751‑M2/56‑1 (MQ=255)                                                           ctatgatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAAt                                         >  2:710757‑M2/42‑100 (MQ=255)                                                              tgatggcggtggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTGCTGAGGAAAGGAATGtc                                      >  2:329392‑M2/39‑100 (MQ=255)                                                               gatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGtct                                     >  1:115030‑M2/38‑100 (MQ=255)                                                               gatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGtct                                     >  1:173971‑M2/38‑100 (MQ=255)                                                               gatgacgctggcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGtct                                     >  1:2233365‑M2/38‑100 (MQ=255)                                                                         gcgcagtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGGCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACg                           <  2:479955‑M2/73‑1 (MQ=255)                                                                             agtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGc                       <  1:1955521‑M2/77‑1 (MQ=255)                                                                             agtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGc                       <  1:256427‑M2/77‑1 (MQ=255)                                                                              gtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCa                      <  2:716866‑M2/78‑1 (MQ=255)                                                                              gtggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCa                      >  1:1189755‑M2/23‑100 (MQ=255)                                                                                ggtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCaga                    <  1:1898040‑M2/80‑1 (MQ=255)                                                                                 gtttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCagag                   >  1:2027438‑M2/20‑100 (MQ=255)                                                                                  tttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAgaga                  >  1:1868004‑M2/19‑100 (MQ=255)                                                                                   ttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGAt                 <  2:297667‑M2/83‑1 (MQ=255)                                                                                   ttgaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGAt                 <  2:1565431‑M2/83‑1 (MQ=255)                                                                                     gaagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATg                <  2:1813105‑M2/84‑1 (MQ=255)                                                                                      aagagaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGc              <  1:972537‑M2/86‑1 (MQ=255)                                                                                       agagaaagggattCAGCAGGGAAGACAGGAAGTAAGTAAGGAATTCGCCCAGCGTCTTCTGAGTAAAGAAATGTCCCGGGAAGACGTTGCAGAGATGGCa             >  2:709820‑M2/14‑100 (MQ=255)                                                                                        gagaaagggattCAGCAGGGAAGACAGAAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCaa            >  1:1730134‑M2/13‑100 (MQ=255)                                                                                          gaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAAt          <  1:329395‑M2/90‑1 (MQ=255)                                                                                          gaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAAt          >  1:1951016‑M2/11‑100 (MQ=255)                                                                                          gaaagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAAt          >  1:1100046‑M2/11‑100 (MQ=255)                                                                                            aagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAAttt        >  2:807118‑M2/9‑100 (MQ=255)                                                                                            aagggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAAttt        >  2:570801‑M2/9‑100 (MQ=255)                                                                                              gggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTAc      <  2:1269070‑M2/94‑1 (MQ=255)                                                                                              gggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTAc      >  1:1625681‑M2/7‑100 (MQ=255)                                                                                               tctttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTAcc     <  1:1333172‑M2/95‑1 (MQ=255)                                                                                               ggattCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTAcc     <  2:1697835‑M2/95‑1 (MQ=255)                                                                                               ggattCAGCAGGGAAGAAAGGAAATAAGTCAGGAATACGCACAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTa       <  2:1233661‑M2/93‑1 (MQ=255)                                                                                                 attCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACct    >  2:560680‑M2/4‑99 (MQ=255)                                                                                                  ttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATttaccaga  >  2:1333163‑M2/3‑97 (MQ=255)                                                                                                  ttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACctct  >  2:40307‑M2/3‑100 (MQ=255)                                                                                                  ttCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACctct  >  2:1613157‑M2/3‑100 (MQ=255)                                                                                                    |                                                                                                  GATTTGTTTTACGGTGTGTTGAGAGACAGGGAAACGGGAGGGGAGTCTATGATGACGCTGGCGCAGTGGTTTGAAGAGAAAGGGATTGAGAAGGGGATTCAGCAGGGAAGACAGGAAGTAAGTCAGGAATTCGCCCAGCGTCTTCTGAGTAAAGGAATGTCTCGGGAAGACGTTGCAGAGATGGCAAATTTACCTCT  >  NC_000913/147640‑147836

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: ‑