BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	400,314	A→G	Y204H (TAC→CAC)	ddlA ←	D‑alanine‑‑D‑alanine ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	400,314	0	A	G	100.0%	1092.5 / NA	256	Y204H (TAC→CAC)	ddlA	D‑alanine‑‑D‑alanine ligase A
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (135/121); total (135/121)

Only 100 of 256 total aligned reads displayed.

GTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCA  >  NC_000913/400216‑400411

gTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCt                                                                                                 <  1:592636/101‑1 (MQ=255)

 ttGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTg                                                                                                >  1:592879/1‑101 (MQ=255)

 ttGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTg                                                                                                <  2:592534/101‑1 (MQ=255)

 ttGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTg                                                                                                >  2:592877/1‑101 (MQ=255)

  tGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGt                                                                                               <  2:592435/101‑1 (MQ=255)

  tGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGt                                                                                               >  2:592882/1‑101 (MQ=255)

   gCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGtt                                                                                              >  1:592883/1‑101 (MQ=255)

    cccAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGttc                                                                                             <  1:592648/101‑1 (MQ=255)

       aGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGttcttc                                                                                          <  2:592326/101‑1 (MQ=255)

       aGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGttcttc                                                                                          >  1:592889/1‑101 (MQ=255)

       aGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGttcttc                                                                                          >  1:592887/1‑101 (MQ=255)

       aGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGttcttc                                                                                          <  1:592543/101‑1 (MQ=255)

           cTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTg                                                                                      >  2:592893/1‑101 (MQ=255)

           cTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTg                                                                                      >  1:592894/1‑101 (MQ=255)

           cTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTg                                                                                      <  1:592462/101‑1 (MQ=255)

            tGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTgg                                                                                     >  2:592896/1‑101 (MQ=255)

            tGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTgg                                                                                     <  2:592401/101‑1 (MQ=255)

             gcgcATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACt                                                                                       >  2:592899/1‑98 (MQ=255)

               gcATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTa                                                                                  <  2:592729/101‑1 (MQ=255)

                 aTTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACt                                                                                >  2:592904/1‑101 (MQ=255)

                 aTTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACt                                                                                >  1:592903/1‑101 (MQ=255)

                   tCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACttt                                                                              >  2:592907/1‑101 (MQ=255)

                   tCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACttt                                                                              >  2:592910/1‑101 (MQ=255)

                     gATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGc                                                                            >  1:592911/1‑101 (MQ=255)

                      aTCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCt                                                                           <  1:592524/101‑1 (MQ=255)

                       tctcACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTg                                                                          >  2:592915/1‑101 (MQ=255)

                         tcACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGac                                                                        >  1:592920/1‑101 (MQ=255)

                          cACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGaca                                                                       >  2:592923/1‑101 (MQ=255)

                          cACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGaca                                                                       <  1:592370/101‑1 (MQ=255)

                          cACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGaca                                                                       >  1:592922/1‑101 (MQ=255)

                            cGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACAcc                                                                     >  2:592928/1‑101 (MQ=255)

                             gACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCa                                                                    <  1:592358/101‑1 (MQ=255)

                             gACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCa                                                                    <  2:592516/101‑1 (MQ=255)

                             gACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCa                                                                    >  1:592929/1‑101 (MQ=255)

                              aCCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCa                                                                    >  1:592932/1‑100 (MQ=255)

                                 tttGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTg                                                                                      >  2:592934/1‑79 (MQ=255)

                                   tGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGt                                                                                    <  2:592451/79‑1 (MQ=255)

                                    gATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACagaag                                                             <  1:592418/101‑1 (MQ=255)

                                    gATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACagaag                                                             >  1:592937/1‑101 (MQ=255)

                                       cccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGc                                                          >  2:592941/1‑101 (MQ=255)

                                       cccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGc                                                          <  2:592558/101‑1 (MQ=255)

                                        ccTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGcc                                                         >  1:592943/1‑101 (MQ=255)

                                         cTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc                                                        <  2:592590/101‑1 (MQ=255)

                                          ttGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCt                                                       >  2:592944/1‑101 (MQ=255)

                                           tGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTg                                                      <  2:592611/101‑1 (MQ=255)

                                           tGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTg                                                      <  1:592509/101‑1 (MQ=255)

                                            gCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGa                                                     <  1:592710/101‑1 (MQ=255)

                                            gCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGa                                                     <  2:592479/101‑1 (MQ=255)

                                             cTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAt                                                    <  1:592602/101‑1 (MQ=255)

                                             cTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAt                                                    >  2:592949/1‑101 (MQ=255)

                                              tCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGccc                                                        >  2:592952/1‑96 (MQ=255)

                                              tCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAtt                                                   >  1:592951/1‑101 (MQ=255)

                                              tCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAtt                                                   <  2:592666/101‑1 (MQ=255)

                                               cAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTa                                                  <  2:592484/101‑1 (MQ=255)

                                                aaCGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAg                                                 <  2:592545/101‑1 (MQ=255)

                                                   gATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCg                                              >  2:592957/1‑101 (MQ=255)

                                                   gATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCg                                              <  2:592526/101‑1 (MQ=255)

                                                    aTCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCgg                                             >  2:592959/1‑101 (MQ=255)

                                                     tCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGt                                            >  1:592960/1‑101 (MQ=255)

                                                      cACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGtt                                           <  2:592688/101‑1 (MQ=255)

                                                       aCTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCt                                                                           >  2:592965/1‑68 (MQ=255)

                                                       aCTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGttt                                          >  1:592963/1‑101 (MQ=255)

                                                        cTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGtttt                                         >  1:592970/1‑101 (MQ=255)

                                                        cTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGtttt                                         <  2:592722/101‑1 (MQ=255)

                                                        cTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGtttt                                         <  1:592591/101‑1 (MQ=255)

                                                         tttATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGcc                                               >  2:592971/1‑94 (MQ=255)

                                                            aTGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACaa                                     >  1:592973/1‑101 (MQ=255)

                                                               aTCGAACTCGAACGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAt                                                    <  2:592505/83‑1 (MQ=255)

                                                                tcgaactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGAt                                                    <  2:592571/82‑1 (MQ=255)

                                                                tcgaactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAg                                 <  1:592576/101‑1 (MQ=255)

                                                                  gaactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTg                               <  1:592557/101‑1 (MQ=255)

                                                                  gaactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTg                               >  1:592977/1‑101 (MQ=255)

                                                                   aactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTgg                              >  2:592980/1‑101 (MQ=255)

                                                                   aactcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTgg                              >  1:592979/1‑101 (MQ=255)

                                                                     ctcgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTgg                              >  2:592984/1‑99 (MQ=255)

                                                                       cgaacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAAc                          >  1:592987/1‑101 (MQ=255)

                                                                         aacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAAccc                        >  1:592989/1‑101 (MQ=255)

                                                                         aacGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAAccc                        <  2:592585/101‑1 (MQ=255)

                                                                           cGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCa                      >  2:592991/1‑101 (MQ=255)

                                                                           cGCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCa                      <  2:592603/101‑1 (MQ=255)

                                                                            gCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAg                     <  1:592681/101‑1 (MQ=255)

                                                                            gCCAGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAg                     >  1:592993/1‑101 (MQ=255)

                                                                               aGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGttt                  <  1:592406/101‑1 (MQ=255)

                                                                               aGATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGttt                  <  1:592628/101‑1 (MQ=255)

                                                                                gATCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTa                 <  2:592741/101‑1 (MQ=255)

                                                                                 aTCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAg                <  2:592643/101‑1 (MQ=255)

                                                                                 aTCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAg                >  1:592997/1‑101 (MQ=255)

                                                                                  tCGACGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGa               <  1:592800/101‑1 (MQ=255)

                                                                                     aCGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTc            <  2:592707/101‑1 (MQ=255)

                                                                                      cGGCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTcc           >  1:593000/1‑101 (MQ=255)

                                                                                        gCAATTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAg                >  2:593001/1‑94 (MQ=255)

                                                                                          aaTTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACtt       <  1:592753/101‑1 (MQ=255)

                                                                                          aaTTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACtt       >  2:593006/1‑101 (MQ=255)

                                                                                          aaTTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACtt       >  1:593005/1‑101 (MQ=255)

                                                                                           aTTGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTc      <  2:592821/101‑1 (MQ=255)

                                                                                            ttGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCg     >  1:593010/1‑101 (MQ=255)

                                                                                            ttGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCg     >  2:593012/1‑101 (MQ=255)

                                                                                             tGCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCgg    >  1:593014/1‑101 (MQ=255)

                                                                                              gCGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGc   <  1:592820/101‑1 (MQ=255)

                                                                                               cGTGCTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCa  <  2:592616/101‑1 (MQ=255)

GTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATCGAACTCGAACGCCAGATCGACGGCAATTGCGTACTGTTCTTCACTGGTTACTTTGCTGACACCAACAGAAGAGCCCTGATTAGCCGGTTTTACAAACAGTGGTAACCCCAGTTTAGACTCCACTTCGGCA  >  NC_000913/400216‑400411

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: `‑`