BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,607,710	A→G	L332P (CTT→CCT)	uxaB ←	altronate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,607,710	0	A	G	100.0%	364.2 / NA	115	L332P (CTT→CCT)	uxaB	altronate oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (52/63); total (52/63)

Only 100 of 115 total aligned reads displayed.

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

   aGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                                                                                      <  1:473594/93‑1 (MQ=255)

   aGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                                                                                      <  1:2181279/93‑1 (MQ=255)

    gcgcGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAAt                                                                                     <  1:843375/93‑1 (MQ=255)

      gcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc                                                                                   >  2:2288230/1‑93 (MQ=255)

      gcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc                                                                                   <  1:3453435/93‑1 (MQ=255)

      gcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc                                                                                   <  2:1706052/93‑1 (MQ=255)

       cgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCa                                                                                  <  1:2565426/93‑1 (MQ=255)

        gATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg                                                                                 >  2:4476090/1‑93 (MQ=255)

        gATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg                                                                                 >  2:166078/1‑93 (MQ=255)

         aTAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt                                                                                >  2:4041294/1‑93 (MQ=255)

         aTAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt                                                                                <  1:435938/93‑1 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               >  1:1006793/1‑93 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               >  1:603512/1‑93 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         >  2:3611584/1‑93 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         <  1:709461/93‑1 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         <  2:1769901/93‑1 (MQ=255)

                 cagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaat                                                                        <  1:1978845/93‑1 (MQ=255)

                 cagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaat                                                                        <  1:1772699/93‑1 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  2:3435284/1‑93 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      <  2:3051602/93‑1 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      <  1:451755/93‑1 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  2:4036880/1‑93 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     >  2:3587722/1‑93 (MQ=255)

                     tGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCACCACGAGGCAAATCCAGTACCGGAATAAtt                                                                    >  1:2200207/1‑93 (MQ=255)

                      gATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt                                                                   >  1:2004222/1‑93 (MQ=255)

                      gATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt                                                                   >  2:4183527/1‑93 (MQ=255)

                        tGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttct                                                                 <  1:1201455/93‑1 (MQ=255)

                          cTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttcttc                                                               <  2:3766076/93‑1 (MQ=255)

                          cTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttcttc                                                               <  2:3548265/93‑1 (MQ=255)

                          cTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttcttc                                                               <  1:351544/93‑1 (MQ=255)

                             aaTGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTa                                                            <  1:2933936/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          >  1:4241117/1‑93 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  2:645461/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          >  1:4536192/1‑93 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  2:634959/93‑1 (MQ=255)

                                gTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGAt                                                         <  1:1971482/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  2:3796776/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        >  2:4265457/1‑93 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  1:2393978/93‑1 (MQ=255)

                                   cGGGTTACGGAAGCGCACGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGc                                                      <  1:570884/93‑1 (MQ=255)

                                   cGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGc                                                      <  1:2809656/93‑1 (MQ=255)

                                    gggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCt                                                     <  1:2700045/93‑1 (MQ=255)

                                     ggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtt                                                    >  2:611099/1‑93 (MQ=255)

                                     ggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtt                                                    >  1:71290/1‑93 (MQ=255)

                                      gTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttt                                                   >  1:4089161/1‑93 (MQ=255)

                                      gTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttt                                                   >  1:4471957/1‑93 (MQ=255)

                                      gTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttt                                                   >  1:3530404/1‑93 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  <  2:3153845/93‑1 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  <  2:2639454/93‑1 (MQ=255)

                                        tACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttttt                                                 <  2:4058981/93‑1 (MQ=255)

                                           ggAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCaa                                              <  1:1647544/93‑1 (MQ=255)

                                            gAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAAc                                             >  1:4595289/1‑93 (MQ=255)

                                              aaCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGa                                           <  2:2136150/93‑1 (MQ=255)

                                               aCGCCCGGTGACTGCACTGGCGAAACATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          <  1:4528722/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  2:142661/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  2:3912326/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  2:3926595/93‑1 (MQ=255)

                                                 ccccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATg                                        <  2:4477669/92‑1 (MQ=255)

                                                   ccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGCATgcg                                      >  2:786574/1‑93 (MQ=255)

                                                   ccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcg                                      >  2:674795/1‑93 (MQ=255)

                                                   ccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcg                                      >  1:2803302/1‑93 (MQ=255)

                                                   ccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcg                                      >  1:2990882/1‑93 (MQ=255)

                                                     ggTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCa                                    >  2:2273301/1‑93 (MQ=255)

                                                     ggTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCa                                    >  2:1229618/1‑93 (MQ=255)

                                                      gTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCaa                                   >  1:3500066/1‑93 (MQ=255)

                                                       tGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCaaa                                  >  2:4373082/1‑93 (MQ=255)

                                                       tGACTGCACTGGCGAAAGATTCAAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCaaa                                  <  1:909622/93‑1 (MQ=255)

                                                           tGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTc                              <  1:906583/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  2:2600135/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  1:3415323/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  1:4334411/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  1:677870/93‑1 (MQ=255)

                                                             cACTGGCTAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAg                            <  1:972094/93‑1 (MQ=255)

                                                              aCTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGc                           >  2:1204087/1‑93 (MQ=255)

                                                               cTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCa                          <  2:2645065/93‑1 (MQ=255)

                                                               cTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCa                          <  1:37042/93‑1 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       >  1:4345362/1‑93 (MQ=255)

                                                                   cGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGt                      <  2:649646/93‑1 (MQ=255)

                                                                   cGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGt                      <  2:3935829/93‑1 (MQ=255)

                                                                   cGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGt                      <  2:448028/93‑1 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTc                    <  2:4019488/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   <  1:2861215/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   <  2:2823672/93‑1 (MQ=255)

                                                                       aGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCAt                  >  2:2834208/1‑93 (MQ=255)

                                                                       aGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCAt                  >  1:3573509/1‑93 (MQ=255)

                                                                        gATTCCAGTTCATCACGAGGCAAATCCAGTACCGGGATAATTTCTTCGTAGATAGCTTTTTCAACGAATGAGCAAATTTAAGCATCGTTCATc                 >  2:373225/1‑93 (MQ=255)

                                                                        gATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 >  2:319591/1‑93 (MQ=255)

                                                                         aTTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCg                >  1:2112321/1‑93 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               <  1:3593484/93‑1 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               <  1:4015820/93‑1 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               <  2:1719132/93‑1 (MQ=255)

                                                                           tCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCt              <  2:3085448/93‑1 (MQ=255)

                                                                              aGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCa           >  2:3978899/1‑93 (MQ=255)

                                                                              aGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCa           >  2:2973584/1‑93 (MQ=255)

                                                                               gTTCATCAGGAGGCAAATCCTGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAc          >  2:4583594/1‑93 (MQ=255)

                                                                               gTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAc          >  1:4141381/1‑93 (MQ=255)

                                                                                ttCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAcc         <  1:2627341/93‑1 (MQ=255)

                                                                                 tcatcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCt        >  2:2781866/1‑93 (MQ=255)

                                                                                  catcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTa       <  2:2072434/93‑1 (MQ=255)

                                                                                       cGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTa  >  2:2763989/1‑93 (MQ=255)

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

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