Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999,514	T→G	intergenic (‑34/+71)	dcyD ← / ← fliY	D‑cysteine desulfhydrase, PLP‑dependent/cystine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,514	0	T	G	100.0%	261.1 / NA	67	intergenic (‑34/+71)	dcyD/fliY	D‑cysteine desulfhydrase, PLP‑dependent/cystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/41);  total (26/41)

GATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCA  >  NC_000913/1999370‑1999657                                                                                                                                                 |                                                                                                                                                gatgaAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTgattat                                                                                                                                            <  1:558924/150‑1 (MQ=255) gatgaAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATg                                                                                                                                           >  1:109682/1‑151 (MQ=255)  atgaAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGc                                                                                                                                          >  2:756714/1‑151 (MQ=255)    gaAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCAc                                                                                                                                        <  1:542655/151‑1 (MQ=255)           tCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGaa                                                                                                                                  <  2:477140/150‑1 (MQ=255)           tCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGaa                                                                                                                                  >  2:660893/1‑150 (MQ=255)           tCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGACTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGa                                                                                                                                   >  2:602514/1‑149 (MQ=255)            cccGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGaaa                                                                                                                                 >  1:208839/1‑150 (MQ=255)              cGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAAcg                                                                                                                              <  2:744463/151‑1 (MQ=255)                    aGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGta                                                                                                                          <  1:409007/149‑1 (MQ=255)                     gATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGtaaat                                                                                                                       <  2:794366/151‑1 (MQ=255)                     gATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGtaaat                                                                                                                       <  1:297483/151‑1 (MQ=255)                      aTAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGtaaat                                                                                                                       <  2:346339/150‑1 (MQ=255)                        aaTCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGtaaataa                                                                                                                     <  2:582388/150‑1 (MQ=255)                          tCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGTAAACGCGTAAATaaaa                                                                                                                   >  2:626891/1‑150 (MQ=255)                           cAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATaaaaaa                                                                                                                 <  2:65535/151‑1 (MQ=255)                              aaaaGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAgg                                                                                                               <  2:655279/150‑1 (MQ=255)                                  gggcgGCAGATATCCGATCGGCGTTGGCGCGCCGATAAACGCCAGCCGTCGAAAACGAGTTAAACTATGCAGTGGCATAACAGCCTCCGATGTTTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTa                                                                                                          <  2:603642/149‑1 (MQ=255)                                   cgcgGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAg                                                                                                         <  2:696684/151‑1 (MQ=255)                                   cgcgGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAg                                                                                                         <  1:705302/151‑1 (MQ=255)                                         aGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAg                                                                                                   <  1:85578/151‑1 (MQ=255)                                                 gAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCttt                                                                                           <  2:379142/151‑1 (MQ=255)                                                 gAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCttt                                                                                           >  1:864805/1‑151 (MQ=255)                                                  aGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCtttt                                                                                          <  1:315412/151‑1 (MQ=255)                                                   gcggcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCttttt                                                                                         >  1:740373/1‑151 (MQ=255)                                                   gcggcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCttttt                                                                                         >  2:494392/1‑151 (MQ=255)                                                     ggcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTg                                                                                       <  2:436279/151‑1 (MQ=255)                                                     ggcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTg                                                                                       >  2:214111/1‑151 (MQ=255)                                                      gcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGt                                                                                      >  1:544654/1‑151 (MQ=255)                                                       cgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTc                                                                                     <  1:353222/151‑1 (MQ=255)                                                         ttGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCAt                                                                                   >  1:521515/1‑151 (MQ=255)                                                           ggCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTa                                                                                 <  1:776569/151‑1 (MQ=255)                                                             cgcgcCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATg                                                                               <  2:326982/151‑1 (MQ=255)                                                                gcCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTg                                                                            >  1:294874/1‑151 (MQ=255)                                                                gcCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTg                                                                            <  2:895501/151‑1 (MQ=255)                                                                 ccGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTg                                                                            >  1:320916/1‑150 (MQ=255)                                                                 ccGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTg                                                                            >  2:347705/1‑150 (MQ=255)                                                                     tAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGatta                                                                        <  2:590118/150‑1 (MQ=255)                                                                              gCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTcaca                                                              <  1:450650/151‑1 (MQ=255)                                                                                    gAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCa                                                        <  2:204659/151‑1 (MQ=255)                                                                                     aaaaCGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCAc                                                       >  1:460419/1‑151 (MQ=255)                                                                                         cGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCaa                                                    >  1:236828/1‑150 (MQ=255)                                                                                           ggTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAAcc                                                 >  1:81703/1‑151 (MQ=255)                                                                                               aaaTTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCAttt                                             <  2:266039/151‑1 (MQ=255)                                                                                                 aTTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCAttttt                                           <  1:129435/151‑1 (MQ=255)                                                                                                 aTTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCAttttt                                           <  1:660908/151‑1 (MQ=255)                                                                                                   tATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTc                                          <  2:418357/150‑1 (MQ=255)                                                                                                    aTGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCgg                                        <  2:811554/151‑1 (MQ=255)                                                                                                    aTGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCg                                         <  2:632665/150‑1 (MQ=255)                                                                                                       cAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCGAACCATTTTTCGGaaa                                     <  2:594371/151‑1 (MQ=255)                                                                                                           ggCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAggg                                  <  1:789204/150‑1 (MQ=255)                                                                                                             cATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCtt                               >  2:764962/1‑151 (MQ=255)                                                                                                              aTAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTg                              >  2:129125/1‑151 (MQ=255)                                                                                                                aaCAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCa                            <  1:466925/151‑1 (MQ=255)                                                                                                                  cAGCCTCCGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCaga                          <  1:507118/151‑1 (MQ=255)                                                                                                                        ccGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGcc                     >  2:39922/1‑150 (MQ=255)                                                                                                                         cGATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCAt                   >  1:446938/1‑151 (MQ=255)                                                                                                                          gATGTGTGTTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATc                  <  2:7162/151‑1 (MQ=255)                                                                                                                             gtgtgtTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCttt               <  2:109681/151‑1 (MQ=255)                                                                                                                               gtgtTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCttttt             >  1:429852/1‑151 (MQ=255)                                                                                                                               gtgtTGTTGTGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCtttt              <  2:668589/150‑1 (MQ=255)                                                                                                                                 gttgttgtGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGc           >  1:293327/1‑151 (MQ=255)                                                                                                                                    gttgtGATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCAtt        <  2:619495/151‑1 (MQ=255)                                                                                                                                      tgtgATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCAttt       <  1:513891/150‑1 (MQ=255)                                                                                                                                       gtgATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTcc     <  2:463085/151‑1 (MQ=255)                                                                                                                                         gATTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGc   <  1:60788/151‑1 (MQ=255)                                                                                                                                          aTTTTCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGc   >  1:343968/1‑150 (MQ=255)                                                                                                                                           ttttCTGATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCa  >  2:334508/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                GATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCA  >  NC_000913/1999370‑1999657

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