Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	69 [0]	[0] 70	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	63 (1.090)	52/268	0.1	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979141‑1979421                                                                                                                                           |                                                                                                                                               ccttatacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                            <  1:291233‑M2/5‑1 (MQ=255)  tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacgtaaatgTATAAg                                                                                                                                           >  1:384503‑M2/138‑143 (MQ=255)  tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgAATAAg                                                                                                                                           >  2:536888‑M2/138‑143 (MQ=255)   ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactgcattacttaacctaagtgtctgagc                                                                                                                                          >  1:503765‑M2/137‑137 (MQ=255)       acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacatccatgTATAAGTCATa                                                                                                                                      >  1:740034‑M2/133‑143 (MQ=255)        cttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                     <  1:1048965‑M2/12‑1 (MQ=255)        cttatcgaggattttaccgaaatcgacgcaactgtacccgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcacttacttaacataaatgTATAAGTCATAc                                                                                                                                     <  1:480214‑M2/12‑1 (MQ=255)            tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattgtcaatagtgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                 <  2:132596‑M2/16‑1 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                               >  2:450824‑M2/126‑143 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                               >  1:684987‑M2/126‑143 (MQ=255)               gacactgttcctaaatcgacgcaactgtcctcttcactacacgcacatacaacggatgggcgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                              <  2:584646‑M2/19‑1 (MQ=255)                  atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                           <  1:336452‑M2/22‑1 (MQ=255)                   tttttcctaaatcgagggaacagtacgcgtcacgacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  1:47132‑M2/23‑1 (MQ=255)                     tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                        <  1:1087923‑M2/25‑1 (MQ=255)                          taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatcttcaataatgagtcatgcagatcacaaaaaacaatcaattacttaacataaatgTATAAGACATACTTtggatttggtagaatt                                                                                                                   >  2:74116‑M2/114‑136 (MQ=255)                              tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATACGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                                                               >  1:951405‑M2/110‑143 (MQ=255)                                   gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                          <  1:1090667‑M2/39‑1 (MQ=255)                                    caactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                         >  1:675466‑M2/104‑143 (MQ=255)                                               gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                              <  2:666436‑M2/51‑1 (MQ=255)                                               gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                              <  1:445801‑M2/51‑1 (MQ=255)                                               gtcactacacgcacatacaacggaggggggctgcgatattaaaccatagtttgtgcagatcacacacaatacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                              <  1:465243‑M2/51‑1 (MQ=255)                                                tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                             >  1:483200‑M2/92‑143 (MQ=255)                                                  actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCACTCTGTCACAAAGTTTTTcgc                                                                                           >  2:1018626‑M2/90‑143 (MQ=255)                                                    tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaaactccattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                         >  1:1035245‑M2/88‑143 (MQ=255)                                                      cacgcacatccaacgcaggggggctgtgattttcaataatgcgtgatgcagatcacaccaaacactcaatttcttaacataaatgTATAAGTAATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                       <  2:337643‑M2/58‑1 (MQ=255)                                                       acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                      >  1:1038568‑M2/85‑143 (MQ=255)                                                       acgcacagacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                      <  1:899330‑M2/59‑1 (MQ=255)                                                         gcacatacaacggaggtgggctgcgattttcaataatgcgtgatgcagatcactcaaaacactctattacttaacttaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                    <  2:13131‑M2/61‑1 (MQ=255)                                                         gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                    <  1:285132‑M2/61‑1 (MQ=255)                                                         gcacatacaacggaggggggctgcgattttcaataatgcgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTGGTTTTGGGTGTATTACAATTACATGAAAAGAttttcgctacacga                                                                                    >  2:141994‑M2/83‑139 (MQ=255)                                                           acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                  <  2:704544‑M2/63‑1 (MQ=255)                                                             atgcaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                <  2:236089‑M2/65‑1 (MQ=255)                                                              tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                               >  2:55938‑M2/78‑143 (MQ=255)                                                                 aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacatatatgTATAAGTCATACTTTTGCCTTGCCTCTCACTTTATCTGTTAAAAAGTTTTTCGCTACgctagcaaggca                                                                            >  2:122723‑M2/75‑140 (MQ=255)                                                                   cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTCAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                          >  1:1024886‑M2/73‑143 (MQ=255)                                                                       gggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTATATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                      <  2:86110‑M2/75‑1 (MQ=255)                                                                               gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                              <  1:131530‑M2/83‑1 (MQ=255)                                                                                cgattttcaaaaatgcgtgatgcagattacatcaaactgtttttgtctgaacagaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                             <  2:834774‑M2/84‑1 (MQ=255)                                                                                 gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacatagatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                            >  1:577404‑M2/59‑143 (MQ=255)                                                                                 gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                            >  2:1111139‑M2/59‑143 (MQ=255)                                                                                    tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                         >  1:26585‑M2/56‑143 (MQ=255)                                                                                     ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTTAATCAGTTAAATAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                        >  1:1106019‑M2/55‑143 (MQ=255)                                                                                      tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTGTGTATTGGGTGTATTTCAATCTGTTAAAAAGTTGTTCGCTACGGTGGCAAGCAAAACTGCAACAGGAGTTATCGa                                                       >  2:994089‑M2/54‑143 (MQ=255)                                                                                        aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATACTCGaaa                                                     >  2:63887‑M2/52‑143 (MQ=255)                                                                                          taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg                                                   <  2:860977‑M2/94‑1 (MQ=255)                                                                                           aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAGATgg                                                  >  1:850637‑M2/49‑143 (MQ=255)                                                                                              gcgtgatgcagataacacaaaacactgtattacttaacataaatgTATAAGTCATACGTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGTTTTTCACTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                               <  2:503765‑M2/98‑1 (MQ=255)                                                                                                gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACTTTTGGTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAAAGGAAACAGGAATATTCGAAATGGGATGt                                             >  2:32699‑M2/44‑143 (MQ=255)                                                                                                 tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                            <  2:316102‑M2/101‑1 (MQ=255)                                                                                                 tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAACAATGAAACAGGAATAATCGAAATGGGATGtt                                            <  2:675466‑M2/101‑1 (MQ=255)                                                                                                     gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcgc                                        >  1:1010124‑M2/39‑143 (MQ=255)                                                                                                     gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAACCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAACTGGCCTGTCgcgc                                        >  1:486716‑M2/39‑143 (MQ=255)                                                                                                      cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                       >  2:313051‑M2/38‑143 (MQ=255)                                                                                                      cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGCTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGGAACAGGAATAATCGAAATGGGATGTTGCGca                                       >  2:805556‑M2/38‑143 (MQ=255)                                                                                                         agctcacaaaacgcgcaattacttaacttaagtgTATAAGTCATACTTCTGTGGTGGGTGTATGTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                    <  2:840767‑M2/109‑1 (MQ=255)                                                                                                                 aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGTTAGCAAGCAAAATTGAAACAGGAATAATCGAAATGGGATGTTGAGCACAGtcaaaaaa                            >  2:86718‑M2/27‑141 (MQ=255)                                                                                                                  aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                           <  2:483200‑M2/118‑1 (MQ=255)                                                                                                                      ctcaattacttaacataagtgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                       <  2:431645‑M2/122‑1 (MQ=255)                                                                                                                       tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTCTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                      <  1:1020717‑M2/123‑1 (MQ=255)                                                                                                                         aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGATTAATCGAAATGGGTTGTTGCGCACAGTCAAAATAACTCACCg                    >  1:473637‑M2/19‑143 (MQ=255)                                                                                                                               ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata              <  1:782077‑M2/131‑1 (MQ=255)                                                                                                                                  acatatatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           <  2:858086‑M2/134‑1 (MQ=255)                                                                                                                                  acataaatgTATGAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           <  1:74116‑M2/134‑1 (MQ=255)                                                                                                                                   cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatca          <  1:1094040‑M2/135‑1 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCGCACCGTCTTATTAACTCACCGTATAtaatcct         >  1:691859‑M2/8‑141 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat         <  2:674936‑M2/136‑1 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat         >  1:698390‑M2/8‑143 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAACTGTTGAGCTAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTAAAAATAACTCACTCTAATTAatcat         >  2:840905‑M2/8‑143 (MQ=255)                                                                                                                                      aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt       <  2:577404‑M2/138‑1 (MQ=255)                                                                                                                                           tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  >  1:756056/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                               GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979141‑1979421

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