BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,288,760	C→T	P111L (CCG→CTG)	ychK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,288,760	0	C	T	99.5%	757.6 / ‑4.0	204	P111L (CCG→CTG)	ychK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (101/102); minor base G (1/0); total (102/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 209 total aligned reads displayed.

GGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGC  >  NC_000913/1288670‑1288848

 ggATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGa                                                                                       <  2:4431026/93‑1 (MQ=255)

     gTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACa                                                                                   >  2:1984445/1‑93 (MQ=255)

     gTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACa                                                                                   <  2:2371166/93‑1 (MQ=255)

     gTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACa                                                                                   >  2:5628291/1‑93 (MQ=255)

     gTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACa                                                                                   >  2:5889590/1‑93 (MQ=255)

      ttttACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACag                                                                                  <  2:5821763/93‑1 (MQ=255)

       tttACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACaga                                                                                 <  1:5050070/93‑1 (MQ=255)

         tACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAgaga                                                                               <  2:4359595/93‑1 (MQ=255)

         tACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAgaga                                                                               <  1:1328764/93‑1 (MQ=255)

           cGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATc                                                                             >  1:5667022/1‑93 (MQ=255)

            gCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCg                                                                            <  2:6249748/93‑1 (MQ=255)

            gCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCg                                                                            <  2:463953/93‑1 (MQ=255)

             tcTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGa                                                                           >  2:1637991/2‑93 (MQ=255)

              cTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGaa                                                                          <  1:2697394/93‑1 (MQ=255)

               tGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGaaa                                                                         >  1:243232/1‑93 (MQ=255)

               tGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGaaa                                                                         >  2:3970732/1‑93 (MQ=255)

                    gATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGt                                                                    >  1:429563/1‑93 (MQ=255)

                      tctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTc                                                                  <  1:3516957/93‑1 (MQ=255)

                      tctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTc                                                                  <  2:3777454/93‑1 (MQ=255)

                       ctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTcc                                                                 <  2:1438305/93‑1 (MQ=255)

                       ctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTcc                                                                 <  2:6205270/93‑1 (MQ=255)

                       ctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTcc                                                                 >  1:3129885/1‑93 (MQ=255)

                       ctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTcc                                                                 >  1:1124535/1‑93 (MQ=255)

                        tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTccc                                                                <  1:5850153/93‑1 (MQ=255)

                        tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTccc                                                                <  1:2144718/93‑1 (MQ=255)

                        tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTccc                                                                <  1:4049221/93‑1 (MQ=255)

                        tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTccc                                                                <  1:1055351/93‑1 (MQ=255)

                          tcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCgc                                                              >  1:6219263/1‑92 (MQ=255)

                           ccTGGCAGCGCGGTGGGGTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTc                                                             >  1:632819/1‑93 (MQ=255)

                               gCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCtt                                                         <  2:2449962/93‑1 (MQ=255)

                                cAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCttt                                                        >  1:258437/1‑93 (MQ=255)

                                 aGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTg                                                       >  2:5871591/1‑93 (MQ=255)

                                  gcgcGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGc                                                      <  2:2208293/93‑1 (MQ=255)

                                  gcgcGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGc                                                      <  1:5221293/93‑1 (MQ=255)

                                     cgGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGc                                                   <  2:4456684/93‑1 (MQ=255)

                                     cgGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGc                                                   >  2:5053011/1‑93 (MQ=255)

                                     cgGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGc                                                   <  1:84640/93‑1 (MQ=255)

                                      ggTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCt                                                  >  2:5774854/1‑93 (MQ=255)

                                      ggTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCt                                                  >  1:1572406/1‑93 (MQ=255)

                                        tGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGt                                                <  1:4779033/93‑1 (MQ=255)

                                         gggTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGtt                                               <  2:2622430/93‑1 (MQ=255)

                                         gggTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGtt                                               <  1:2826517/93‑1 (MQ=255)

                                            ttccTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGcc                                            <  2:6083122/90‑1 (MQ=255)

                                            ttACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGcc                                            <  2:3042621/93‑1 (MQ=255)

                                              aCTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGccac                                          <  2:3344615/93‑1 (MQ=255)

                                              aCTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGccac                                          >  1:4503260/1‑93 (MQ=255)

                                              aCTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGATGTTGccac                                          >  2:2296330/1‑93 (MQ=255)

                                                 gcgcGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCaa                                       <  1:2274814/93‑1 (MQ=255)

                                                  cgcgGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCTGCTGTTGCCACCAAt                                      <  2:4231238/93‑1 (MQ=255)

                                                   gcggcgAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAAtt                                     >  2:4008863/1‑93 (MQ=255)

                                                     ggcgAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTa                                   >  1:5289061/1‑93 (MQ=255)

                                                     ggcgAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTa                                   >  1:1875931/1‑93 (MQ=255)

                                                        gAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGACAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGt                                <  1:4446386/93‑1 (MQ=255)

                                                        gAGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGt                                >  2:1411756/1‑93 (MQ=255)

                                                         aGCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTa                               >  1:3763530/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  2:3585983/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  2:4117485/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  2:1004252/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  2:1787721/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  2:1052949/1‑93 (MQ=255)

                                                          gCGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                              >  1:68599/1‑93 (MQ=255)

                                                           cGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACg                             <  1:2865540/93‑1 (MQ=255)

                                                           cGTGTCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACg                             <  1:5071230/93‑1 (MQ=255)

                                                            gtgtCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACgg                            <  2:3000178/93‑1 (MQ=255)

                                                             tgtCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACggg                           >  1:574982/1‑93 (MQ=255)

                                                              gtCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGa                          >  1:111602/1‑93 (MQ=255)

                                                               tCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                         >  1:923030/1‑93 (MQ=255)

                                                               tCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                         >  1:1764370/1‑93 (MQ=255)

                                                               tCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                         >  1:78570/1‑93 (MQ=255)

                                                               tCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                         <  1:986074/93‑1 (MQ=255)

                                                               tCTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                         >  1:3395403/1‑93 (MQ=255)

                                                                cTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACg                        >  2:3111562/1‑93 (MQ=255)

                                                                cTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACg                        >  2:3322067/1‑93 (MQ=255)

                                                                cTTCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACg                        >  1:5893751/1‑93 (MQ=255)

                                                                 ttCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGt                       <  2:1255486/93‑1 (MQ=255)

                                                                 ttCAATCAATACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGt                       <  2:3799919/93‑1 (MQ=255)

                                                                  tcaatcaaTACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTg                      >  2:5854795/1‑93 (MQ=255)

                                                                  tcaatcaaTACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTg                      >  1:4172775/1‑93 (MQ=255)

                                                                  tcaatcaaTACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTg                      >  2:1290737/1‑93 (MQ=255)

                                                                    aatcaatACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTgca                    >  1:4708409/1‑91 (MQ=255)

                                                                    aatcaatACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGaa                    <  1:2172731/93‑1 (MQ=255)

                                                                      tcaatACCGCGAAATAATGCTGGAAACAGCGATCGAAAATTGTTCCTGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGCTGTGAAtt                  >  1:4586668/1‑93 (MQ=255)

                                                                       caatACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTa                 <  2:1951398/93‑1 (MQ=255)

                                                                        aatACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTAt                >  2:4327056/1‑93 (MQ=255)

                                                                         atACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATg               >  1:2848767/1‑93 (MQ=255)

                                                                          tACCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATgg              <  2:3548908/93‑1 (MQ=255)

                                                                           aCCGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGt             <  2:5585721/93‑1 (MQ=255)

                                                                            ccGCGAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGtt            <  2:2473496/93‑1 (MQ=255)

                                                                             cgcgAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGttt           <  1:2049890/93‑1 (MQ=255)

                                                                                gAAATAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACt        >  1:5526898/1‑93 (MQ=255)

                                                                                 aaaTAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTg       >  2:4945635/1‑93 (MQ=255)

                                                                                 aaaTAATGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTg       >  2:2876658/1‑93 (MQ=255)

                                                                                  aataatGCGGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGa      >  2:6136245/1‑93 (MQ=255)

                                                                                   ataatGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGaa     <  2:2776703/93‑1 (MQ=255)

                                                                                   ataatGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGaa     <  1:1411756/93‑1 (MQ=255)

                                                                                    taatGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAg    >  2:2238405/1‑93 (MQ=255)

                                                                                     aatGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAgg   <  1:250805/93‑1 (MQ=255)

                                                                                     aatGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAgg   <  1:4369218/93‑1 (MQ=255)

                                                                                      atGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGc  <  1:6023846/93‑1 (MQ=255)

                                                                                      atGCTGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGc  <  1:5296912/93‑1 (MQ=255)

GGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGC  >  NC_000913/1288670‑1288848

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