BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,466,831	T→C	I96T (ATT→ACT)	intS →	CPS‑53 (KpLE1) prophage; prophage CPS‑53 integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,466,831	0	T	C	100.0%	1226.8 / NA	280	I96T (ATT→ACT)	intS	CPS‑53 (KpLE1) prophage; prophage CPS‑53 integrase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (153/127); total (153/127)

Only 100 of 285 total aligned reads displayed.

AAGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGATTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCATT  >  NC_000913/2466733‑2466928

aaGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGActt                                                                                                 >  1:3974891/1‑101 (MQ=255)

 aGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTa                                                                                                >  2:3974893/1‑101 (MQ=255)

     aTAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAg                                                                                            >  2:3974899/1‑101 (MQ=255)

     aTAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAg                                                                                            >  1:3974898/1‑101 (MQ=255)

      tAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAgg                                                                                           >  1:3974903/1‑101 (MQ=255)

        aaGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGaa                                                                                         >  2:3974907/1‑101 (MQ=255)

         aGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAAt                                                                                        <  2:3974569/101‑1 (MQ=255)

          gCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATg                                                                                       <  1:3974340/101‑1 (MQ=255)

          gCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATg                                                                                       <  1:3974454/101‑1 (MQ=255)

          gCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATg                                                                                       >  2:3974909/1‑101 (MQ=255)

          gCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATg                                                                                       <  1:3974574/101‑1 (MQ=255)

            aTGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGt                                                                                     <  1:3974426/101‑1 (MQ=255)

            aTGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGt                                                                                     >  1:3974911/1‑101 (MQ=255)

              ggACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTAc                                                                                   <  1:3974521/101‑1 (MQ=255)

              ggACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTAc                                                                                   <  1:3974766/101‑1 (MQ=255)

               gACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACg                                                                                  >  2:3974915/1‑101 (MQ=255)

                aCTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGa                                                                                 >  1:3974916/1‑101 (MQ=255)

                 cTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAg                                                                                <  2:3974702/101‑1 (MQ=255)

                     gAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATgg                                                                                      >  2:3974924/1‑91 (MQ=255)

                     gAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGcaca                                                                            >  1:3974922/1‑101 (MQ=255)

                      aGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACaa                                                                           >  1:3974925/1‑101 (MQ=255)

                        gAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACaaga                                                                         >  2:3974928/1‑101 (MQ=255)

                          aaGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACaagaag                                                                       >  2:3974932/1‑101 (MQ=255)

                          aaGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACaaga                                                                         >  2:3974935/1‑99 (MQ=255)

                              cATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAg                                                                   <  1:3974611/101‑1 (MQ=255)

                               aTCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGt                                                                  <  1:3974480/101‑1 (MQ=255)

                                  tcGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATg                                                               <  1:3974376/101‑1 (MQ=255)

                                  tcGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATg                                                               >  1:3974940/1‑101 (MQ=255)

                                   cGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATgg                                                                                      >  2:3974941/1‑77 (MQ=255)

                                   cGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATgg                                                              <  1:3974643/101‑1 (MQ=255)

                                     gTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTc                                                            >  2:3974947/1‑101 (MQ=255)

                                     gTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTc                                                            >  1:3974945/1‑101 (MQ=255)

                                      ttGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCa                                                           <  1:3974608/101‑1 (MQ=255)

                                       tGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAg                                                          >  1:3974950/1‑101 (MQ=255)

                                        ggCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAATAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGt                                                         <  1:3974540/101‑1 (MQ=255)

                                         gCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATgg                                                              <  2:3974483/95‑1 (MQ=255)

                                         gCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTa                                                        >  2:3974955/1‑101 (MQ=255)

                                         gCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTa                                                        >  1:3974954/1‑101 (MQ=255)

                                           aTCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAgg                                                      >  2:3974960/1‑101 (MQ=255)

                                            tCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAggg                                                     >  1:3974961/1‑101 (MQ=255)

                                              gATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGc                                                                      >  1:3974966/1‑82 (MQ=255)

                                              gATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTa                                                   >  1:3974964/1‑101 (MQ=255)

                                                tCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATg                                                 >  2:3974970/1‑101 (MQ=255)

                                                  cTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGGAGCAAGTATGGTCAGTAGGGTATGCa                                               >  1:3974974/1‑101 (MQ=255)

                                                  cTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCa                                               <  1:3974504/101‑1 (MQ=255)

                                                   tGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCaa                                              <  1:3974598/101‑1 (MQ=255)

                                                    gTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAAc                                             >  2:3974975/1‑101 (MQ=255)

                                                    gTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAAc                                             <  1:3974528/101‑1 (MQ=255)

                                                     tAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACt                                            <  1:3974618/101‑1 (MQ=255)

                                                     tAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACt                                            <  1:3974808/101‑1 (MQ=255)

                                                        aGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGaa                                         >  2:3974981/1‑101 (MQ=255)

                                                          gCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAg                                                                   <  1:3974471/73‑1 (MQ=255)

                                                           cGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACtt                                      <  1:3974678/101‑1 (MQ=255)

                                                           cGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACtt                                      >  1:3974985/1‑101 (MQ=255)

                                                           cGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACtt                                      >  2:3974983/1‑101 (MQ=255)

                                                            gAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTg                                     <  2:3974722/101‑1 (MQ=255)

                                                            gAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTg                                     <  2:3974658/101‑1 (MQ=255)

                                                            gAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTg                                     >  2:3974987/1‑101 (MQ=255)

                                                             aaaaaGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCa                                               >  2:3974991/1‑90 (MQ=255)

                                                             aaaaaGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGc                                    <  2:3974559/101‑1 (MQ=255)

                                                               aaaGGCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCa                                  <  2:3974606/101‑1 (MQ=255)

                                                                  ggCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCaaaa                               >  2:3974994/1‑101 (MQ=255)

                                                                   gCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAAt                              >  1:3974998/1‑101 (MQ=255)

                                                                   gCTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAAt                              <  2:3974088/101‑1 (MQ=255)

                                                                   gCTTCGTCTAACAACAATTCCTTTAGTGCCACTTACAAGGAATGGTACGAGCACaagaa                                                                        <  2:3974516/59‑1 (MQ=255)

                                                                    cTTCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATg                             >  2:3974999/1‑101 (MQ=255)

                                                                     ttCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGt                            <  2:3974624/101‑1 (MQ=255)

                                                                     ttCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGt                            >  2:3975002/1‑101 (MQ=255)

                                                                     ttCGTCTAACAACAATTCCTTTAGTGCAACTTACAAGGAATGGTACGAGca                                                                              <  2:3974343/51‑1 (MQ=255)

                                                                      tCGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGtt                           <  2:3974533/101‑1 (MQ=255)

                                                                       cGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGttt                          <  1:3974705/101‑1 (MQ=255)

                                                                       cGTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGttt                          <  1:3974656/101‑1 (MQ=255)

                                                                        gTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTg                         <  1:3974549/101‑1 (MQ=255)

                                                                        gTCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTg                         >  2:3975006/1‑101 (MQ=255)

                                                                         tCTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTga                        >  1:3975009/1‑101 (MQ=255)

                                                                          cTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgac                       >  1:3975013/1‑101 (MQ=255)

                                                                          cTAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgac                       >  2:3975011/1‑101 (MQ=255)

                                                                           tAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacg                      >  1:3975018/1‑101 (MQ=255)

                                                                           tAACAACAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacg                      >  2:3975017/1‑101 (MQ=255)

                                                                            aacaacAATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacga                     <  1:3974592/101‑1 (MQ=255)

                                                                             acaacaATTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgac                    <  2:3974642/101‑1 (MQ=255)

                                                                              caacaaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgacg                   >  1:3975022/1‑101 (MQ=255)

                                                                               aacaaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgacga                  <  2:3974698/101‑1 (MQ=255)

                                                                               aacaaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgacga                  >  2:3975024/1‑101 (MQ=255)

                                                                                acaaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgacgac                 <  2:3974590/101‑1 (MQ=255)

                                                                                acaaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTgacgacgac                 <  1:3974494/101‑1 (MQ=255)

                                                                                  aaTTCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACAt               >  1:3975028/1‑101 (MQ=255)

                                                                                     tCCTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACAtttt            <  2:3974735/101‑1 (MQ=255)

                                                                                      ccTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTa           >  1:3975031/1‑101 (MQ=255)

                                                                                       cTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTAc          >  1:3975036/1‑101 (MQ=255)

                                                                                       cTTTAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTAc          >  2:3975035/1‑101 (MQ=255)

                                                                                          tAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTa       <  2:3974781/101‑1 (MQ=255)

                                                                                          tAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTa       >  1:3975042/1‑101 (MQ=255)

                                                                                          tAGTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTa       >  2:3975041/1‑101 (MQ=255)

                                                                                            gTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATc     <  1:3974780/101‑1 (MQ=255)

                                                                                            gTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATc     <  1:3974713/101‑1 (MQ=255)

                                                                                            gTGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATc     <  1:3974563/101‑1 (MQ=255)

                                                                                             tGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCt        >  2:3975046/1‑97 (MQ=255)

                                                                                             tGCGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCa    <  1:3974577/101‑1 (MQ=255)

                                                                                               cGACTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCAtt  <  2:3974683/101‑1 (MQ=255)

AAGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAAAGGCGAAAAAGGCTTCGTCTAACAACAATTCCTTTAGTGCGATTTACAAGGAATGGTACGAGCACAAGAAGCAAGTATGGTCAGTAGGGTATGCAACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCATT  >  NC_000913/2466733‑2466928

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