BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	401,670	Δ728 bp		[phoA]	[phoA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	401670	402397	728	37 [3]	[7] 41	[phoA]	[phoA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 401669	3 (0.070)	34 (0.810)	32/470	0.3	88.4%	intergenic (+23/‑78)	iraP/phoA	anti‑adaptor protein IraP/alkaline phosphatase
*	?	NC_000913	402398 =	6 (0.140)	34 (0.810)	32/470	0.3	88.4%	coding (652/1416 nt)	phoA	alkaline phosphatase

AAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAAT  >  NC_000913/402153‑402646

caaaaagaagaagagtcgaaagaactgtgtgcgcaggtagaagctttggagattatcgtcactgcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAAc                                                                                                                                                                                                                                                     >  1:437659‑M2/246‑251 (MQ=255)

     agaagaagagtcgaaagaactgtgtgcgcaggtagaagctttggagattatcgtcactgcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCg                                                                                                                                                                                                                                                <  1:284‑M2/11‑1 (MQ=255)

                          gtgtgcgcaggtagaagctttggagattatcgtcactgcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcatcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTggcgg                                                                                                                                                                                                                           >  1:382280‑M2/220‑251 (MQ=255)

                                 caggtagaagctttggagattatcgtcactgcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCa                                                                                                                                                                                                                    <  1:411156‑M2/39‑1 (MQ=255)

                                      agaagctttggagattatcgtcactgcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAAc                                                                                                                                                                                                               <  1:449356‑M2/44‑1 (MQ=255)

                                             ttggagattagcgtcaatgcaatgctacgcaatatggcgcaaaatgaccaacagcggttgattgaacaggtagagggggcgctgtacgaggtaaagcccgatgccagaattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcttcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCGTTGc                                                                                                                                                                                                         <  2:10182‑M2/50‑1 (MQ=255)

                                                         gtcactgcaaagcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAAcc                                                                                                                                                                                            <  2:335222‑M2/63‑1 (MQ=255)

                                                          ggAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAAc                                                                                                                                                                                             <  2:277228/249‑1 (MQ=255)

                                                               gcaatgcttcgcaatatggcgcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGt                                                                                                                                                                                      <  1:265886‑M2/69‑1 (MQ=255)

                                                                  atgcttcgcaatatggagcaaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaaggccgatgccagcattcctgaagacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaagagttaatcttttcaacagatgCTTAACGCTCGTACCGACGATACGCGTGGCTGAGGAGCAAACACCTTTGCTGAACCGGACaccgctcgggc                                                                                                                                                                                    >  2:163650‑M2/180‑245 (MQ=255)

                                                                            atatggggcaaaaggtacaacagctgttgattgactatgtagagggggagctgtaagaggtaaataccgatgccagcattcctgacgacgatacggatctgctgcgctattacgtaaagaagttattgaagcatcctcgtcagtaaaaagctaatcttttcaacagctgCTTAACGTGCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTGGCTGAAACGGCAACCGCTGGTGAATGGCAGGGaa                                                                                                                                                                         <  2:18656‑M2/82‑1 (MQ=255)

                                                                                    caaaatgaccaacagcggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTg                                                                                                                                                                 <  1:49540‑M2/90‑1 (MQ=255)

                                                                                           accaacatcggttgattgatcaggtagagggggcgctgtacgagggagtccccgatgccttcattcctgacgacgttacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgTTTAACGCTCGTGCCGACGTTACTCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCGGGTGAATGGCAGGGAAATACGCGGCGTGAac                                                                                                                                                          <  2:278727‑M2/97‑1 (MQ=255)

                                                                                              cGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAATCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAacagg                                                                                                                                                       >  1:45754/1‑251 (MQ=255)

                                                                                                   cggttgattgatcaggtagagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCCGGGAAAAACGCTGCGTGaacaggcaa                                                                                                                                                    >  2:349368‑M2/147‑248 (MQ=255)

                                                                                                                    agagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGtt                                                                                                                                 <  2:437659‑M2/122‑1 (MQ=255)

                                                                                                                    agagggggcgctgtacgaggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGtt                                                                                                                                 <  1:94584‑M2/122‑1 (MQ=255)

                                                                                                                           cGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAAGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGc                                                                                                                          >  2:201555/1‑251 (MQ=255)

                                                                                                                                    cgctgGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGcc                                                                                                                 >  1:87384/1‑251 (MQ=255)

                                                                                                                                     aggtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCt                                                                                                                <  1:184609‑M2/139‑1 (MQ=255)

                                                                                                                                       gtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgcagcctcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGGGCCGCCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCGGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCa                                                                                                              <  1:313884‑M2/141‑1 (MQ=255)

                                                                                                                                       gtaaagcccgatgccagcattcctgacgacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCa                                                                                                              <  1:463076‑M2/141‑1 (MQ=255)

                                                                                                                                                        ccgcgCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAATCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGa                                                                                             <  2:45754/250‑1 (MQ=255)

                                                                                                                                                                  gacgatacggagctgctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTAACTGAATTCGGTGACGGAAGCGAATc                                                                                     >  1:298761‑M2/84‑249 (MQ=255)

                                                                                                                                                                  gCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATc                                                                                     >  1:124668/1‑249 (MQ=255)

                                                                                                                                                                                gctgcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTg                                                                       <  1:418083‑M2/180‑1 (MQ=255)

                                                                                                                                                                                   gcgcgattacgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTg                                                                   <  1:166798‑M2/184‑1 (MQ=255)

                                                                                                                                                                                           acgtaaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTg                                                          >  1:49413‑M2/59‑251 (MQ=255)

                                                                                                                                                                                              taaagaagttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTg                                                       >  1:185316‑M2/56‑251 (MQ=255)

                                                                                                                                                                                                     gttattgaagcatcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCGGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGAGGCTGCCTCACCGACGTCGGGGACGGAAGCGAATCAGCAAAAACccctgcttggcctggtgtctgtaggca                                                 >  2:357463‑M2/49‑237 (MQ=255)

                                                                                                                                                                                                      ttattgaagcatactcgtcagtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAat                                               >  1:95435‑M2/48‑251 (MQ=255)

                                                                                                                                                                                                       aattgaagcatcctcgacagtaaaaagttagtgttttcaacagctgCTTAACGCTCGTGCCGACGTTACGGTTGGCGGCGGGGAAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAATCGCTGCGTTAACAGGCACAGGCGCGTGGTTATCAGGTGGTGAGCGATGCTGCCTCACTGAATTCGGTGATGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAata                                              <  1:181700‑M2/205‑1 (MQ=255)

                                                                                                                                                                                                        attgaagcatcctcgtcagtacaaagttactcttttcaacagctgCTTAACGCTCGTGCCGACGTTCCGCTTGGCGGCGGCGCAAAAACCTTTGCTGAACAGGCAACCGCTGGTGACTGGCAGGGAAAAACGCTCCGTGAACAGGACCAGACGCGTGGTTATCAGTTGGTGAGCGGTGCTGCCTCACTGAATTCGGGGACGGAAGCGCATCCGCAAAAACCACTGCTTGGCCTGGTTGCTGACGGc                                                  >  2:284872‑M2/46‑246 (MQ=255)

                                                                                                                                                                                                        attgaagcaccctcgtaattaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGCCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAACTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGc                                                  <  1:284872‑M2/201‑1 (MQ=255)

                                                                                                                                                                                                           gaagcatcctagtcagtaaaaagttaatcttttcaacagctgCTTAACGCCCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGc                                           <  1:70923‑M2/208‑1 (MQ=255)

                                                                                                                                                                                                                 tcctcgtcagtaaaaagttaatcttttcaacagctgCTTAACGTCCGGGCCGACTTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgc                                     <  1:184451‑M2/214‑1 (MQ=255)

                                                                                                                                                                                                                   ctcgtgagtaacaagttactcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGAGACTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGCTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCcagtgagcg                                  >  2:194671‑M2/35‑247 (MQ=255)

                                                                                                                                                                                                                         agtaaaaagttaatcttttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAg                            <  1:17927‑M2/223‑1 (MQ=255)

                                                                                                                                                                                                                             aaaaggtaatcttttcaccagctcCTTAACGCTCTTGCCGCCGTTACGCTTCGCGGCGGCGCAAAATCCTTTGCTTTGACCGCACCCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTTTTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGTTGCCGGAAGCGAATCAGCAAAATCCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAcc                        <  1:188673‑M2/227‑1 (MQ=255)

                                                                                                                                                                                                                                     atcatttcaacagctgCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGGGAACAGGCACAGGCGCGTGGTTAGCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAACCCCCTGCTTGGCCTGTTTGCTGACGGCAATATGTCAGTGCGCTGGCTAGGACCGAAAGCaa                <  2:61013‑M2/235‑1 (MQ=255)

                                                                                                                                                                                                                                                    gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGCCGGCCATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGGACCATGGCAat  >  2:58760/1‑250 (MQ=255)

AAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAAT  >  NC_000913/402153‑402646

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