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%	1031.2 / ‑5.8	278	R245C (CGC→TGC)	intD	prophage DLP12 integrase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (139/138); total (140/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 289 total aligned reads displayed.

TGGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTTC  >  NC_000913/564379‑564558

tGGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGt                                                                                         >  1:5594172/1‑93 (MQ=255)

 tgTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGtt                                                                                        <  2:2049479/92‑1 (MQ=255)

 ggTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGtt                                                                                        <  1:4193902/93‑1 (MQ=255)

   tATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTg                                                                                      >  2:1019246/1‑93 (MQ=255)

    atatACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTGCAGGCGGTAGCGTTCAGCGCCACACCAATGGCGCAGTTTGa                                                                                     <  1:3042733/90‑1 (MQ=255)

    aTGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGa                                                                                     >  2:613621/1‑93 (MQ=255)

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

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

         cAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGc                                                                                <  2:1572547/93‑1 (MQ=255)

         cAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGc                                                                                <  2:5700011/93‑1 (MQ=255)

           aaCACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctc                                                                              >  2:2962505/1‑93 (MQ=255)

            acacCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctct                                                                             >  2:1257108/1‑93 (MQ=255)

              acCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTAGAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctct                                                                           >  2:4706933/1‑93 (MQ=255)

              acCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctct                                                                           >  1:942494/1‑93 (MQ=255)

                 cATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtcttct                                                                        >  1:2337803/1‑93 (MQ=255)

                 cATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtcttct                                                                        <  1:3816067/93‑1 (MQ=255)

                 cATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtcttct                                                                        >  2:3912768/1‑93 (MQ=255)

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

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

                    ttGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGa                                                                     <  1:5365789/93‑1 (MQ=255)

                     tgtgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGAt                                                                    <  1:82534/93‑1 (MQ=255)

                     tgtgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGAt                                                                    >  1:4258506/1‑93 (MQ=255)

                     tgtgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGAt                                                                    <  1:3270989/93‑1 (MQ=255)

                      gtgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGAtt                                                                   <  2:4803416/93‑1 (MQ=255)

                       tgATGCTTGCCTATTTGATCACGCACCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGAGCAGTCTGATTTGCTCTCTTCTGGATTc                                                                  >  1:3955381/1‑93 (MQ=255)

                       tgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTc                                                                  >  1:3504515/1‑93 (MQ=255)

                          tgtttGCTTATTGGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAcc                                                               <  1:4111255/90‑1 (MQ=255)

                          tGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTTTGGATTCAcc                                                               >  1:473937/1‑93 (MQ=255)

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

                           gCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAccc                                                              >  1:3440308/1‑93 (MQ=255)

                           gCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCAccc                                                              >  1:965631/1‑93 (MQ=255)

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

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

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

                               gCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCACGc                                                          >  1:54818/1‑93 (MQ=255)

                                ccTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCa                                                         >  1:5656142/1‑93 (MQ=255)

                                 cTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCaa                                                        >  2:86037/1‑93 (MQ=255)

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

                                  tATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAAc                                                       <  1:1983394/93‑1 (MQ=255)

                                  tATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAAc                                                       <  1:1756610/93‑1 (MQ=255)

                                  caTTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACCCCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAAc                                                       <  1:6638674/92‑1 (MQ=255)

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

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

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

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

                                             gCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCa                                            >  1:2096483/1‑93 (MQ=255)

                                              cAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCAt                                           <  1:602766/93‑1 (MQ=255)

                                              cAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCAt                                           >  1:3897089/1‑93 (MQ=255)

                                              cAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCAt                                           <  2:6350938/93‑1 (MQ=255)

                                               aaCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATg                                          >  1:5115904/1‑93 (MQ=255)

                                               aaCACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATg                                          >  1:491603/1‑93 (MQ=255)

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

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

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

                                                 cacTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTc                                        <  2:4053961/93‑1 (MQ=255)

                                                  acTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCg                                       >  1:4437197/1‑93 (MQ=255)

                                                  acTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCg                                       >  2:1431524/1‑93 (MQ=255)

                                                    tttACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAt                                     >  2:1975549/1‑93 (MQ=255)

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

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

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

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

                                                         aGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATttgtt                                >  2:972632/1‑93 (MQ=255)

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

                                                          ggCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATttgttg                               >  2:4465505/1‑93 (MQ=255)

                                                            cGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTTGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGcc                             <  1:2213046/93‑1 (MQ=255)

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

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

                                                                aTCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTc                         >  1:2021065/1‑93 (MQ=255)

                                                                aTCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTc                         >  1:1494060/1‑93 (MQ=255)

                                                                 tCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTcc                        >  2:5423705/1‑93 (MQ=255)

                                                                 tCGTTCAGCGCCACACCAAAGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGATGCCAtttc                        >  1:1875995/1‑91 (MQ=255)

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

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

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

                                                                      cAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGAtt                   <  1:3647940/93‑1 (MQ=255)

                                                                        gcgcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTa                 <  2:2627808/93‑1 (MQ=255)

                                                                         cgcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTat                <  1:672859/93‑1 (MQ=255)

                                                                         cgcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTat                <  2:1488173/93‑1 (MQ=255)

                                                                           ccACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatga              >  1:1383541/1‑93 (MQ=255)

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

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

                                                                            cacaCCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatgat             <  1:6235391/93‑1 (MQ=255)

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

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

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

                                                                                ccAATGACGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTc         <  1:1013568/93‑1 (MQ=255)

                                                                                 cAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCg        >  1:504833/1‑93 (MQ=255)

                                                                                  aaTGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGa       >  2:685991/1‑93 (MQ=255)

                                                                                   aTGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGAc      <  2:1210022/93‑1 (MQ=255)

                                                                                   aTGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGAc      <  2:5355254/93‑1 (MQ=255)

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

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

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

                                                                                    tGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACt     <  1:5452209/93‑1 (MQ=255)

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

                                                                                      gcgcAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACttt   <  1:1085267/93‑1 (MQ=255)

                                                                                       cgcaGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTtc  <  2:6565969/93‑1 (MQ=255)

                                                                                       cgcaGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTtc  <  2:4344090/93‑1 (MQ=255)

                                                                                       cgcaGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTtc  <  1:2622348/93‑1 (MQ=255)

TGGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTTC  >  NC_000913/564379‑564558

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