Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	515,808	C→T	intergenic (‑35/‑111)	qmcA ← / → fetA	PHB domain membrane‑anchored putative protease/iron export ABC transporter ATPase, peroxide resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	515,808	0	C	T	100.0%	239.5 / NA	74	intergenic (‑35/‑111)	qmcA/fetA	PHB domain membrane‑anchored putative protease/iron export ABC transporter ATPase, peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/40);  total (34/40)

CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTT  >  NC_000913/515680‑515942                                                                                                                                 |                                                                                                                                        cTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                               <  2:1198852/139‑1 (MQ=255)       tCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGAt                                                                                                                         <  2:50190/139‑1 (MQ=255)        ccACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGAtt                                                                                                                        >  2:767263/1‑139 (MQ=255)         cACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTa                                                                                                                       >  1:269752/1‑139 (MQ=255)          aCTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                      <  2:885726/139‑1 (MQ=255)          aCTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                      <  2:213721/139‑1 (MQ=255)           cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     <  2:24141/139‑1 (MQ=255)           cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     <  2:130621/139‑1 (MQ=255)           cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     <  1:1801901/139‑1 (MQ=255)            tGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAAc                                                                                                                    >  1:879684/1‑139 (MQ=255)                 gCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     <  2:550300/133‑1 (MQ=255)                 gCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     >  1:550300/1‑133 (MQ=255)                 gCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAAtt                                                                                                               >  2:647701/1‑139 (MQ=255)                  cccTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATta                                                                                                              <  2:305512/139‑1 (MQ=255)                   ccTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATtat                                                                                                             >  1:1304392/1‑139 (MQ=255)                    cTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATtata                                                                                                            <  2:535891/139‑1 (MQ=255)                      gCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATAcc                                                                                                          >  2:1694262/1‑139 (MQ=255)                        ggTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGGTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTg                                                                                                        >  2:1245707/1‑139 (MQ=255)                            cGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTa                                                                                                    >  1:285972/1‑139 (MQ=255)                                    acacCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     <  1:1491779/114‑1 (MQ=255)                                    acacCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                     >  2:1491779/1‑114 (MQ=255)                                    acacCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAaga                                                                                            <  1:589778/139‑1 (MQ=255)                                       cccGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAagagag                                                                                         <  1:286864/139‑1 (MQ=255)                                        ccGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTg                                                                                                 <  1:633785/130‑1 (MQ=255)                                        ccGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTg                                                                                                 >  2:633785/1‑130 (MQ=255)                                           cgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGATAGTTcc                                                                                     >  2:228004/1‑139 (MQ=255)                                           cgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTcc                                                                                     >  2:1652620/1‑139 (MQ=255)                                                    tGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTc                                                                           <  1:1652620/139‑1 (MQ=255)                                                    tGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTc                                                                           >  1:297168/1‑139 (MQ=255)                                                    atacccacgcgAAAAAAAAGGGAAACGGGATAAAAATAAAGATTCGAAAAACCCTCTGGTTGAACGGCCAATATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTc                                                                           <  2:172137/132‑1 (MQ=255)                                                        cAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGc                                                                       <  1:750229/139‑1 (MQ=255)                                                         aGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCg                                                                      <  1:767263/139‑1 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                      <  2:191734/123‑1 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                      >  1:191734/1‑123 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGa                                                                     <  1:647701/139‑1 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGa                                                                     >  2:1043052/1‑139 (MQ=255)                                                           cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                      >  1:178172/1‑106 (MQ=255)                                                           cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                      <  2:178172/106‑1 (MQ=255)                                                            gcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAg                                                                   >  2:884002/1‑139 (MQ=255)                                                            gcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCAT‑ATCAATAAAATTGCTGGTTTGTTAAAAAAATTATTACTTGTTAATTAAAGAAGGTTCCCCCTTTATTCTTGCGagg                                                                   >  1:172137/1‑137 (MQ=255)                                                                cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATa                                                               >  1:788810/1‑139 (MQ=255)                                                                         gaATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGtttt                                                      <  2:285972/139‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGa                                                                                                                          >  2:737690/1‑71 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGa                                                                                                                          <  1:737690/71‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACa                                                                                                                   <  1:846774/78‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACa                                                                                                                   >  2:846774/1‑78 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGc                                                                       <  2:767539/121‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGc                                                                       >  1:767539/1‑121 (MQ=255)                                                                              gggATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                      >  1:125953/1‑87 (MQ=255)                                                                              gggATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                      <  2:125953/87‑1 (MQ=255)                                                                                        tAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaa                                       <  2:1475719/139‑1 (MQ=255)                                                                                         aaGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaa                                      <  2:269752/139‑1 (MQ=255)                                                                                          aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                       <  2:219546/90‑1 (MQ=255)                                                                                          aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                       >  1:219546/1‑90 (MQ=255)                                                                                          aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                       >  1:554449/1‑90 (MQ=255)                                                                                          aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                       <  2:554449/90‑1 (MQ=255)                                                                                           gCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                   <  2:582034/77‑1 (MQ=255)                                                                                           gCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                   >  1:582034/1‑77 (MQ=255)                                                                                              tCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTaaaa                                              <  1:387805/126‑1 (MQ=255)                                                                                              tCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTaaaa                                              >  2:387805/1‑126 (MQ=255)                                                                                                      aCCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGa                         <  2:788810/139‑1 (MQ=255)                                                                                                        ctcctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                       <  2:1723842/139‑1 (MQ=255)                                                                                                          cctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCa                     <  1:1694262/139‑1 (MQ=255)                                                                                                           ctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg                    <  1:1043052/139‑1 (MQ=255)                                                                                                                   ccGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAg            >  1:631149/1‑139 (MQ=255)                                                                                                                     gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg                    >  1:962450/1‑129 (MQ=255)                                                                                                                     gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg                    <  2:962450/129‑1 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTcc         >  2:205581/1‑139 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaag                                    >  2:1420683/1‑109 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaag                                    <  1:1420683/110‑2 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGAGAAGGATAAACTGTTTTTAGTAAATATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt       >  1:572700/1‑139 (MQ=255)                                                                                                                           aaTCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg                    >  1:69425/1‑123 (MQ=255)                                                                                                                           aaTCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg                    <  2:69425/123‑1 (MQ=255)                                                                                                                             tCAGTAAAATTTCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCtt  <  1:205581/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                        CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCC‑TTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCTT  >  NC_000913/515680‑515942

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