BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,487,208	(A)_6→5	coding (20/1185 nt)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,487,203	0	A	.	100.0%	1077.6 / NA	232	coding (25/1185 nt)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (108/124); total (108/124)

Only 100 of 256 total aligned reads displayed.

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

gAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGa                                                                                          <  1:5291756/93‑1 (MQ=255)

gAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTGTCATTATTTGTCATTTGCTTATGCCTTGa                                                                                          >  2:6191457/1‑93 (MQ=255)

 aaaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaa                                                                                         >  1:1480590/1‑93 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        >  2:589990/1‑93 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        <  2:4026805/93‑1 (MQ=255)

  aaGGTCCAGTAAGGACATGAGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        >  2:6156753/1‑93 (MQ=255)

   aGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaa                                                                                       >  1:4179290/1‑93 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                      <  2:5304771/93‑1 (MQ=255)

     gTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaag                                                                                     >  1:4851885/1‑92 (MQ=255)

     gTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaag                                                                                     >  2:1765775/1‑92 (MQ=255)

     gTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaag                                                                                     >  2:2832607/1‑92 (MQ=255)

       ccATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagcc                                                                                   >  2:92746/1‑90 (MQ=255)

         aTTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGa                                                                                <  1:2141521/93‑1 (MQ=255)

         aTTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGa                                                                                <  1:3320492/93‑1 (MQ=255)

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                               <  2:5140769/93‑1 (MQ=255)

            aaGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAgg                                                                             >  1:3530566/1‑93 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                            <  1:3443653/93‑1 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          >  2:4367501/1‑93 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          <  1:5809706/93‑1 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          <  1:4659305/93‑1 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                       <  2:5578154/93‑1 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                       <  1:3793875/93‑1 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                   >  2:242707/1‑93 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                   <  1:6305272/93‑1 (MQ=255)

                       tCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacac                                                                  <  1:2541085/93‑1 (MQ=255)

                       tCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacac                                                                  >  2:5458533/1‑93 (MQ=255)

                       tCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacac                                                                  <  2:605879/93‑1 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                <  1:6052478/93‑1 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                >  1:3147365/1‑93 (MQ=255)

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTc                                                               >  1:2357353/1‑93 (MQ=255)

                           atCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCg                                                              <  2:3028518/93‑1 (MQ=255)

                           atCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCg                                                              <  1:435773/93‑1 (MQ=255)

                            tCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCgg                                                             <  1:4885307/93‑1 (MQ=255)

                              gATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCa                                                           >  2:213243/1‑93 (MQ=255)

                              gATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCa                                                           >  2:2520839/1‑93 (MQ=255)

                               aTAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAt                                                          <  1:3941892/93‑1 (MQ=255)

                                tAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAtt                                                         <  1:425353/93‑1 (MQ=255)

                                 aaCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGGAttt                                                        >  2:2873998/1‑93 (MQ=255)

                                 aaCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAttt                                                        <  1:4097971/93‑1 (MQ=255)

                                 aaCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAttt                                                        <  2:134665/93‑1 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       <  2:2594540/93‑1 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       <  1:3840533/93‑1 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCAAAGGTGAAAAGACACTTCGGCATTTa                                                       <  2:6302201/93‑1 (MQ=255)

                                    taataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATg                                                     <  1:5909255/93‑1 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    <  1:3292406/93‑1 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    >  1:764956/1‑93 (MQ=255)

                                      ataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTa                                                   >  1:2250/1‑93 (MQ=255)

                                      ataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTa                                                   >  2:4273076/1‑93 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                 >  1:1048726/1‑93 (MQ=255)

                                         aGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaaa                                                <  2:1118196/93‑1 (MQ=255)

                                           ccTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATc                                              >  2:3814016/1‑93 (MQ=255)

                                           ccTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATc                                              <  2:3742304/93‑1 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           <  1:6223072/93‑1 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           >  1:1788648/1‑93 (MQ=255)

                                               cAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATa                                          >  1:335666/1‑93 (MQ=255)

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                         >  2:1923862/1‑93 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaaa                                        >  2:760747/1‑92 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAAc                                        >  1:3056278/1‑93 (MQ=255)

                                                  aCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                       >  2:1470391/1‑93 (MQ=255)

                                                   cGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACAt                                      <  2:1254933/93‑1 (MQ=255)

                                                   cGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACAt                                      <  1:4364847/93‑1 (MQ=255)

                                                    ggCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                     <  1:820390/93‑1 (MQ=255)

                                                    ggCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                     <  1:892559/93‑1 (MQ=255)

                                                      cccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGat                                   <  1:5861156/93‑1 (MQ=255)

                                                       cccttGCATTCATTAATTTTCATTTGGTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACTCTTCGGCATTTATGTAAAATCAATAACATGatt                                  <  1:2860149/90‑1 (MQ=255)

                                                       ccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatt                                  >  1:2766393/1‑93 (MQ=255)

                                                       ccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatt                                  <  1:2650949/93‑1 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                 <  1:476547/93‑1 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                 >  1:1945698/1‑93 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                 <  1:5746994/93‑1 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                >  1:308700/1‑93 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                             <  1:1432995/93‑1 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                             >  2:3929237/1‑93 (MQ=255)

                                                             cTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTAt                            >  2:468303/1‑93 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                        <  2:1265136/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     <  2:5424501/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     <  2:4892052/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     <  1:6106564/93‑1 (MQ=255)

                                                                     attTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAg                    <  1:1569889/93‑1 (MQ=255)

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                   <  1:514777/93‑1 (MQ=255)

                                                                       ttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTc                  <  1:1688268/93‑1 (MQ=255)

                                                                         gTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGa                >  2:1835073/1‑93 (MQ=255)

                                                                          tCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGAc               >  1:1142800/1‑93 (MQ=255)

                                                                            aTTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTg             >  1:4009292/1‑93 (MQ=255)

                                                                             tttGCTTTTGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACtaa            <  1:5056894/93‑3 (MQ=255)

                                                                              ttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAc           >  2:3621754/1‑93 (MQ=255)

                                                                              ttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAc           <  1:4994464/93‑1 (MQ=255)

                                                                               tGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACa          <  1:2688488/93‑1 (MQ=255)

                                                                               tGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACa          <  2:2659419/93‑1 (MQ=255)

                                                                                gCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACAt         >  1:3503909/1‑93 (MQ=255)

                                                                                gCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACAt         >  1:3634828/1‑93 (MQ=255)

                                                                                 cTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATa        >  2:1065210/1‑93 (MQ=255)

                                                                                  ttATGCCTTG‑AAAAAGCCGAAGGGGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaa       >  1:2518053/1‑93 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa      <  1:551686/93‑1 (MQ=255)

                                                                                    aTGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAAt     >  2:2286683/1‑93 (MQ=255)

                                                                                    aTGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAAt     <  1:5988353/93‑1 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc    <  2:4135554/93‑1 (MQ=255)

                                                                                      gCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATcc   <  2:3416974/93‑1 (MQ=255)

                                                                                       ccTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCca  <  2:5697102/93‑1 (MQ=255)

                                                                                       ccTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCca  <  2:60147/93‑1 (MQ=255)

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`