BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	564,469	G→A	R245C (CGC→TGC)	intD ←	prophage DLP12 integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	564,469	0	G	A	99.6%	1008.0 / NA	270	R245C (CGC→TGC)	intD	prophage DLP12 integrase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (133/136); minor base C (1/0); total (134/136)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTTC  >  NC_000913/564380‑564558

 gTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGttt                                                                                       >  2:5876073/1‑93 (MQ=255)

 gTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGttt                                                                                       <  1:1881729/93‑1 (MQ=255)

  tATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGGGCAGTTTg                                                                                      <  2:2334234/93‑1 (MQ=255)

  tATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTg                                                                                      <  1:1451181/93‑1 (MQ=255)

   aTGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGa                                                                                     <  1:102707/93‑1 (MQ=255)

    tGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGAt                                                                                    >  1:4790689/1‑93 (MQ=255)

     gTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGAtt                                                                                   <  1:551500/93‑1 (MQ=255)

       aCAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTg                                                                                 >  2:5788528/1‑93 (MQ=255)

        cAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGc                                                                                >  2:3125468/1‑93 (MQ=255)

         taaCCCCAATTTGAGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCATCGCCACACCAATGGCGCAGTTTGATTTGct                                                                               >  2:111776/2‑93 (MQ=255)

         aaaCACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGct                                                                               >  1:465128/1‑93 (MQ=255)

           acacCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctct                                                                             <  2:2852860/93‑1 (MQ=255)

           acacCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctct                                                                             <  2:4638770/93‑1 (MQ=255)

             acCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctct                                                                           <  1:987097/93‑1 (MQ=255)

                 aTTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTg                                                                       <  2:2226422/93‑1 (MQ=255)

                 aTTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTg                                                                       <  1:5859033/93‑1 (MQ=255)

                  tttGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTgg                                                                      >  2:4476349/1‑93 (MQ=255)

                  tttGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTgg                                                                      >  1:3238436/1‑93 (MQ=255)

                   ttGTGATGCTTGCCTATTTGATCACGCAACACTTAACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTgaa                                                                     <  2:4273930/93‑3 (MQ=255)

                     gtgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGAtt                                                                   >  1:1690529/1‑93 (MQ=255)

                       gATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTcc                                                                 >  1:5119752/1‑92 (MQ=255)

                       gATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCa                                                                 <  2:160211/93‑1 (MQ=255)

                       gATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCa                                                                 >  2:537708/1‑93 (MQ=255)

                       gATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCa                                                                 >  1:3040921/1‑93 (MQ=255)

                        aTGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAc                                                                >  2:1994377/1‑93 (MQ=255)

                        aTGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAc                                                                >  1:710424/1‑93 (MQ=255)

                        aTGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAc                                                                <  1:4046265/93‑1 (MQ=255)

                         tGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAcc                                                               <  2:249489/93‑1 (MQ=255)

                         tGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAcc                                                               >  2:3948299/1‑93 (MQ=255)

                          gCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAccc                                                              <  1:5273235/93‑1 (MQ=255)

                          gCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAccc                                                              <  2:5192908/93‑1 (MQ=255)

                           cTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCa                                                             <  2:3375024/93‑1 (MQ=255)

                           cTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCa                                                             >  2:5675348/1‑93 (MQ=255)

                           cTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCa                                                             <  1:5854077/93‑1 (MQ=255)

                            ttGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAg                                                            <  2:5502656/93‑1 (MQ=255)

                            ttGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAg                                                            >  2:6044177/1‑93 (MQ=255)

                            ttGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAg                                                            <  1:3899018/93‑1 (MQ=255)

                              gCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAgga                                                          >  2:2404222/1‑92 (MQ=255)

                              gCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGc                                                          >  2:1906033/1‑93 (MQ=255)

                              gCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGc                                                          >  1:3496571/1‑93 (MQ=255)

                                 tATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCCGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAAc                                                       >  1:1195689/1‑93 (MQ=255)

                                 tATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAAc                                                       >  1:4836274/1‑93 (MQ=255)

                                  aTTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACt                                                      <  1:4242080/93‑1 (MQ=255)

                                  aTTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACt                                                      >  2:1448851/1‑93 (MQ=255)

                                   tttGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtc                                                     <  1:4132015/93‑1 (MQ=255)

                                    ttGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcg                                                    <  1:3691601/93‑1 (MQ=255)

                                     tGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcgt                                                   <  1:4580143/93‑1 (MQ=255)

                                      gATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcgtc                                                  >  1:6233460/1‑93 (MQ=255)

                                       aTCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcgtcg                                                 <  2:4929454/93‑1 (MQ=255)

                                        tCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGc                                                <  2:5718801/93‑1 (MQ=255)

                                          aCGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTg                                              >  2:3763789/1‑93 (MQ=255)

                                          aCGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTg                                              >  2:1820696/1‑93 (MQ=255)

                                            gCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCa                                            <  1:2068191/93‑1 (MQ=255)

                                              aaCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATg                                          <  1:2791069/93‑1 (MQ=255)

                                              aaCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATg                                          <  2:856945/93‑1 (MQ=255)

                                              aaCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATg                                          <  1:3476238/93‑1 (MQ=255)

                                               acacTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGt                                         >  1:4522607/1‑93 (MQ=255)

                                               acacTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGt                                         <  1:5020122/93‑1 (MQ=255)

                                               acacTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGt                                         >  1:2334530/1‑93 (MQ=255)

                                                  cTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGa                                      >  1:4916996/1‑93 (MQ=255)

                                                    ttACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAtt                                    <  1:365925/93‑1 (MQ=255)

                                                     tACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAttt                                   >  1:5782776/1‑93 (MQ=255)

                                                     tACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAttt                                   >  2:663675/1‑93 (MQ=255)

                                                       cAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATttgt                                 >  1:5068821/1‑93 (MQ=255)

                                                       cAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATttgt                                 <  1:6099877/93‑1 (MQ=255)

                                                        aGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATttgtt                                <  2:4129144/93‑1 (MQ=255)

                                                          gCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGc                              <  1:3345640/93‑1 (MQ=255)

                                                          gCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGc                              <  1:2582474/93‑1 (MQ=255)

                                                            ggTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCa                            >  1:5462765/1‑93 (MQ=255)

                                                             gTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCAt                           <  1:4377442/93‑1 (MQ=255)

                                                             gTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACACAGGCAACTCGTCGCTGCATGTCGATTTGTTGACAt                           >  1:2883006/1‑93 (MQ=255)

                                                                tCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTcc                        <  2:4691162/93‑1 (MQ=255)

                                                                  gTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAg                      >  1:3177876/1‑93 (MQ=255)

                                                                   ttCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGa                     <  1:4425307/93‑1 (MQ=255)

                                                                    tCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGAt                    <  2:3346475/93‑1 (MQ=255)

                                                                     taGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGAtt                   <  1:3775366/92‑1 (MQ=255)

                                                                      aGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGAttt                  >  1:3426260/1‑93 (MQ=255)

                                                                      aGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGAttt                  <  1:530401/93‑1 (MQ=255)

                                                                       gcgcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTCGCCATTCCAGATTTa                 <  2:352356/93‑1 (MQ=255)

                                                                        ggcCACACCAATGGCGCAGTTTGATTTGATCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTat                <  1:3114401/92‑1 (MQ=255)

                                                                         gcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatg               >  1:224787/1‑93 (MQ=255)

                                                                         gcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatg               >  1:670826/1‑93 (MQ=255)

                                                                          ccACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatga              <  2:338693/93‑1 (MQ=255)

                                                                          ccACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatga              <  2:1228180/93‑1 (MQ=255)

                                                                           cacaCCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgat             >  2:3879407/1‑93 (MQ=255)

                                                                           cacaCCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgat             >  2:5353965/1‑93 (MQ=255)

                                                                            ccacCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgatg            <  1:2500282/92‑1 (MQ=255)

                                                                            acacCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgatg            >  1:906494/1‑93 (MQ=255)

                                                                            acacCAACGGCGCAGTGTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgatg            <  1:1177270/93‑1 (MQ=255)

                                                                             cacCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGt           >  2:107397/1‑93 (MQ=255)

                                                                             cacCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGt           <  2:6172393/93‑1 (MQ=255)

                                                                              acCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGtt          <  1:1812823/93‑1 (MQ=255)

                                                                               ccAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTc         >  2:4596440/1‑93 (MQ=255)

                                                                               ccAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTc         >  1:6142846/1‑93 (MQ=255)

                                                                               ccAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTc         >  1:3618181/1‑93 (MQ=255)

                                                                                 aaTGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGa       <  2:4008170/93‑1 (MQ=255)

                                                                                  aTGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGAc      >  1:1971036/1‑93 (MQ=255)

                                                                                   tGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACt     >  2:4655174/1‑93 (MQ=255)

                                                                                    ggCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACtt    >  2:2994642/1‑93 (MQ=255)

                                                                                    ggCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACtt    <  2:2652998/93‑1 (MQ=255)

GGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTTC  >  NC_000913/564380‑564558

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