BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	510,451	C→T	Q51Q (CAG→CAA)	ybaR ←	putative ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	510,451	0	C	T	99.6%	943.3 / NA	255	Q51Q (CAG→CAA)	ybaR	putative ATPase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (139/115); total (140/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

TTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTT  >  NC_000913/510361‑510540

ttCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTg                                                                                         <  2:3583594/93‑1 (MQ=255)

 tCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGt                                                                                        <  1:5514039/93‑1 (MQ=255)

 tCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGt                                                                                        <  1:1015628/93‑1 (MQ=255)

   cGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTc                                                                                      <  1:3204158/93‑1 (MQ=255)

   cGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTc                                                                                      <  1:5784399/93‑1 (MQ=255)

   cGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTc                                                                                      >  2:231233/1‑93 (MQ=255)

    gACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCt                                                                                     >  2:4461973/1‑93 (MQ=255)

       gggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCTCTTGTTTGATGGTTTCGATCAGTTGTTCTGCa                                                                                  <  1:7022/93‑1 (MQ=255)

       gggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCa                                                                                  >  2:1248521/1‑93 (MQ=255)

       gggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCa                                                                                  >  2:652686/1‑93 (MQ=255)

        ggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACTGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCAc                                                                                 >  1:2706877/1‑93 (MQ=255)

        ggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCAc                                                                                 >  2:437536/1‑93 (MQ=255)

        ggATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCAc                                                                                 <  1:3951744/93‑1 (MQ=255)

         gATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTTGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACt                                                                                <  1:2036480/93‑1 (MQ=255)

         gATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACt                                                                                <  1:3518126/93‑1 (MQ=255)

         gATTGATGACTCCGCCAGCGGTTTAGCCTTGGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACt                                                                                <  2:3101029/93‑1 (MQ=255)

           ttGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTgg                                                                              <  2:3985616/93‑1 (MQ=255)

             gatgaCTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCa                                                                            <  2:2238119/93‑1 (MQ=255)

               tgaCTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGt                                                                          <  2:4407939/93‑1 (MQ=255)

                gaCTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTc                                                                         <  1:844791/93‑1 (MQ=255)

                  cTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTccc                                                                       >  1:1022515/1‑93 (MQ=255)

                  cTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTccc                                                                       >  2:28900/1‑93 (MQ=255)

                  cTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTccc                                                                       >  1:6067874/1‑93 (MQ=255)

                  cTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTccc                                                                       <  1:4872622/93‑1 (MQ=255)

                    ccGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCgg                                                                     <  2:652014/93‑1 (MQ=255)

                       ccNGCGGTNTNGCCTTTGGGTGGCTTACNGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTaa                                                                  <  2:4025573/93‑1 (MQ=255)

                          gCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGt                                                               <  2:5234462/93‑1 (MQ=255)

                           cGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTg                                                              <  2:3698184/93‑1 (MQ=255)

                            ggTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTgc                                                             >  2:2138043/1‑93 (MQ=255)

                             gTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTgcg                                                            >  1:4388559/1‑93 (MQ=255)

                              tttAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTgcgc                                                           >  1:2555196/1‑93 (MQ=255)

                               ttAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCt                                                          <  1:5388096/93‑1 (MQ=255)

                               ttAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCt                                                          <  1:5108748/93‑1 (MQ=255)

                                tAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCtt                                                         >  2:1804153/1‑93 (MQ=255)

                                 aGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTc                                                        >  2:4290456/1‑93 (MQ=255)

                                 aGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTc                                                        >  2:1701897/1‑93 (MQ=255)

                                 aGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTc                                                        >  1:451858/1‑93 (MQ=255)

                                    cTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGt                                                     >  2:3046718/1‑93 (MQ=255)

                                    cTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGt                                                     >  1:3916349/1‑93 (MQ=255)

                                     tttGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTg                                                    >  1:5644120/1‑93 (MQ=255)

                                      ttGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGa                                                   <  1:4664074/93‑1 (MQ=255)

                                      ttGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGa                                                   <  1:4826528/93‑1 (MQ=255)

                                      ttGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGa                                                   <  2:3776809/93‑1 (MQ=255)

                                        gggTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATa                                                 <  1:2928350/93‑1 (MQ=255)

                                        gggTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATa                                                 >  1:587719/1‑93 (MQ=255)

                                         ggTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAg                                                <  2:5302512/93‑1 (MQ=255)

                                         ggTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAg                                                >  1:3707450/1‑93 (MQ=255)

                                          gTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGa                                               >  2:3838392/1‑93 (MQ=255)

                                          gTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGa                                               <  2:1920128/93‑1 (MQ=255)

                                           tGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGac                                              >  1:2450031/1‑93 (MQ=255)

                                            ggCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGaca                                             <  2:2810599/93‑1 (MQ=255)

                                            ggCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGGCACGCGCGCTTCAGTGATAGaca                                             >  2:5326818/1‑93 (MQ=255)

                                             gCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGacac                                            <  2:2230593/93‑1 (MQ=255)

                                              cTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGAcaca                                           >  2:3320253/1‑93 (MQ=255)

                                               ttACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACAt                                          >  1:1529112/1‑93 (MQ=255)

                                                tACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATc                                         <  1:1131886/93‑1 (MQ=255)

                                                tACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATc                                         >  2:487310/1‑93 (MQ=255)

                                                     aTGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCt                                    <  2:5754553/93‑1 (MQ=255)

                                                       gCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGc                                  >  1:1371722/1‑93 (MQ=255)

                                                        cGTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCt                                 >  2:5416435/1‑93 (MQ=255)

                                                         gTCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTc                                >  2:450236/1‑93 (MQ=255)

                                                          tCATAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCa                               <  1:3032850/93‑1 (MQ=255)

                                                            aTAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAAc                             <  1:3320253/93‑1 (MQ=255)

                                                            aTAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAAc                             >  1:1088697/1‑93 (MQ=255)

                                                             tAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACa                            <  1:158998/93‑1 (MQ=255)

                                                             tAACCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACa                            >  2:3446416/1‑93 (MQ=255)

                                                               aCCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATc                          >  2:264110/1‑93 (MQ=255)

                                                               aCCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATc                          >  2:2707795/1‑93 (MQ=255)

                                                               aCCCGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATc                          <  1:1175396/93‑1 (MQ=255)

                                                                  cGCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCgg                       <  1:4272073/93‑1 (MQ=255)

                                                                   gCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGa                      <  1:4347030/93‑1 (MQ=255)

                                                                   gCTTGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGa                      >  1:2076592/1‑93 (MQ=255)

                                                                     ttGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACg                    >  2:622729/1‑93 (MQ=255)

                                                                     ttGTTTGATGGTTTCCATCAGTTGTTCTGCACTGGGAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACg                    <  1:1770669/93‑1 (MQ=255)

                                                                      tGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGc                   >  1:1309260/1‑93 (MQ=255)

                                                                      tGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGc                   >  1:2093946/1‑93 (MQ=255)

                                                                      tGTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGc                   >  2:5970947/1‑93 (MQ=255)

                                                                       gTTTGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCt                  >  2:2533171/1‑93 (MQ=255)

                                                                        tttGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTg                 >  1:6211697/1‑93 (MQ=255)

                                                                        tttGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTg                 <  2:4696612/93‑1 (MQ=255)

                                                                        tttGATGGTTTCGATCACTTGTTCTGCACCGCCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTg                 <  2:165297/93‑1 (MQ=255)

                                                                          tGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGtt               >  2:546396/1‑93 (MQ=255)

                                                                          tGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGtt               <  2:758327/93‑1 (MQ=255)

                                                                          tGATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGtt               >  2:2436104/1‑93 (MQ=255)

                                                                           gATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTc              >  1:5282532/1‑93 (MQ=255)

                                                                           gATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTc              >  1:4709045/1‑93 (MQ=255)

                                                                           gATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTc              >  1:2772384/1‑93 (MQ=255)

                                                                           gATGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTc              >  2:6211508/1‑93 (MQ=255)

                                                                            aTGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCa             >  2:1871668/1‑93 (MQ=255)

                                                                            aTGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCa             >  2:1043586/1‑93 (MQ=255)

                                                                            aTGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCa             <  1:4165306/93‑1 (MQ=255)

                                                                             tGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCaa            >  2:5001322/1‑93 (MQ=255)

                                                                             tGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCaa            >  2:1467425/1‑93 (MQ=255)

                                                                             tGGTTTCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCaa            <  1:2866029/93‑1 (MQ=255)

                                                                                tttCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGCTAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGAc         <  1:3375497/93‑1 (MQ=255)

                                                                                tttCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGAc         >  1:1949486/1‑93 (MQ=255)

                                                                                 ttCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGAct        <  2:5667170/93‑1 (MQ=255)

                                                                                  tCGATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActt       >  1:3786649/1‑93 (MQ=255)

                                                                                    gATCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttc     >  1:460002/1‑93 (MQ=255)

                                                                                      tCAGTTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttctt   <  2:182399/93‑1 (MQ=255)

TTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTT  >  NC_000913/510361‑510540

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