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.6%	882.6 / ‑5.5	263	P56P (CCG→CCA)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (160/101); total (161/102)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 270 total aligned reads displayed.

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

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

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

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

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

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

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

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

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

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

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

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

         aTAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCa                                                                                <  1:6150307/93‑1 (MQ=255)

         aTAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCa                                                                                >  1:5673117/1‑93 (MQ=255)

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

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

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

             ccACTTACTACGGCATCACAGGTGTGTAATTAATCCATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGACATCACGCAGCCa                                                                            <  2:8925794/93‑1 (MQ=255)

             ccACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCa                                                                            >  1:7786016/1‑93 (MQ=255)

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

               aCTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGt                                                                          >  1:4481545/1‑93 (MQ=255)

               aCTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGt                                                                          >  2:6424699/1‑93 (MQ=255)

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

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

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

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

                   aaTACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAt                                                                      >  1:4371663/1‑93 (MQ=255)

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

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

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

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

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

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

                          cATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGt                                                               >  2:6766480/1‑93 (MQ=255)

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

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

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

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTACATGGCTGGGTGGCATCACGCAGCCCGTCCATTTGCCGTAt                                                             >  1:1027740/1‑93 (MQ=255)

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

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                             <  2:501933/93‑1 (MQ=255)

                             cacaGGTTTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                            <  2:602945/93‑1 (MQ=255)

                             cacaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                            <  1:963779/93‑1 (MQ=255)

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

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

                               caGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTgg                                                          <  1:1442994/93‑1 (MQ=255)

                                aGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTggg                                                         >  1:9879094/1‑93 (MQ=255)

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTgtgt                                                        >  2:5119289/1‑93 (MQ=255)

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

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

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

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

                                     tgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGa                                                    >  1:6948173/1‑93 (MQ=255)

                                     tgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGa                                                    >  1:8250442/1‑93 (MQ=255)

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

                                       tAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGGGTGAtt                                                  >  2:7332911/1‑93 (MQ=255)

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

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

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

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

                                            tATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaaa                                             >  2:6528375/1‑93 (MQ=255)

                                            tATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaaa                                             >  2:5840279/1‑93 (MQ=255)

                                               cAATGGCGGCAATGCCTTGCACAATTTCGGTTACATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAtt                                          >  2:6368811/1‑93 (MQ=255)

                                               cAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAtt                                          >  1:4128306/1‑93 (MQ=255)

                                                aaTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCCTCACGCAGCCAGTCCCTTTGCCGTACTGGGTGTGATTAGAAAATTg                                         >  2:5206568/1‑93 (MQ=255)

                                                 aTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGTGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGa                                        <  1:8015094/93‑1 (MQ=255)

                                                 aTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGGGTGATTAGAAAATTGa                                        >  2:312396/1‑93 (MQ=255)

                                                  tGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTTGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGaa                                       <  1:3338553/93‑1 (MQ=255)

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

                                                       gcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtg                                  <  1:4662307/93‑1 (MQ=255)

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

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

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

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

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

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

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACACAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                              <  1:8355089/93‑1 (MQ=255)

                                                            gCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAg                             >  2:4122215/1‑93 (MQ=255)

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

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATCAGAAAATTGAACGGTGTTCAGCg                           >  2:5504417/1‑93 (MQ=255)

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

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

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

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

                                                                      aTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACg                   >  1:5379984/1‑93 (MQ=255)

                                                                      aTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACg                   >  2:5564372/1‑93 (MQ=255)

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

                                                                        ttCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGtt                 >  2:6441451/1‑93 (MQ=255)

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

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

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

                                                                           ggTTAAATGGCCGGGTGGCATCACGCAGCCATTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGATGTTCAGCGGCCAGACGTTcgc              <  2:8711236/93‑1 (MQ=255)

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

                                                                                 aTGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAg        <  1:5532279/93‑1 (MQ=255)

                                                                                  tGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAgg       <  1:3057564/93‑1 (MQ=255)

                                                                                  tGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAgg       >  1:5191696/1‑93 (MQ=255)

                                                                                  tGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAgg       <  2:2460895/93‑1 (MQ=255)

                                                                                  tGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAgg       <  1:1699095/93‑1 (MQ=255)

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

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

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

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

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

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