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%	648.7 / NA	192	P56P (CCG→CCA)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (109/82); minor base G (1/0); total (110/82)
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 197 total aligned reads displayed.

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

cgATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtgg                                                                                         >  1:7291551/1‑93 (MQ=255)

  aTCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCa                                                                                       >  1:1697925/1‑93 (MQ=255)

  aTCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCa                                                                                       >  1:1271501/1‑93 (MQ=255)

    cccAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATc                                                                                     >  2:4926643/1‑93 (MQ=255)

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

    cccAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATc                                                                                     >  1:7055410/1‑93 (MQ=255)

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

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

     ccAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCa                                                                                    >  2:6293353/1‑93 (MQ=255)

     ccACATAGCCACTTAATACGGCCTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCa                                                                                    <  1:2902324/93‑1 (MQ=255)

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

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

        gATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGc                                                                                 >  1:4828687/1‑93 (MQ=255)

          tAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAg                                                                               >  1:3321063/1‑93 (MQ=255)

           aGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGc                                                                              <  2:4736435/93‑1 (MQ=255)

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

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

            gCCACCTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                             >  2:7856843/1‑93 (MQ=255)

              cACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAg                                                                           <  1:2032915/93‑1 (MQ=255)

              cACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAg                                                                           >  1:6671271/1‑93 (MQ=255)

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

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                         <  2:7871895/93‑1 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                         <  2:3546440/93‑1 (MQ=255)

                 ttAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTcc                                                                        >  2:5708323/1‑93 (MQ=255)

                  tAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCa                                                                       <  1:1877726/93‑1 (MQ=255)

                  tAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCa                                                                       <  2:2300679/93‑1 (MQ=255)

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

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

                     tACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAttt                                                                    >  2:749394/1‑93 (MQ=255)

                       cGGCATCACAGGTGTGTAATTTATCAATGGCGGCAAGGCCTTGCACAATTTCGGTTCACTGGCTGGGTGGAACCACACGACCCGTCAATTTGc                                                                  >  1:6384534/1‑93 (MQ=255)

                        ggCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGcc                                                                 >  2:6621811/1‑93 (MQ=255)

                           aTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTa                                                              >  1:4199856/1‑93 (MQ=255)

                           aTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTa                                                              >  2:6676830/1‑93 (MQ=255)

                           aTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTa                                                              >  2:5462611/1‑93 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                             >  1:6844476/1‑93 (MQ=255)

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

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

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

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgt                                                        >  2:7067754/1‑93 (MQ=255)

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgt                                                        >  1:371361/1‑93 (MQ=255)

                                  gtgtgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtg                                                       >  2:6952431/1‑93 (MQ=255)

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

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

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

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

                                      gtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGAt                                                   >  2:6675328/1‑93 (MQ=255)

                                        aaTTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTa                                                 >  1:1604766/1‑93 (MQ=255)

                                          tttATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGa                                               >  2:5596070/1‑93 (MQ=255)

                                             aTCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaaaa                                            <  2:2785708/93‑1 (MQ=255)

                                              tCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAt                                           <  2:7799420/93‑1 (MQ=255)

                                              tCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAt                                           >  1:1130712/1‑93 (MQ=255)

                                              tCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATCTGCCGTATTGGGTGTGATTAGAAAAt                                           >  1:582099/1‑93 (MQ=255)

                                                aaTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTg                                         >  2:4846026/1‑93 (MQ=255)

                                                 aTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGa                                        >  1:4755828/1‑93 (MQ=255)

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

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

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

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

                                                    gcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                     >  2:5922666/1‑93 (MQ=255)

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

                                                     cggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACgg                                    <  2:7649009/93‑1 (MQ=255)

                                                     cggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACgg                                    >  2:7832255/1‑93 (MQ=255)

                                                      ggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgt                                   >  1:7481448/1‑93 (MQ=255)

                                                      ggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgt                                   <  2:196286/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                               ttGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCgg                          <  2:7378593/93‑1 (MQ=255)

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

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

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

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

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

                                                                   acaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAg                      >  2:1936283/1‑93 (MQ=255)

                                                                   acaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAg                      >  2:366290/1‑93 (MQ=255)

                                                                    caATTTCGGTTAAATGGCTGGGTGGCATAACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGa                     <  1:7761927/93‑1 (MQ=255)

                                                                     aaTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGAc                    >  1:6148500/1‑93 (MQ=255)

                                                                     aaTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGAc                    <  1:5899485/93‑1 (MQ=255)

                                                                       tttCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGt                  >  1:5990533/1‑93 (MQ=255)

                                                                        ttCGGTTAAATTGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGtt                 <  1:2649120/93‑1 (MQ=255)

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

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

                                                                          cGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcg               <  2:1528893/93‑1 (MQ=255)

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

                                                                             ttAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCgcgc            <  1:4302615/93‑1 (MQ=255)

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

                                                                              tAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGCGTTCAGCGGCCAGACGTTCGCGcc           >  2:5553801/1‑93 (MQ=255)

                                                                               aaaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGGACGGTGTTCAGCGGCCAGACGTTCGCGccc          >  2:4150482/1‑93 (MQ=255)

                                                                                aaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCa         >  1:5558684/1‑93 (MQ=255)

                                                                                aaTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAACTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCa         <  1:3004845/93‑1 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc      >  2:5963681/1‑93 (MQ=255)

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

                                                                                       gggTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGgcggcg  >  1:5660036/1‑93 (MQ=255)

                                                                                       gggTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGgcggcg  <  2:1043124/93‑1 (MQ=255)

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

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