BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,841,925	T→C	H324R (CAC→CGC)	b1762 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,841,925	0	T	C	84.9%	541.6 / 85.3	212	H324R (CAC→CGC)	b1762	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (19/13); new base C (85/95); total (104/108)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 216 total aligned reads displayed.

CCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAT  >  NC_000913/1841834‑1842017

ccccGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTg                                                                                             >  2:1692695/1‑93 (MQ=255)

 cccGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCTCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAgcgc                                                                                            >  1:294077/1‑93 (MQ=255)

 cccGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAgcgc                                                                                            >  2:4703716/1‑93 (MQ=255)

    ggCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTcc                                                                                         <  2:5108073/93‑1 (MQ=255)

     gCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTccc                                                                                        <  1:153261/93‑1 (MQ=255)

       aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCaa                                                                                      <  2:4267361/93‑1 (MQ=255)

       aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCaa                                                                                      >  1:5666577/1‑93 (MQ=255)

       aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCaa                                                                                      >  1:331097/1‑93 (MQ=255)

        aaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCaaa                                                                                     >  1:4525512/1‑93 (MQ=255)

        aaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCaaa                                                                                     >  2:2393120/1‑93 (MQ=255)

         aaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCaaaa                                                                                    >  1:4220787/1‑93 (MQ=255)

             tCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGc                                                                                <  1:4409644/93‑1 (MQ=255)

              cAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCt                                                                               <  2:3152700/93‑1 (MQ=255)

              cAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCt                                                                               <  1:4330739/93‑1 (MQ=255)

               aGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCtt                                                                              <  1:2186894/93‑1 (MQ=255)

                  ggTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTc                                                                           <  2:5284889/92‑1 (MQ=255)

                  cGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTc                                                                           <  1:1264952/93‑1 (MQ=255)

                   gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCg                                                                          <  1:1094283/93‑1 (MQ=255)

                    tttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGt                                                                         <  1:4947177/93‑1 (MQ=255)

                     ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTa                                                                        <  2:4219892/93‑1 (MQ=255)

                     ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTa                                                                        >  1:4835133/1‑93 (MQ=255)

                     ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTa                                                                        >  2:3921700/1‑93 (MQ=255)

                        tGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACg                                                                     >  2:2095457/1‑93 (MQ=255)

                        tGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACGGTAGCGCTCCCAAAACTGCTTGTCGTAACg                                                                     >  2:648978/1‑93 (MQ=255)

                         gAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGa                                                                    <  1:3418712/93‑1 (MQ=255)

                           gTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATa                                                                  >  1:3665590/1‑93 (MQ=255)

                             aTTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAac                                                                <  2:5681984/93‑1 (MQ=255)

                             aTTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAac                                                                <  2:938449/93‑1 (MQ=255)

                             aTTTAACGAATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAac                                                                <  2:1385374/93‑1 (MQ=255)

                              tttAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAaca                                                               <  1:3854230/93‑1 (MQ=255)

                              tttAACGGATGAAGCTCCCCGCGGAACATTTCTTATGGGCCACCGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAaca                                                               <  2:2286800/93‑1 (MQ=255)

                                tAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACg                                                             >  1:1503798/1‑93 (MQ=255)

                                tAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACg                                                             <  2:5526747/93‑1 (MQ=255)

                                 aaCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGc                                                            >  2:2900627/1‑93 (MQ=255)

                                  aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCt                                                           >  2:5722558/1‑93 (MQ=255)

                                    ggATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTc                                                         <  2:3014926/93‑1 (MQ=255)

                                     gatgatGGTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                        >  1:6315725/1‑93 (MQ=255)

                                      atgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCaa                                                       <  2:5614652/93‑1 (MQ=255)

                                       tgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAg                                                      >  1:5521292/1‑93 (MQ=255)

                                       tgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAg                                                      <  1:2588072/93‑1 (MQ=255)

                                       tgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAg                                                      >  2:4961990/1‑93 (MQ=255)

                                           gCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCa                                                  <  2:1111518/93‑1 (MQ=255)

                                            cTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAg                                                 >  1:6613162/1‑93 (MQ=255)

                                            cTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAg                                                 <  2:3145894/93‑1 (MQ=255)

                                              ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                               >  2:3214448/1‑93 (MQ=255)

                                               cccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                              <  1:4903003/93‑1 (MQ=255)

                                               cccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                              >  2:1682849/1‑93 (MQ=255)

                                                ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATc                                             >  2:1018989/1‑93 (MQ=255)

                                                ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATc                                             >  2:4850790/1‑93 (MQ=255)

                                                 cACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATcc                                            <  1:813440/93‑1 (MQ=255)

                                                    ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                         <  1:1911747/93‑1 (MQ=255)

                                                    ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                         <  2:214532/93‑1 (MQ=255)

                                                    ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                         >  1:5758137/1‑93 (MQ=255)

                                                     gAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTa                                        >  1:5341744/1‑93 (MQ=255)

                                                     gAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTa                                        >  1:202638/1‑93 (MQ=255)

                                                       aCATTTCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAAc                                      >  1:5118151/1‑93 (MQ=255)

                                                         aTTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACtt                                    <  1:3513995/93‑1 (MQ=255)

                                                          tttCTTATGGGCCAACGGCCTTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACttt                                   >  1:6464922/1‑93 (MQ=255)

                                                          tttCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACttt                                   >  2:5870775/1‑93 (MQ=255)

                                                            tCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCt                                 <  2:1830119/93‑1 (MQ=255)

                                                            tCTTATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCt                                 >  1:2501630/1‑93 (MQ=255)

                                                               tATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCg                              >  2:5710698/1‑93 (MQ=255)

                                                               tATGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCg                              >  1:4661402/1‑93 (MQ=255)

                                                                aTGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGa                             >  1:171676/1‑93 (MQ=255)

                                                                aTGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGa                             <  2:3481981/93‑1 (MQ=255)

                                                                aTGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGa                             <  1:4160583/93‑1 (MQ=255)

                                                                 tGGGCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                            <  2:5981265/93‑1 (MQ=255)

                                                                  gggCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                           <  2:1948589/93‑1 (MQ=255)

                                                                  gggCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGCCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                           >  1:371850/1‑93 (MQ=255)

                                                                   ggCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATag                          <  2:2471551/93‑1 (MQ=255)

                                                                   ggCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATag                          <  2:1665591/93‑1 (MQ=255)

                                                                    gCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagc                         <  2:3274076/93‑1 (MQ=255)

                                                                    gCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagc                         <  2:2793782/93‑1 (MQ=255)

                                                                    gCCAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagc                         >  2:896609/1‑93 (MQ=255)

                                                                     ccAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagca                        >  2:3652546/1‑93 (MQ=255)

                                                                     ccAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagca                        <  1:2435154/93‑1 (MQ=255)

                                                                      cAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcag                       <  2:3575085/93‑1 (MQ=255)

                                                                      cAACGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcag                       >  1:6163182/1‑93 (MQ=255)

                                                                       aaCGGCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcagc                      >  1:655939/1‑93 (MQ=255)

                                                                         cGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGcagcag                    >  1:5046898/1‑93 (MQ=255)

                                                                          ggCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGa                   >  2:4928175/1‑93 (MQ=255)

                                                                           gCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACTCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat                  >  2:766510/1‑93 (MQ=255)

                                                                           gCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat                  >  2:3668063/1‑93 (MQ=255)

                                                                           gCATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat                  <  2:4742295/93‑1 (MQ=255)

                                                                            cATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                 <  2:4560220/93‑1 (MQ=255)

                                                                            cATTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                 <  2:4454330/93‑1 (MQ=255)

                                                                             aTTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGatat                >  2:1998206/1‑93 (MQ=255)

                                                                             aTTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGatat                >  2:1597481/1‑93 (MQ=255)

                                                                             aTTTCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGatat                >  1:3008972/1‑93 (MQ=255)

                                                                              tttCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg               >  1:2666981/1‑93 (MQ=255)

                                                                              tttCTTACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg               >  2:652017/1‑93 (MQ=255)

                                                                                  ttACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAg           <  1:4873931/93‑1 (MQ=255)

                                                                                   tACTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAga          >  1:3092778/1‑92 (MQ=255)

                                                                                     cTGTAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaa        <  2:721160/93‑1 (MQ=255)

                                                                                      tGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaaa       <  2:6126677/93‑1 (MQ=255)

                                                                                       gTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaaaa      >  1:1888029/1‑92 (MQ=255)

                                                                                       gTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAg      >  1:526363/1‑93 (MQ=255)

                                                                                        tAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGa     <  1:5940322/93‑1 (MQ=255)

                                                                                        tAGCGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGa     >  1:2321024/1‑93 (MQ=255)

                                                                                           tGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt  <  2:4044080/93‑1 (MQ=255)

CCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAT  >  NC_000913/1841834‑1842017

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