BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,560	G→A	A144T (GCC→ACC)	yjfJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,560	0	G	A	99.7%	1253.8 / ‑5.9	361	A144T (GCC→ACC)	yjfJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (183/176); total (183/178)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

Only 100 of 368 total aligned reads displayed.

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACGGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        <  1:970554/93‑1 (MQ=255)

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        >  1:5173587/1‑93 (MQ=255)

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        <  2:6019047/93‑1 (MQ=255)

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        <  1:3857695/93‑1 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       >  2:7501176/1‑93 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  1:4851032/93‑1 (MQ=255)

  ggAAGGCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  2:6850669/93‑1 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      <  2:3956480/93‑1 (MQ=255)

       tctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAg                                                                                  >  1:3895946/1‑93 (MQ=255)

        ctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGc                                                                                 <  1:7386950/93‑1 (MQ=255)

        ctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGc                                                                                 <  2:1391894/93‑1 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                <  1:4487754/93‑1 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              <  1:6309655/93‑1 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              <  2:4968677/93‑1 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              >  1:6429597/1‑93 (MQ=255)

             gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGc                                                                            >  2:6317886/1‑93 (MQ=255)

             gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGc                                                                            >  1:5760182/1‑93 (MQ=255)

               gATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAg                                                                          >  2:5668176/1‑93 (MQ=255)

                aTGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGc                                                                         <  2:7426604/93‑1 (MQ=255)

                aTGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGc                                                                         <  2:1560/93‑1 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       <  1:7065249/93‑1 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       <  2:1158887/93‑1 (MQ=255)

                    cGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGc                                                                     >  2:351474/1‑93 (MQ=255)

                    cGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGc                                                                     >  2:1021724/1‑93 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    <  2:7621113/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    <  2:5201080/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  1:2107996/1‑93 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  1:5996244/1‑93 (MQ=255)

                      tGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcac                                                                   >  2:2154525/1‑93 (MQ=255)

                        gAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaa                                                                 <  1:637564/93‑1 (MQ=255)

                         aaaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCCCTGAAGCCATGCAGCGTGCacaac                                                                <  2:3175051/93‑1 (MQ=255)

                          aaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               >  2:7699158/1‑93 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              >  2:4837872/1‑93 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              >  2:6448044/1‑93 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              >  2:3332116/1‑93 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              <  1:3585519/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           <  2:3646806/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           >  1:6859759/1‑93 (MQ=255)

                               cAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTg                                                          >  1:4165850/1‑93 (MQ=255)

                               cAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGAACAACAGGCTg                                                          >  1:2544977/1‑93 (MQ=255)

                                 gTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTa                                                        <  1:6600150/93‑1 (MQ=255)

                                  ttACAGCGACAGCGCAGCGGATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaa                                                       >  2:810973/1‑93 (MQ=255)

                                   tACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaac                                                      <  2:2758/93‑1 (MQ=255)

                                     cagcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaa                                                    <  1:7022154/93‑1 (MQ=255)

                                      agcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaac                                                   >  1:1020244/1‑93 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  1:2021276/93‑1 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               <  1:1352646/93‑1 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               >  1:2934071/1‑93 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  2:2585202/93‑1 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  1:3623710/93‑1 (MQ=255)

                                            agcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAc                                             >  2:2046126/1‑93 (MQ=255)

                                            agcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAc                                             <  2:5481188/93‑1 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            >  2:109860/1‑93 (MQ=255)

                                               gcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGt                                          <  2:7227850/93‑1 (MQ=255)

                                                cagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACCGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGtt                                         <  2:5359916/93‑1 (MQ=255)

                                                 agcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTg                                        <  1:511647/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  2:3374013/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  1:4982431/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  2:6197128/93‑1 (MQ=255)

                                                   cgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgc                                      <  2:7929628/93‑1 (MQ=255)

                                                    gTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcg                                     >  2:6787740/1‑93 (MQ=255)

                                                    gTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcg                                     >  2:5947991/1‑93 (MQ=255)

                                                      aTCTCTCAGTTTGAGCAGCAAATGGAAGTCGCTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCa                                   <  1:7819086/93‑1 (MQ=255)

                                                         gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                <  1:5914434/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  2:3834804/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  1:3023294/93‑1 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            <  2:1301062/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           >  1:7512154/1‑93 (MQ=255)

                                                               tttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                          <  1:4372382/93‑1 (MQ=255)

                                                                ttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCCTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCg                         <  1:5466126/93‑1 (MQ=255)

                                                                  gAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGtt                       <  2:5884722/93‑1 (MQ=255)

                                                                   agcagcAAATGAAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGttt                      <  2:2987710/93‑1 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAATCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     <  2:4276748/93‑1 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     <  1:1043331/93‑1 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     >  1:6211084/1‑93 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     >  2:4604022/1‑93 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     <  1:5572050/93‑1 (MQ=255)

                                                                      agcaAATGGAAGTCGTTAAAACCACTGAAGCCATGGAGCGTGCACAACAGGATGTCACAACGTCTACCGTTGGCGCCTCTTCAAGCGTTTCGa                   >  1:4130191/1‑93 (MQ=255)

                                                                      agcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGa                   >  1:4742558/1‑93 (MQ=255)

                                                                      agcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGa                   >  2:3075166/1‑93 (MQ=255)

                                                                      agcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGa                   <  1:6105818/93‑1 (MQ=255)

                                                                       gcaAATGGAAGTCGTTAAAACTACTGAAGCCATGCAGCGTGCACAACGGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                  >  1:354521/1‑93 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                <  1:1226122/93‑1 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACAGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                >  2:6947925/1‑93 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               <  1:3414739/93‑1 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               <  1:3155543/93‑1 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               >  1:5058969/1‑93 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              <  1:7243095/93‑1 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              <  1:5252878/93‑1 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              >  2:6080683/1‑93 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              <  2:4453764/93‑1 (MQ=255)

                                                                              gAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCa           <  1:2737466/93‑1 (MQ=255)

                                                                               aaGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAg          >  1:2773649/1‑93 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTAGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        >  2:6264914/1‑93 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       <  1:1897742/93‑1 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       <  1:7721913/93‑1 (MQ=255)

                                                                                    gTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATcc     <  1:1493275/93‑1 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  2:2254645/1‑93 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  2:1755841/1‑93 (MQ=255)

                                                                                       aaaaCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTa  <  2:2532092/93‑1 (MQ=255)

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

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