BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,140,586	G→A	Q72* (CAA→TAA)	hybE ←	hydrogenase 2‑specific chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,140,586	0	G	A	100.0%	1604.7 / NA	373	Q72* (CAA→TAA)	hybE	hydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (197/176); total (197/176)

Only 100 of 373 total aligned reads displayed.

TTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACC  >  NC_000913/3140488‑3140683

ttGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAAt                                                                                                 <  1:5793227/101‑1 (MQ=255)

ttGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAAt                                                                                                 >  2:5793824/1‑101 (MQ=255)

    gAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATccgg                                                                                             <  1:5792971/101‑1 (MQ=255)

    gAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATccgg                                                                                             >  1:5793827/1‑101 (MQ=255)

       acacCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATccggg                                                                                            <  1:5792613/99‑1 (MQ=255)

        cacCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATccgggccg                                                                                         <  2:5793377/101‑1 (MQ=255)

         acCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATccgggccgg                                                                                        >  2:5793832/1‑101 (MQ=255)

            gTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGa                                                                                     >  1:5793839/1‑101 (MQ=255)

            gTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGa                                                                                     >  2:5793838/1‑101 (MQ=255)

              ccAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAg                                                                                   <  1:5793289/101‑1 (MQ=255)

                aGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGAt                                                                                 <  1:5792906/101‑1 (MQ=255)

                 gTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATa                                                                                >  1:5793845/1‑101 (MQ=255)

                  ttCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATaa                                                                               >  2:5793848/1‑101 (MQ=255)

                   tCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAAc                                                                              >  2:5793852/1‑101 (MQ=255)

                    cGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACt                                                                             >  1:5793854/1‑101 (MQ=255)

                       cAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCa                                                                          >  1:5793860/1‑101 (MQ=255)

                        aaCGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCAc                                                                         >  2:5793861/1‑101 (MQ=255)

                          cGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTc                                                                       >  2:5793869/1‑101 (MQ=255)

                           ggTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCa                                                                      >  2:5793873/1‑101 (MQ=255)

                           ggTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCa                                                                      >  1:5793874/1‑101 (MQ=255)

                              aaaGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCa                                                                   >  2:5793881/1‑101 (MQ=255)

                              aaaGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCa                                                                   >  1:5793882/1‑101 (MQ=255)

                              aaaGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCa                                                                   <  1:5793119/101‑1 (MQ=255)

                                aGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGCTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATc                                                                 <  1:5793250/101‑1 (MQ=255)

                                 ggTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATcc                                                                >  2:5793889/1‑101 (MQ=255)

                                 ggTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATcc                                                                <  2:5793337/101‑1 (MQ=255)

                                 ggTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATcc                                                                <  1:5793126/101‑1 (MQ=255)

                                  gTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCa                                                               <  2:5793388/101‑1 (MQ=255)

                                  gTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCa                                                               >  1:5793890/1‑101 (MQ=255)

                                    cATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACg                                                             >  1:5793896/1‑101 (MQ=255)

                                     aTAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgg                                                            >  2:5793901/1‑101 (MQ=255)

                                     aTAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgg                                                            <  1:5793315/101‑1 (MQ=255)

                                       aGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACgggg                                                          <  2:5793217/101‑1 (MQ=255)

                                        gTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGt                                                         <  2:5793169/101‑1 (MQ=255)

                                           ccATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGAt                                                      >  1:5793904/1‑101 (MQ=255)

                                           ccATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGAt                                                      <  2:5793628/101‑1 (MQ=255)

                                             aTACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATca                                                    >  1:5793909/1‑101 (MQ=255)

                                              tACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcac                                                   >  2:5793911/1‑101 (MQ=255)

                                               aCGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCTAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGc                                                                    <  2:5792950/83‑1 (MQ=255)

                                               aCGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATcc                                                                <  2:5793321/87‑1 (MQ=255)

                                               aCGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcaca                                                  >  1:5793914/1‑101 (MQ=255)

                                               aCGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcaca                                                  <  2:5793300/101‑1 (MQ=255)

                                                cGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcacac                                                 >  2:5793917/1‑101 (MQ=255)

                                                 ggCAGTCGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcacaca                                                <  1:5793756/101‑1 (MQ=255)

                                                   cAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcc                                              <  1:5793452/101‑1 (MQ=255)

                                                   cAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcc                                              >  2:5793923/1‑101 (MQ=255)

                                                    aGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAccc                                             >  1:5793926/1‑101 (MQ=255)

                                                     gTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcccc                                            >  2:5793928/1‑101 (MQ=255)

                                                     gTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcccc                                            <  2:5793283/101‑1 (MQ=255)

                                                     gTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACAcccc                                            <  2:5793072/101‑1 (MQ=255)

                                                      ttGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGTGTGATcacaca                                                <  1:5793548/96‑1 (MQ=255)

                                                      ttGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCg                                           <  2:5793703/101‑1 (MQ=255)

                                                       tGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATcacaca                                                >  2:5793931/1‑95 (MQ=255)

                                                         cAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACggg                                                           >  1:5793938/1‑82 (MQ=255)

                                                         cAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTcc                                        <  1:5793270/101‑1 (MQ=255)

                                                          aGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCa                                       <  2:5793515/101‑1 (MQ=255)

                                                          aGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCa                                       >  2:5793939/1‑101 (MQ=255)

                                                           gACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCAc                                      <  1:5793035/101‑1 (MQ=255)

                                                            aCCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACt                                     <  2:5793355/101‑1 (MQ=255)

                                                              cAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGa                                   >  2:5793946/1‑101 (MQ=255)

                                                              cAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGa                                   >  1:5793947/1‑101 (MQ=255)

                                                                aTTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGAcc                                 >  1:5793953/1‑101 (MQ=255)

                                                                aTTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGAcc                                 >  2:5793954/1‑101 (MQ=255)

                                                                 ttttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCAc                                                              <  1:5793369/71‑1 (MQ=255)

                                                                 ttttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCt                                >  2:5793958/1‑101 (MQ=255)

                                                                  tttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCtt                               <  1:5793526/101‑1 (MQ=255)

                                                                   ttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                              <  1:5793237/101‑1 (MQ=255)

                                                                   ttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                              <  2:5793492/101‑1 (MQ=255)

                                                                   ttttCACTCACTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTc                              >  1:5793963/1‑101 (MQ=255)

                                                                      tcactcacTTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGaa                           >  2:5793966/1‑101 (MQ=255)

                                                                        actcactTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACa                         <  2:5793418/101‑1 (MQ=255)

                                                                        actcactTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACa                         >  1:5793971/1‑101 (MQ=255)

                                                                          tcactTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGc                       <  1:5793615/101‑1 (MQ=255)

                                                                           cactTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCg                      >  1:5793975/1‑101 (MQ=255)

                                                                            actTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGt                     >  2:5793978/1‑101 (MQ=255)

                                                                             ctTTGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGCGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTa                    >  1:5793981/1‑101 (MQ=255)

                                                                              tttGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaa                   >  1:5793987/1‑101 (MQ=255)

                                                                               ttGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaa                  <  1:5793392/101‑1 (MQ=255)

                                                                               ttGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaa                  >  2:5793992/1‑101 (MQ=255)

                                                                               ttGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaa                  >  1:5793991/1‑101 (MQ=255)

                                                                                tGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCa                                                                      >  1:5793995/1‑48 (MQ=255)

                                                                                tGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaaa                 <  2:5793140/101‑1 (MQ=255)

                                                                                tGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaaa                 <  2:5793201/101‑1 (MQ=255)

                                                                                tGCGCAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTaaaa                 <  2:5793595/101‑1 (MQ=255)

                                                                                   gcAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAAtta              <  1:5793336/101‑3 (MQ=255)

                                                                                   gcAGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCa              >  2:5793999/1‑101 (MQ=255)

                                                                                     aGCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGa            >  1:5794003/1‑101 (MQ=255)

                                                                                      gCGGCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGaa           <  2:5793474/101‑1 (MQ=255)

                                                                                        ggCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAAc         >  1:5794010/1‑101 (MQ=255)

                                                                                        ggCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAAc         >  2:5794009/1‑101 (MQ=255)

                                                                                        ggCCAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAAc         <  1:5793203/101‑1 (MQ=255)

                                                                                          ccAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACAt       <  1:5793347/101‑1 (MQ=255)

                                                                                           cAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCAc                                      <  2:5793457/69‑1 (MQ=255)

                                                                                           cAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATa      >  2:5794018/1‑101 (MQ=255)

                                                                                           cAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATa      <  2:5793245/101‑1 (MQ=255)

                                                                                           cAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATa      <  2:5792824/101‑1 (MQ=255)

                                                                                           cAGAGTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATa      >  1:5794017/1‑101 (MQ=255)

                                                                                            agagTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATac     <  1:5793069/101‑1 (MQ=255)

                                                                                             gagTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGGGATCACACAcc                                              >  1:5794023/1‑59 (MQ=255)

                                                                                               gTTAATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACAcc  <  2:5793576/101‑1 (MQ=255)

TTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGATCACACACCCCGTCCACTGACCTTCGAACAGCGTAAAATCAGAAACATACACC  >  NC_000913/3140488‑3140683

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: `‑`