BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,115,180	T→C	Q416R (CAG→CGG)	wcaK ←	putative galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,115,180	0	T	C	100.0%	353.2 / NA	112	Q416R (CAG→CGG)	wcaK	putative galactokinase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (57/55); total (57/55)

Only 100 of 116 total aligned reads displayed.

GTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCA  >  NC_000913/2115090‑2115269

gTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGc                                                                                         <  2:1197237/93‑1 (MQ=255)

gTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGc                                                                                         >  2:357113/1‑93 (MQ=255)

 ttAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCa                                                                                        <  1:1401398/93‑1 (MQ=255)

 ttAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCa                                                                                        >  1:1835560/1‑93 (MQ=255)

  tAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCAc                                                                                       >  1:3152953/1‑93 (MQ=255)

    aGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCa                                                                                     >  1:1189379/1‑93 (MQ=255)

     ggACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCAt                                                                                    <  1:569917/93‑1 (MQ=255)

     ggACGAAGGTCTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTGCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCAt                                                                                    <  2:2055935/93‑1 (MQ=255)

       aCGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCt                                                                                  <  2:2631173/93‑1 (MQ=255)

         gAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGc                                                                                <  2:2700613/93‑1 (MQ=255)

          aaGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCa                                                                               <  1:2509705/93‑1 (MQ=255)

             gTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTc                                                                            >  1:3492457/1‑93 (MQ=255)

              tttGTGACGACAGGGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCGGATGCGCTCAAGCACAGACCGCACCATCTGCATTcc                                                                           <  1:527023/93‑1 (MQ=255)

               ttCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCt                                                                          <  1:1849214/93‑1 (MQ=255)

                tCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTg                                                                         >  1:2688132/1‑93 (MQ=255)

                 cTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGt                                                                        <  1:922311/93‑1 (MQ=255)

                  tGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTc                                                                       >  1:3360808/1‑93 (MQ=255)

                  tGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTc                                                                       <  1:621517/93‑1 (MQ=255)

                   gacgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCt                                                                      <  2:298431/93‑1 (MQ=255)

                     cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGa                                                                    <  2:1103398/93‑1 (MQ=255)

                     cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGa                                                                    <  2:2304380/93‑1 (MQ=255)

                          gCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTc                                                               >  1:4299789/1‑93 (MQ=255)

                          gCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTc                                                               >  2:108574/1‑93 (MQ=255)

                          gCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTc                                                               >  1:488579/1‑93 (MQ=255)

                           cGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCTCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGGCGCTCa                                                              <  1:1778995/93‑1 (MQ=255)

                             ggAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACg                                                            >  1:2907617/1‑93 (MQ=255)

                              gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                           <  1:446630/93‑1 (MQ=255)

                              gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                           <  2:4085776/93‑1 (MQ=255)

                               aaTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGAc                                                          <  1:2433304/93‑1 (MQ=255)

                               aaTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGAc                                                          <  2:2865857/93‑1 (MQ=255)

                                 tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTg                                                        >  1:261776/1‑93 (MQ=255)

                                 tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTg                                                        >  1:1765481/1‑93 (MQ=255)

                                  ttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACAGCACCATCTGCATTCCTGTCTGACGCTCACGACTgc                                                       >  1:1778798/1‑92 (MQ=255)

                                   tCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                      <  1:1297744/93‑1 (MQ=255)

                                    cAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACg                                                     <  2:2683664/93‑1 (MQ=255)

                                     aGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACgg                                                    >  1:4653133/1‑93 (MQ=255)

                                         aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTc                                                >  1:3971322/1‑93 (MQ=255)

                                         aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTc                                                <  1:3113563/93‑1 (MQ=255)

                                          aaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCa                                               >  1:2948535/1‑93 (MQ=255)

                                           aagaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCat                                              >  1:2281221/1‑92 (MQ=255)

                                           aagaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc                                              >  1:4107594/1‑93 (MQ=255)

                                             gaagCCGCCCTTCATTTCGCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCCTCACTa                                            <  2:1450764/93‑1 (MQ=255)

                                             gaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTa                                            >  1:367052/1‑93 (MQ=255)

                                             gaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTa                                            >  2:4227205/1‑93 (MQ=255)

                                              gagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           <  1:2143874/92‑1 (MQ=255)

                                              aagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           <  1:3569272/93‑1 (MQ=255)

                                              aagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           <  2:4348366/93‑1 (MQ=255)

                                              aagCCGACCTTCATTCCGCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           <  1:1655914/93‑1 (MQ=255)

                                               agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAg                                          <  2:2267281/93‑1 (MQ=255)

                                                gCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGt                                         >  1:4192508/1‑93 (MQ=255)

                                                 ccGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTc                                        <  2:2903782/93‑1 (MQ=255)

                                                  cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGCCTGACTGCTTCACTAAGTcg                                       >  1:2894217/1‑93 (MQ=255)

                                                  cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcg                                       >  1:4337380/1‑93 (MQ=255)

                                                  cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcg                                       <  1:2320187/93‑1 (MQ=255)

                                                   gACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgc                                      <  2:1602802/93‑1 (MQ=255)

                                                   gACCTTCATGTCTCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTTACGCTCACGACGGACGGCTTCTCTAAGTcgc                                      <  2:4579169/93‑1 (MQ=255)

                                                    aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg                                     >  1:2953711/1‑93 (MQ=255)

                                                    aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg                                     >  1:2384129/1‑93 (MQ=255)

                                                     ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCgcgc                                    >  2:3070143/1‑93 (MQ=255)

                                                      cTTCATTTCTCCCCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                   <  2:143767/93‑1 (MQ=255)

                                                      cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                   <  2:2167260/93‑1 (MQ=255)

                                                        tCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCAtt                                 <  1:3883368/93‑1 (MQ=255)

                                                         cATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                                >  1:270131/1‑93 (MQ=255)

                                                         cATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                                >  2:2408639/1‑93 (MQ=255)

                                                          aTTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCa                               <  2:2208830/93‑1 (MQ=255)

                                                          aTTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCa                               >  1:1461609/1‑93 (MQ=255)

                                                           tttCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAg                              <  2:1835560/93‑1 (MQ=255)

                                                            ttCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgc                             <  2:1984061/93‑1 (MQ=255)

                                                              ccccTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcgc                           <  2:2951383/91‑1 (MQ=255)

                                                              cACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcgc                           <  2:3822193/93‑1 (MQ=255)

                                                               aCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAAGGCTTCACTAAGCCGCGCATTCAGCGcc                          >  1:2376129/1‑93 (MQ=255)

                                                                ccTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCg                         >  1:1055313/1‑93 (MQ=255)

                                                                  tCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGa                       >  1:351202/1‑93 (MQ=255)

                                                                  tCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGa                       >  2:3058720/1‑93 (MQ=255)

                                                                  tCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGa                       >  1:2386002/1‑93 (MQ=255)

                                                                  tCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGa                       >  1:4046809/1‑93 (MQ=255)

                                                                   ccccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGaa                      <  2:3684130/93‑1 (MQ=255)

                                                                   ccccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGaa                      <  1:3834399/93‑1 (MQ=255)

                                                                   ccccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGaa                      <  2:1394811/93‑1 (MQ=255)

                                                                    cccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAg                     >  2:1202258/1‑93 (MQ=255)

                                                                     ccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGc                    <  2:4016763/93‑1 (MQ=255)

                                                                     ccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGc                    <  2:1439161/93‑1 (MQ=255)

                                                                     ccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGc                    >  2:2790790/1‑93 (MQ=255)

                                                                       gATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTg                  <  1:199520/93‑1 (MQ=255)

                                                                       gATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTg                  >  1:485680/1‑93 (MQ=255)

                                                                       gATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTg                  >  2:1815461/1‑93 (MQ=255)

                                                                        aTGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTgg                 >  1:3975655/1‑93 (MQ=255)

                                                                        aTGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTgg                 >  1:623339/1‑93 (MQ=255)

                                                                        aTGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTgc                 >  1:1387947/1‑92 (MQ=255)

                                                                              tCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAg           <  2:3568019/93‑1 (MQ=255)

                                                                              tCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAg           >  1:3905335/1‑93 (MQ=255)

                                                                               cAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAgg          <  2:2063459/93‑1 (MQ=255)

                                                                                aggCACAGCCCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGA‑GGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        <  1:3613820/91‑1 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         <  2:1104541/93‑1 (MQ=255)

                                                                                 aGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        >  1:3012529/1‑93 (MQ=255)

                                                                                 aGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        <  2:2408473/93‑1 (MQ=255)

                                                                                 aGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        <  2:398414/93‑1 (MQ=255)

                                                                                    acaGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTCAGGTATcc     >  2:2708452/1‑93 (MQ=255)

                                                                                    acaGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATcc     >  2:3361531/1‑93 (MQ=255)

                                                                                       gACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCa  >  1:305638/1‑93 (MQ=255)

GTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCA  >  NC_000913/2115090‑2115269

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