BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,198,118	A→G	D87G (GAC→GGC)	hupA →	DNA‑binding protein HU‑alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,198,118	0	A	G	100.0%	1254.0 / NA	399	D87G (GAC→GGC)	hupA	DNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (173/226); total (173/226)

Only 100 of 404 total aligned reads displayed.

AGCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198028‑4198207

 gCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGgcgc                                                                                        >  1:5703137/1‑93 (MQ=255)

    tACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGt                                                                                     >  1:7077837/1‑93 (MQ=255)

    tACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGt                                                                                     <  2:3153019/93‑1 (MQ=255)

      cTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTta                                                                                   <  1:2783362/93‑1 (MQ=255)

       tGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaa                                                                                  <  1:5891076/93‑1 (MQ=255)

        ggCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaag                                                                                 >  2:1175656/1‑93 (MQ=255)

        ggCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaag                                                                                 <  2:1540619/93‑1 (MQ=255)

         gCCGCATCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAGCGTACCTTCATTTGTTTCGGGCAAGGCACTGAAAGGCGCAGTtaagt                                                                                <  2:299326/93‑1 (MQ=38)

         gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagt                                                                                <  1:7564279/93‑1 (MQ=255)

          ccGCAACCCGCAGACCGGTAAAGCAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagta                                                                               <  2:4357788/93‑1 (MQ=255)

          ccGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagta                                                                               <  1:4070577/93‑1 (MQ=255)

           cGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaa                                                                              >  2:4022214/1‑93 (MQ=255)

            gCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaag                                                                             >  2:5196873/1‑93 (MQ=255)

             cAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGa                                                                            <  1:666280/93‑1 (MQ=255)

                  cGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGTCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCg                                                                       <  2:7797793/93‑1 (MQ=255)

                  cGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCg                                                                       <  2:7153508/93‑1 (MQ=255)

                   gCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGt                                                                      >  1:2504550/1‑93 (MQ=255)

                     aGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTgg                                                                    >  2:6852767/1‑93 (MQ=255)

                     aGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTgg                                                                    >  1:7266481/1‑93 (MQ=255)

                      gACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGc                                                                   <  2:6306895/93‑1 (MQ=255)

                      gACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGc                                                                   <  2:6992425/93‑1 (MQ=255)

                      gACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGc                                                                   >  2:497135/1‑93 (MQ=255)

                      gACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGc                                                                   >  1:5351307/1‑93 (MQ=255)

                       aCCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCa                                                                  <  2:564477/93‑1 (MQ=255)

                        ccGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAg                                                                 >  1:5812891/1‑93 (MQ=255)

                          ggTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTg                                                               >  2:5250912/1‑93 (MQ=255)

                          ggTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTg                                                               >  2:4939046/1‑93 (MQ=255)

                           gTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGa                                                              >  1:2071585/1‑93 (MQ=255)

                            tAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGaa                                                             <  2:1245481/93‑1 (MQ=255)

                            tAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGaa                                                             >  1:3077259/1‑93 (MQ=255)

                             aaaGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAAc                                                            <  1:5174377/93‑1 (MQ=255)

                              aaGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACa                                                           <  1:1529812/93‑1 (MQ=255)

                               aGAAATCAAAATTGCCGCAGCTAACGTACCGGCATGTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAg                                                          <  1:1106607/93‑1 (MQ=255)

                                gAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGt                                                         <  2:2397041/93‑1 (MQ=255)

                                 aaaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtt                                                        <  2:1913327/93‑1 (MQ=255)

                                 aaaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtt                                                        >  2:4578218/1‑93 (MQ=255)

                                 aaaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtt                                                        <  2:5715527/93‑1 (MQ=255)

                                  aaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGttt                                                       <  2:2450917/93‑1 (MQ=255)

                                    tCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTa                                                     >  1:2913865/1‑93 (MQ=255)

                                     cAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTaa                                                    >  1:2999487/1‑93 (MQ=255)

                                     cAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTaa                                                    <  2:4847327/93‑1 (MQ=255)

                                     cAAAATTGCCGCAGCT‑ACGTACCGGCATTTGTTTCTGG‑AAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACg                                                  >  1:1249138/1‑93 (MQ=255)

                                      aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   >  1:3739926/1‑93 (MQ=255)

                                      aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   >  1:1748190/1‑93 (MQ=255)

                                      aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   >  2:7339687/1‑93 (MQ=255)

                                      aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   >  1:5502782/1‑93 (MQ=255)

                                       aaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACg                                                  >  2:306200/1‑93 (MQ=255)

                                        aaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGa                                                 <  2:5094976/93‑1 (MQ=255)

                                         aTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGaa                                                <  1:1896774/93‑1 (MQ=255)

                                         aTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGaa                                                <  1:3594393/93‑1 (MQ=255)

                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               <  1:6980676/93‑1 (MQ=255)

                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              >  2:2147350/1‑93 (MQ=255)

                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:802817/93‑1 (MQ=255)

                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  1:5570207/93‑1 (MQ=255)

                                             ccGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgggg                                            <  1:6028455/93‑1 (MQ=255)

                                               gCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtg                                          <  2:5541590/93‑1 (MQ=255)

                                                 aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtggt                                        <  2:137233/93‑1 (MQ=255)

                                                      aCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAc                                   <  2:2627320/93‑1 (MQ=255)

                                                       cGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcc                                  <  2:3325497/93‑1 (MQ=255)

                                                         tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                <  2:5382911/93‑1 (MQ=255)

                                                         tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                <  2:3406491/93‑1 (MQ=255)

                                                          aCCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCt                               <  2:5935029/93‑1 (MQ=255)

                                                           ccGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtt                              <  2:4482476/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAGCAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:7252642/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:1020976/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  1:4438496/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:6645551/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:7639355/93‑1 (MQ=255)

                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  1:7678303/93‑1 (MQ=255)

                                                              gCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTg                           <  2:367853/93‑1 (MQ=255)

                                                              gCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTg                           <  1:4147434/93‑1 (MQ=255)

                                                              gCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGGTTCAACCCTTTTg                           >  1:4359670/1‑93 (MQ=255)

                                                                 tttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCt                        >  1:2295765/1‑93 (MQ=255)

                                                                  ttGTTTCTGGCAGGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCtt                       >  1:7415186/1‑93 (MQ=255)

                                                                  ttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCtt                       >  1:2654700/1‑93 (MQ=255)

                                                                  ttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCtt                       <  1:451462/93‑1 (MQ=255)

                                                                    gTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTc                     >  2:2696194/1‑93 (MQ=255)

                                                                    gTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTc                     <  2:1808572/93‑1 (MQ=255)

                                                                     tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    >  1:5004045/1‑93 (MQ=255)

                                                                      ttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTg                   >  2:6440032/1‑93 (MQ=255)

                                                                        cTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGCGGTTTCACTCCTTTTGTGTTTCTGGc                 >  1:1493540/1‑93 (MQ=255)

                                                                         tGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCg                >  1:4694246/1‑93 (MQ=255)

                                                                              aGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAACGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGAt           >  1:2112990/1‑93 (MQ=255)

                                                                               ggCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATc          <  1:1670339/93‑1 (MQ=255)

                                                                               ggCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATc          <  1:6220882/93‑1 (MQ=255)

                                                                                gCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa         >  2:4620324/1‑93 (MQ=255)

                                                                                gCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa         <  1:1596377/93‑1 (MQ=255)

                                                                                gCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa         <  1:3898721/93‑1 (MQ=255)

                                                                                gCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa         <  1:6386243/93‑1 (MQ=255)

                                                                                  aCTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAtt       <  1:4908837/93‑1 (MQ=255)

                                                                                  aCTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAtt       >  1:6705797/1‑93 (MQ=255)

                                                                                  aCTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAtt       <  2:65051/93‑1 (MQ=255)

                                                                                   cTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg      >  1:3292525/1‑93 (MQ=255)

                                                                                   cTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg      <  2:4310244/93‑1 (MQ=255)

                                                                                    tGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa     >  1:7490117/1‑93 (MQ=255)

                                                                                    tGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa     >  1:888801/1‑93 (MQ=255)

                                                                                      aaaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATg   <  2:3892991/93‑1 (MQ=255)

                                                                                       aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  2:3810737/93‑1 (MQ=255)

                                                                                       aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  >  1:3506910/1‑93 (MQ=255)

                                                                                       aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  2:4413405/93‑1 (MQ=255)

AGCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198028‑4198207

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