Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,185,131	Δ4,735 bp		ygiL–[yqiG]	ygiL, yqiG, insC5, insD5, [yqiG]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3185131	3189865	4735	42 [0]	[1] 44	ygiL–[yqiG]	ygiL,yqiG,insC5,insD5,[yqiG]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3185130	0 (0.000)	42 (0.650)	40/472	0.2	100%	coding (291/291 nt)	ubiK	ubiquinone biosynthesis accessory factor UbiK
	?	NC_000913	3189866 =	0 (0.000)					intergenic (+1/‑15)	yqiG/yqiH	putative outer membrane usher protein YqiG/putative fimbrial chaperone YqiH

ATTCGCCTGGCAGATGGTAGTCATCCACCACTGGGGATTAGCGTAAAAGATGAAACCAGCCACAAAGAATTAGGACTGGTTGCTGATGGCGGCTTTGTATACCTCAACGGCATTCAGGATGATAACAAACTTGCTTTACGCTGGGGTGACAAATCTTGTTTTATTCAACCACCCAATAGCAGCAACTTAACCACCGGAACGGCTATTTTACCGTGTATTAGCCAAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCC  >  NC_000913/3189635‑3190116                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           ggaatcagggagttcggggaagatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATgcg                                                                                                                                                                                                                                         <  2:149736‑M2/20‑1 (MQ=255) ggaatcagggagttcggggaagatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATgcg                                                                                                                                                                                                                                         >  2:389951‑M2/232‑251 (MQ=255)            gttcggggaagatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagctcacaaagcagccagatccagcgactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAAt                                                                                                                                                                                                                               >  2:419998‑M2/221‑250 (MQ=255)               cggggaagatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaactcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtcaTTCAGGCTGCAACTAATGCGCTACTTAAATACa                                                                                                                                                                                                                            >  2:56446‑M2/218‑250 (MQ=255)                 ggtaagatgtggagaaaaaaacccgccaaaccctacaagcgcatctgacgcgcctcgaccttgtaagccgtgaatagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaagcagcgaacgggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATCCAAAg                                                                                                                                                                                                                         <  2:30068‑M2/36‑1 (MQ=255)                 gggaagatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAg                                                                                                                                                                                                                         <  1:230993‑M2/36‑1 (MQ=255)                      gatgtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAatat                                                                                                                                                                                                                    >  1:15544‑M2/211‑251 (MQ=255)                         gtggagaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCAt                                                                                                                                                                                                                 <  1:56446‑M2/44‑1 (MQ=255)                              gaaaaaaatccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTgc                                                                                                                                                                                                              >  1:485503‑M2/203‑249 (MQ=255)                                  aaaagccgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAgg                                                                                                                                                                                                        <  2:484794‑M2/53‑1 (MQ=255)                                   aaatcagccaaacccaaaaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtacaaaagcaagtcctgttacggactcgtgacaaactggcaatgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAAAAAATGCGCTACTTAAATACAAAGAAtatcattgaagaagg                                                                                                                                                                                                        >  2:460305‑M2/198‑243 (MQ=255)                                        cgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAgg                                                                                                                                                                                                                                                       >  1:450041‑M2/193‑198 (MQ=255)                                        cgccaaaccctacaagcgcagctgacgcgcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAgg                                                                                                                                                                                                                                                       <  2:450041‑M2/6‑1 (MQ=255)                                                  cacaagcgcagctgacgcgcctcgaccttgtnagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgctttaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGATTATTACTCTCCTGt                                                                                                                                                                                        <  2:233361‑M2/69‑1 (MQ=255)                                                                    gcctcgaccttgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCCCCTGTATGAGCAGCATTGCATgg                                                                                                                                                                      >  1:242989‑M2/165‑251 (MQ=255)                                                                              tgtaagccgtgaagagttcgacgtccaaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCcatc                                                                                                                                                             >  2:65549‑M2/155‑250 (MQ=255)                                                                                        gaagagttcgacgtccaaacgcaagtcctggtacggactcgtgaaaaacgggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATc                                                                                                                                                   <  2:389936‑M2/106‑1 (MQ=255)                                                                                                        aaacgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaacgaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTACTGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGTGTGCCATCCTCCCGGATCGTACCCGTATCATaatg                                                                                                                                  >  2:378355‑M2/129‑251 (MQ=255)                                                                                                          acgcaagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACGAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACGCTCCTGTATGAGCAGCATTGCATGGGGTGCCATCATCCCGGATCTTACCCGTATCATaatgaa                                                                                                                                <  2:419656‑M2/125‑1 (MQ=255)                                                                                                              aagtcctgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAg                                                                                                                            <  2:153984‑M2/129‑1 (MQ=255)                                                                                                                    tgttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGa                                                                                                                       >  1:274831‑M2/117‑250 (MQ=255)                                                                                                                     gttacggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAgac                                                                                                                      <  1:68756‑M2/135‑2 (MQ=255)                                                                                                                          ggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCAt                                                                                                                                                                             >  2:233344‑M2/111‑190 (MQ=255)                                                                                                                          ggactcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCAt                                                                                                                                                                             <  1:233344‑M2/80‑1 (MQ=255)                                                                                                                          ggactcgtgaaaaactggcgctgcttgaacaaagaagcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctctccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATAGCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCAGCCCGGATCGTACCCGTACCATAATGAATTAGTCAGACAAAGGc                                                                                                                <  2:49015‑M2/141‑1 (MQ=255)                                                                                                                            actcgtgaaaaactggcgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGa                                                                                                              <  1:427112‑M2/143‑1 (MQ=255)                                                                                                                                            cgctgcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACt                                                                                              >  1:192153‑M2/93‑251 (MQ=255)                                                                                                                                                gcttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATc                                                                                          <  2:15544‑M2/163‑1 (MQ=255)                                                                                                                                                  ttgaacaaagaatcagcgaactggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                         >  1:335093‑M2/87‑250 (MQ=255)                                                                                                                                                                        ggaaaatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTACAAACCTGCCCTATCTgg                                                                  >  2:129187‑M2/65‑251 (MQ=255)                                                                                                                                                                            aatcgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCa                                                              <  1:389951‑M2/191‑1 (MQ=255)                                                                                                                                                                               cgcagtactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAAc                                                           <  2:335093‑M2/194‑1 (MQ=255)                                                                                                                                                                                 cagtactgagatcaaaaagcagccagccccagagactctccccccaacgctgtaaTTCAGGACGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAACGCGAAGCATTAAAACTTACTAATAAAAGTAAAAACCTGCCCTATCTGGACCAATCAt                                                         >  2:51536‑M2/56‑251 (MQ=255)                                                                                                                                                                                    tactgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATgg                                                       >  2:179862‑M2/53‑250 (MQ=255)                                                                                                                                                                                     actgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTAt                                                                                                                                                                                                   <  1:5647‑M2/58‑1 (MQ=255)                                                                                                                                                                                     actgagatcaaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTAt                                                                                                                                                                                                   >  2:5647‑M2/52‑109 (MQ=255)                                                                                                                                                                                              aaaaagcagccagatccagagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACCAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATa                                             >  2:64105‑M2/43‑250 (MQ=255)                                                                                                                                                                                                aaagcagccagatccagagactctcccaccaacgctgtatTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAGCTTACTAATCAAAGTAAAAAACTGCCCAATCTGGCCCAAACATGGATAGAAGATACg                                           >  2:327061‑M2/41‑250 (MQ=255)                                                                                                                                                                                                              cagagactctcccccgaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAAACTGCCCTATCTGGCCCAAACTTGGATAGAAGATACGAAAGGCAATAAATcg                            >  2:449098‑M2/27‑251 (MQ=255)                                                                                                                                                                                                               agagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTa                                                                       >  2:216374‑M2/26‑207 (MQ=255)                                                                                                                                                                                                                gagactctccccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTAGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCAGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGAAATAAATcgcg                          >  2:416699‑M2/25‑251 (MQ=255)                                                                                                                                                                                                                 ggacgatccccccaacgctgtaattCAGGATGAAACTAATGCGCTACTGAAATTCAAAGAAAATCCTTGCAGCAGGAGTATTACTCTCCCGTATGAGCATGAATGCCTGGGGTGCCATCAGCCCGGATCATACCCGTAACATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCATAGTTAAAACCTGCCCTTTCGCGCCCAAACGTGGATAGAAGATACGAAATGCAATAAATCGCGt                         <  2:210160/231‑1 (MQ=255)                                                                                                                                                                                                                         cccccaacgctgtaaTTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTAt                 >  1:452552‑M2/16‑251 (MQ=255)                                                                                                                                                                                                                                        ttCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTtcctcc  >  2:61096/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           ATTCGCCTGGCAGATGGTAGTCATCCACCACTGGGGATTAGCGTAAAAGATGAAACCAGCCACAAAGAATTAGGACTGGTTGCTGATGGCGGCTTTGTATACCTCAACGGCATTCAGGATGATAACAAACTTGCTTTACGCTGGGGTGACAAATCTTGTTTTATTCAACCACCCAATAGCAGCAACTTAACCACCGGAACGGCTATTTTACCGTGTATTAGCCAAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCC  >  NC_000913/3189635‑3190116

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