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.5%	676.2 / NA	187	R245C (CGC→TGC)	intD	prophage DLP12 integrase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (105/81); minor base C (1/0); total (106/81)
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 192 total aligned reads displayed.

TGGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTT  >  NC_000913/564379‑564557

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

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

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

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

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

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

       tACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGAttt                                                                                 <  1:1629923/93‑1 (MQ=255)

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

        aCAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTg                                                                                <  1:1583169/93‑1 (MQ=255)

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

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

           aaCACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctc                                                                             <  2:845957/93‑1 (MQ=255)

           aaCACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctc                                                                             <  1:4311790/93‑1 (MQ=255)

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

             cccccATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctc                                                                           <  1:665922/91‑1 (MQ=255)

             cacCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctc                                                                           <  1:1080524/93‑1 (MQ=255)

              acCCATTCGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGctctct                                                                          >  2:3021779/1‑93 (MQ=255)

               cccATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtctt                                                                         <  2:3195645/93‑1 (MQ=255)

               cccATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtctt                                                                         >  1:741687/1‑93 (MQ=255)

                ccATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCtcttc                                                                        <  1:562770/93‑1 (MQ=255)

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

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

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

                  aTTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTg                                                                      >  1:2239717/1‑93 (MQ=255)

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

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

                       tgATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTc                                                                 <  1:2673315/93‑1 (MQ=255)

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

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

                             ttGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGTCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAg                                                           <  2:2590055/93‑1 (MQ=255)

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

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

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

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

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

                                ccTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGAGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTAACCCAGGCa                                                        >  1:1449939/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

                                     ttGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcg                                                   >  1:904363/1‑93 (MQ=255)

                                       gATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACtcgtc                                                 <  2:1749419/93‑1 (MQ=255)

                                          cACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCt                                              >  2:3348474/1‑93 (MQ=255)

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

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

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

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

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

                                                acacTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCGCTCTCCTGGATTCACCCAGGCAACTCGTCGCTGCATGt                                        <  2:1235527/93‑1 (MQ=255)

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

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

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

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

                                                      tACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAttt                                  <  2:1804555/93‑1 (MQ=255)

                                                      tACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAttt                                  <  1:101677/93‑1 (MQ=255)

                                                      tACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGAttt                                  <  1:3576073/93‑1 (MQ=255)

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

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

                                                           gCGGTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGc                             >  2:3770505/1‑93 (MQ=255)

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

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

                                                              gTATCGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCAt                          >  1:2951433/1‑93 (MQ=255)

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

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

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

                                                                  cGTTCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCa                      <  1:2995120/93‑1 (MQ=255)

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

                                                                    gtCAGCGCCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGa                    <  1:3690141/92‑1 (MQ=255)

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

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

                                                                        gcgcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTa                >  1:1369209/1‑93 (MQ=255)

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

                                                                          gcCACACCAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTatg              <  1:2932617/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                ccAATGGCGCAGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTc        <  2:1282865/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                      gcgcAGCTTGATTTGCCCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACttt  >  2:2125936/1‑93 (MQ=255)

TGGTATGTACAAACACCCATTTGTGATGCTTGCCTATTTGATCACGCAACACTTTACAGGCGGTATCGTTCAGCGCCACACCAATGGCGCGGTTTGATTTGCTCTCTTCTGGATTCACCCAGGCAACTCGTCGCTGCATGTCGATTTGTTGCCATTCCAGATTTATGATGTTCGACTTT  >  NC_000913/564379‑564557

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