BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,053,140	A→G	W75R (TGG→CGG)	yccM ←	putative 4Fe‑4S membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,140	0	A	G	100.0%	1038.7 / NA	245	W75R (TGG→CGG)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (128/117); total (128/117)

Only 100 of 254 total aligned reads displayed.

AAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTC  >  NC_000913/1053042‑1053237

  aaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCt                                                                                               <  1:1554260/101‑1 (MQ=255)

  aaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCt                                                                                               <  1:1554486/101‑1 (MQ=255)

   aGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCtt                                                                                              >  2:1554883/1‑101 (MQ=255)

   aGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCtt                                                                                              <  1:1554567/101‑1 (MQ=255)

    gCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTc                                                                                             >  1:1554885/1‑101 (MQ=255)

    gCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTc                                                                                             <  2:1554432/101‑1 (MQ=255)

      gACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGTTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGa                                                                                           >  1:1554888/1‑101 (MQ=255)

      gACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGa                                                                                           >  2:1554887/1‑101 (MQ=255)

        ccAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATa                                                                                         >  2:1554891/1‑101 (MQ=255)

        ccAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATa                                                                                         <  1:1554427/101‑1 (MQ=255)

          aaTAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATAcc                                                                                       >  1:1554893/1‑101 (MQ=255)

            tAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTc                                                                                     >  2:1554898/1‑101 (MQ=255)

            tAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTc                                                                                     <  2:1554193/101‑1 (MQ=255)

             aaCAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTcc                                                                                    >  1:1554899/1‑101 (MQ=255)

                aGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCgg                                                                                 >  2:1554902/1‑101 (MQ=255)

                aGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCgg                                                                                 <  2:1554396/101‑1 (MQ=255)

                aGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCgg                                                                                 <  2:1554274/101‑1 (MQ=255)

                 gCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGt                                                                                >  1:1554904/1‑101 (MQ=255)

                   aTCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCt                                                                              >  1:1554906/1‑101 (MQ=255)

                     cgcgGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTgg                                                                            <  2:1554462/101‑1 (MQ=255)

                      gcggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGc                                                                           >  2:1554910/1‑101 (MQ=255)

                      gcggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGc                                                                           <  1:1554607/101‑1 (MQ=255)

                       cggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCg                                                                          >  2:1554912/1‑101 (MQ=255)

                        ggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCt                                                                              <  1:1554349/96‑1 (MQ=255)

                         gcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGAc                                                                        >  2:1554916/1‑101 (MQ=255)

                           ggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACaa                                                                      >  2:1554919/1‑101 (MQ=255)

                           ggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACaa                                                                      >  1:1554918/1‑101 (MQ=255)

                              gTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATg                                                                   >  1:1554921/1‑101 (MQ=255)

                              gTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATg                                                                   <  1:1554653/101‑1 (MQ=255)

                               tGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGt                                                                  >  1:1554924/1‑101 (MQ=255)

                                ggacggacGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTc                                                                 <  1:1554505/101‑1 (MQ=255)

                                ggacggacGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTc                                                                 <  1:1554531/101‑1 (MQ=255)

                                ggacggacGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTc                                                                 >  2:1554927/1‑101 (MQ=255)

                                  acggacggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGt                                                               <  1:1554514/101‑1 (MQ=255)

                                    ggacggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGc                                                             >  2:1554933/1‑101 (MQ=255)

                                    ggacggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGc                                                             >  1:1554931/1‑101 (MQ=255)

                                     gacggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCt                                                            >  1:1554934/1‑101 (MQ=255)

                                      acggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTa                                                           <  1:1554311/101‑1 (MQ=255)

                                       cggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTccc                                                                                   <  2:1554404/76‑1 (MQ=255)

                                        ggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACt                                                         <  2:1554591/101‑1 (MQ=255)

                                        ggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACt                                                         <  2:1554574/101‑1 (MQ=255)

                                          tAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTg                                                       >  2:1554939/1‑101 (MQ=255)

                                            aCTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGcc                                                     >  1:1554940/1‑101 (MQ=255)

                                            aCTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGcc                                                     <  1:1554461/101‑1 (MQ=255)

                                               gACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGtt                                                  >  2:1554943/1‑101 (MQ=255)

                                               gACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGtt                                                  <  1:1554586/101‑1 (MQ=255)

                                               gACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGtt                                                  <  2:1554533/101‑1 (MQ=255)

                                                aCCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGttt                                                 >  1:1554944/1‑101 (MQ=255)

                                                aCCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGttt                                                 <  2:1554494/101‑1 (MQ=255)

                                                 ccTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTc                                                <  1:1554339/101‑1 (MQ=255)

                                                   tGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGt                                              >  2:1554949/1‑101 (MQ=255)

                                                   tGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGt                                              >  1:1554950/1‑101 (MQ=255)

                                                   tGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGt                                              <  1:1554464/101‑1 (MQ=255)

                                                    gTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTa                                             >  1:1554953/1‑101 (MQ=255)

                                                     ttACAAGGCTATATTTCAGATTCATCAGCCCGGTAATCGGTAGCCGGCCTTCGATACCt                                                                                      >  2:1554957/1‑59 (MQ=255)

                                                     ttACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTaa                                            <  1:1554439/101‑1 (MQ=255)

                                                      tACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAAt                                           >  1:1554960/1‑101 (MQ=255)

                                                      tACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAAt                                           >  1:1554958/1‑101 (MQ=255)

                                                      tACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAAt                                           <  2:1554677/101‑1 (MQ=255)

                                                        cAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAg                                         <  1:1554784/101‑1 (MQ=255)

                                                         aaGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAgc                                        <  1:1554397/101‑1 (MQ=255)

                                                         aaGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAgc                                        <  1:1554634/101‑1 (MQ=255)

                                                          aGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAgcg                                       >  2:1554964/1‑101 (MQ=255)

                                                            gCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCa                                     >  1:1554965/1‑101 (MQ=255)

                                                               atatTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCAccc                                  <  1:1554539/101‑1 (MQ=255)

                                                                 atTTCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCaa                                >  2:1554968/1‑101 (MQ=255)

                                                                  tttCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAAt                               <  2:1554453/101‑1 (MQ=255)

                                                                   ttCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATa                              >  1:1554971/1‑101 (MQ=255)

                                                                    tCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATAc                             <  2:1554699/101‑1 (MQ=255)

                                                                    tCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATAc                             <  1:1554700/101‑1 (MQ=255)

                                                                    tCAGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATAc                             >  2:1554972/1‑101 (MQ=255)

                                                                      aGATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCa                           <  2:1554508/101‑1 (MQ=255)

                                                                       gATTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAg                          >  1:1554976/1‑101 (MQ=255)

                                                                        aTTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGa                         >  2:1554979/1‑101 (MQ=255)

                                                                        aTTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGa                         >  1:1554982/1‑101 (MQ=255)

                                                                        aTTCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGa                         <  2:1554423/101‑1 (MQ=255)

                                                                         ttCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGaa                        >  2:1554983/1‑101 (MQ=255)

                                                                         ttCATCAGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGaa                        <  1:1554573/101‑1 (MQ=255)

                                                                          tcatcaGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAAc                       <  1:1554850/101‑1 (MQ=255)

                                                                           catcaGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACg                      >  1:1554985/1‑101 (MQ=255)

                                                                             tcaGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTg                    <  2:1554816/101‑1 (MQ=255)

                                                                               aGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGAt                  <  1:1554541/101‑1 (MQ=255)

                                                                               aGCCCGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGAt                  >  2:1554989/1‑101 (MQ=255)

                                                                                  ccGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGc               >  1:1554992/1‑101 (MQ=255)

                                                                                   cGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCa              <  2:1554587/101‑1 (MQ=255)

                                                                                   cGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCa              >  2:1554994/1‑101 (MQ=255)

                                                                                   cGGCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCa              <  1:1554479/101‑1 (MQ=255)

                                                                                    ggCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGtt                                                  <  2:1554648/64‑1 (MQ=255)

                                                                                    ggCAATCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCaa             >  1:1554996/1‑101 (MQ=255)

                                                                                       aaTCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAatat          >  2:1554998/1‑101 (MQ=255)

                                                                                        aTCGGTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAAtata         <  1:1554549/101‑1 (MQ=255)

                                                                                           ggTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAat      >  2:1555003/1‑101 (MQ=255)

                                                                                           ggTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAat      >  1:1555002/1‑101 (MQ=255)

                                                                                           ggTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAat      <  1:1554400/101‑1 (MQ=255)

                                                                                            gTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCTAATACCAGAACGTGATGGCAATATAAata     <  1:1554559/101‑1 (MQ=255)

                                                                                            gTAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAata     >  1:1555007/1‑101 (MQ=255)

                                                                                             tAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAatat    <  2:1554550/101‑1 (MQ=255)

                                                                                             tAGCCGGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAatat    >  2:1555009/1‑101 (MQ=255)

                                                                                              agccggccTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGt                                              <  1:1554482/58‑1 (MQ=255)

                                                                                              agccggccTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATAtt   <  1:1554524/101‑1 (MQ=255)

AAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTC  >  NC_000913/1053042‑1053237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`