BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,713,746	C→T	P56P (CCG→CCA)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,713,746	0	C	T	99.5%	764.9 / NA	216	P56P (CCG→CCA)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (110/105); minor base G (1/0); total (111/105)
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 229 total aligned reads displayed.

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGC  >  NC_000913/1713656‑1713834

cgATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtgg                                                                                        <  1:2048110/93‑1 (MQ=255)

 gATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtggc                                                                                       >  1:4861791/1‑93 (MQ=255)

 gATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtggc                                                                                       >  2:294330/1‑93 (MQ=255)

 gATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtggc                                                                                       >  1:6237454/1‑93 (MQ=255)

   tCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCAt                                                                                     >  1:3170256/1‑93 (MQ=255)

    cccAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATc                                                                                    <  1:4118057/93‑1 (MQ=255)

     ccAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCa                                                                                   <  1:1404860/93‑1 (MQ=255)

      cAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCAc                                                                                  <  1:3441197/93‑1 (MQ=255)

       aGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACg                                                                                 <  2:1917276/93‑1 (MQ=255)

       aGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACg                                                                                 >  1:4049737/1‑93 (MQ=255)

        gATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGc                                                                                <  2:218239/93‑1 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                            >  2:416303/1‑93 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                            <  1:1633138/93‑1 (MQ=255)

             ccACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCa                                                                           <  2:5860230/93‑1 (MQ=255)

             ccACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCa                                                                           <  2:5490842/93‑1 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                        >  2:5342735/1‑93 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                        >  2:5268062/1‑93 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                        >  2:75561/1‑93 (MQ=255)

                 ttAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTcc                                                                       <  2:5059516/93‑1 (MQ=255)

                 ttAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTcc                                                                       <  1:4912960/93‑1 (MQ=255)

                  tAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCa                                                                      >  1:4206713/1‑93 (MQ=255)

                    aTACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAtt                                                                    >  1:986544/1‑93 (MQ=255)

                    aTACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAtt                                                                    <  2:261486/93‑1 (MQ=255)

                       cGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGc                                                                 <  2:5112640/93‑1 (MQ=255)

                        ggCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGcc                                                                <  1:5924185/93‑1 (MQ=255)

                         gCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCg                                                               >  2:5548006/1‑93 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGACAGTCCATTTGCCGTAt                                                            >  1:5063776/1‑93 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGACTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                            >  1:3885405/1‑93 (MQ=255)

                             cacaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                           >  1:3820037/1‑93 (MQ=255)

                             cacaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCACCACGCAGCCAGTCCATTTGCCGTAtt                                                           <  2:4418527/93‑1 (MQ=255)

                              acaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTg                                                          >  2:1048680/1‑93 (MQ=255)

                              acaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTg                                                          <  1:1000534/93‑1 (MQ=255)

                               caGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTgg                                                         >  2:1502389/1‑93 (MQ=255)

                               caGGTGTGTAATTTATCAATGGCGGCAATGCCATGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTgg                                                         <  2:5760853/93‑1 (MQ=255)

                                aGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTggg                                                        <  1:5364305/93‑1 (MQ=255)

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgt                                                       <  1:2223360/93‑1 (MQ=255)

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgt                                                       <  1:4408925/93‑1 (MQ=255)

                                  gtgtgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtg                                                      <  1:3314960/93‑1 (MQ=255)

                                   tgtgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtgt                                                     >  1:4640996/1‑93 (MQ=255)

                                   tgtgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtgt                                                     <  1:4168428/93‑1 (MQ=255)

                                       tAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGAtt                                                 >  2:717060/1‑93 (MQ=255)

                                        aaTTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTa                                                <  2:3750934/93‑1 (MQ=255)

                                        aaTTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTa                                                <  1:4296084/93‑1 (MQ=255)

                                           ttATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaa                                             >  1:5709705/1‑93 (MQ=255)

                                                 aTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGa                                       <  1:1962949/93‑1 (MQ=255)

                                                  tGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGaa                                      <  1:5970585/93‑1 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAAc                                     <  1:1849736/93‑1 (MQ=255)

                                                    gcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                    <  1:4743726/93‑1 (MQ=255)

                                                    gcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCACTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                    >  2:1240461/1‑93 (MQ=255)

                                                      ggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGGGGCATCACGCAGCCCGTCCATTTGACGTATTAGGTGTGATTAGAAAAGTGAACGgt                                  >  2:3272445/1‑93 (MQ=255)

                                                       gcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtg                                 >  1:3129399/1‑93 (MQ=255)

                                                        cAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtgt                                >  2:5718881/1‑93 (MQ=255)

                                                        cAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtgt                                <  2:2290128/93‑1 (MQ=255)

                                                        cAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtgt                                >  2:2914850/1‑93 (MQ=255)

                                                         aaTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGtt                               >  1:2414381/1‑93 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                              >  2:2377585/1‑93 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                              <  2:6066140/93‑1 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                              <  2:5960352/93‑1 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                              <  1:1088551/93‑1 (MQ=255)

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                             <  1:1288573/93‑1 (MQ=255)

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                             >  2:4513187/1‑93 (MQ=255)

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                             >  1:5583605/1‑93 (MQ=255)

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                             >  2:3014924/1‑93 (MQ=255)

                                                            gCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAg                            <  2:2734431/93‑1 (MQ=255)

                                                             ccTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGc                           <  2:2036617/93‑1 (MQ=255)

                                                             ccTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGc                           <  1:6127064/93‑1 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                          <  1:2598352/93‑1 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                          <  1:5724150/93‑1 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                          >  2:6193616/1‑93 (MQ=255)

                                                               ttGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCgg                         >  2:5021350/1‑93 (MQ=255)

                                                                 gCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGcc                       <  1:5107256/93‑1 (MQ=255)

                                                                 gCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGcc                       >  2:796647/1‑93 (MQ=255)

                                                                  cacaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCa                      <  2:1396908/93‑1 (MQ=255)

                                                                  cacaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCa                      >  2:1693863/1‑93 (MQ=255)

                                                                  cacaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGCGTTCAGCGGCCa                      >  2:2158976/1‑93 (MQ=255)

                                                                   acaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTAAGCGGCCAg                     >  1:1201895/1‑93 (MQ=255)

                                                                    caATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGa                    >  2:4026155/1‑93 (MQ=255)

                                                                         tCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc               <  2:4255689/93‑1 (MQ=255)

                                                                         tCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc               >  2:1629053/1‑93 (MQ=255)

                                                                           ggTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgc             >  1:1818986/1‑93 (MQ=255)

                                                                           ggTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgc             >  1:2293754/1‑93 (MQ=255)

                                                                           ggTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgc             <  1:2452308/93‑1 (MQ=255)

                                                                           ggTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgc             >  1:2694117/1‑93 (MQ=255)

                                                                           ggTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgc             >  1:279538/1‑93 (MQ=255)

                                                                            gTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcgcg            >  2:4039704/1‑93 (MQ=255)

                                                                             ttAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCgcgc           >  2:2667579/1‑93 (MQ=255)

                                                                              tAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGcc          >  1:2991063/1‑93 (MQ=255)

                                                                               aaaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGccc         >  1:4081856/1‑93 (MQ=255)

                                                                               aaaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGccc         <  2:5728354/93‑1 (MQ=255)

                                                                               aaaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGccc         <  1:5543591/93‑1 (MQ=255)

                                                                                aaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCa        <  2:3242327/93‑1 (MQ=255)

                                                                                 aTGGCTGGGGGGCATCACGCAGCCAGGCCATTTTTCGCATTTGGGGTGGTTAGAAAATTGAACGGGGTTCAGCGGGCAGACGGTCcccccccg       >  2:304699/1‑91 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc     <  1:3644775/93‑1 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc     <  2:5200771/93‑1 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc     <  1:1356133/93‑1 (MQ=255)

                                                                                    gCTGGGTGGTAACACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGATCGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcg    <  2:2369236/93‑1 (MQ=255)

                                                                                     cTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcgg   <  1:5217826/93‑1 (MQ=255)

                                                                                      tGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcggc  <  2:5574433/93‑1 (MQ=255)

                                                                                      tGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcggc  <  2:4755777/93‑1 (MQ=255)

                                                                                      tGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcggc  >  2:3453100/1‑93 (MQ=255)

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGC  >  NC_000913/1713656‑1713834

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