BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,302,652	A→G	Y513H (TAT→CAT)	yjcS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,302,652	0	A	G	98.3%	553.7 / ‑2.8	176	Y513H (TAT→CAT)	yjcS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (68/105); total (71/105)
Fisher's exact test for biased strand distribution p-value = 6.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 184 total aligned reads displayed.

TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741

tggtggTGTGGCGCATCTTATGCACCCCTGCGAGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCAATGTGGCGGACTCGGCCTGATGGc                                                                                         <  2:2097642/93‑1 (MQ=255)

 ggtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGcc                                                                                        >  2:132249/1‑93 (MQ=255)

  gtgCCGTGGCTGAACTTATGCCCCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGccc                                                                                       <  2:4404156/93‑1 (MQ=255)

  gtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGccc                                                                                       <  2:101436/93‑1 (MQ=255)

  gtgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGccc                                                                                       >  2:1497517/1‑93 (MQ=255)

   tgCCGTGGCTGAACTTATGCCCCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCa                                                                                      <  1:4007637/93‑1 (MQ=255)

   tgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAGCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCa                                                                                      <  2:2966035/93‑1 (MQ=255)

   tgCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCa                                                                                      >  2:2036000/1‑93 (MQ=255)

     ccGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGc                                                                                    <  2:3203441/93‑1 (MQ=255)

         ggCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGtt                                                                                <  1:346118/93‑1 (MQ=255)

          gCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTc                                                                               >  1:3332185/1‑93 (MQ=255)

           cTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCCGTGAGAAACCGCGCCAGGGGGCGGACTCGGCCTGAGGGCACAGCTgcacc                                                                              >  1:1097732/1‑89 (MQ=255)

           cTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCa                                                                              >  2:4301107/1‑93 (MQ=255)

            tGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaa                                                                             >  1:4471008/1‑93 (MQ=255)

            tGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCaa                                                                             >  1:2978677/1‑93 (MQ=255)

              aaCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAg                                                                           <  2:1956401/93‑1 (MQ=255)

              aaCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAg                                                                           <  2:657043/93‑1 (MQ=255)

              aaCTCATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGCCCTGATGGCCCAGCTGTTCAAAg                                                                           <  2:2982142/93‑1 (MQ=255)

               aCTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGt                                                                          <  1:2166489/93‑1 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        >  2:2254937/1‑93 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  2:308800/93‑1 (MQ=255)

                 ttATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTa                                                                        <  2:3448418/93‑1 (MQ=255)

                   aTGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTAtt                                                                      >  2:77512/1‑93 (MQ=255)

                    tGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTc                                                                     >  2:2157603/1‑93 (MQ=255)

                      cACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGc                                                                   >  1:3445981/1‑93 (MQ=255)

                       cccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCt                                                                  <  1:2359406/92‑1 (MQ=255)

                       cccccTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCt                                                                  <  1:4653517/92‑1 (MQ=255)

                       aCCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCt                                                                  <  2:613353/93‑1 (MQ=255)

                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  2:2448009/93‑1 (MQ=255)

                           cTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCa                                                              <  2:3871185/93‑1 (MQ=255)

                            ttCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAg                                                             >  1:872928/1‑93 (MQ=255)

                            ttCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAg                                                             >  1:1111337/1‑93 (MQ=255)

                               gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAtt                                                          <  1:3076118/93‑1 (MQ=255)

                               gcgcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAtt                                                          <  1:1737920/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttct                                                        <  2:4559655/93‑1 (MQ=255)

                                 gcAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttat                                                        >  1:2673518/1‑91 (MQ=255)

                                  cAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttctt                                                       <  1:3915492/93‑1 (MQ=255)

                                   aGCTTTTTCGCGCGGGTCAGGTAGAAACGGCGCTATGTGGCGGCTTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGAttcttc                                                      <  2:2663557/93‑1 (MQ=255)

                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     >  2:3455993/1‑93 (MQ=255)

                                    gCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcg                                                     <  2:3072568/93‑1 (MQ=255)

                                     ctctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcgc                                                    <  2:2474257/93‑1 (MQ=255)

                                      tctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTcgcg                                                   <  2:3277233/93‑1 (MQ=255)

                                       ctTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGa                                                  >  2:242251/1‑93 (MQ=255)

                                        tttCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGAc                                                 <  1:3616418/93‑1 (MQ=255)

                                         ttCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGAcc                                                <  2:3844881/93‑1 (MQ=255)

                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  1:4311884/93‑1 (MQ=255)

                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  1:3758594/93‑1 (MQ=255)

                                           cgcgCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTg                                              <  1:635576/93‑1 (MQ=255)

                                              gcCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTc                                           >  1:802150/1‑93 (MQ=255)

                                              gcCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTc                                           <  1:1630276/93‑1 (MQ=255)

                                               ccGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCa                                          <  1:758888/93‑1 (MQ=255)

                                                  gTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAccc                                       <  2:414197/93‑1 (MQ=255)

                                                  gTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCAccc                                       <  1:2017485/93‑1 (MQ=255)

                                                  gTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTATTCGCGACCTGGTCAccc                                       <  1:1061914/93‑1 (MQ=255)

                                                   tCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCg                                      <  1:3095517/93‑1 (MQ=255)

                                                   tCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCg                                      <  1:4700350/93‑1 (MQ=255)

                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   >  2:3851086/1‑93 (MQ=255)

                                                      ggTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGt                                   >  2:178067/1‑93 (MQ=255)

                                                       gTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGtt                                  <  1:1014164/93‑1 (MQ=255)

                                                         aGAAACCGCGCCATGTGGCGGACTCTGCCTGATAGCCCAGCTGTTCAAAGGTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTgg                                >  1:2717151/1‑93 (MQ=255)

                                                          gAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggc                               <  1:1298454/93‑1 (MQ=255)

                                                          gAAACCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggc                               <  2:2282216/93‑1 (MQ=255)

                                                             aCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcggc                            <  1:2827457/93‑1 (MQ=255)

                                                             aCCGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTggcggc                            >  1:4288834/1‑93 (MQ=255)

                                                              ccGCGCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCa                           >  1:3700240/1‑93 (MQ=255)

                                                               cgcgCCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCaa                          <  2:2340171/93‑1 (MQ=255)

                                                                gcgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAAt                         <  2:1967344/93‑1 (MQ=255)

                                                                 cgcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATc                        >  2:2717747/1‑93 (MQ=255)

                                                                  gcCATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCa                       <  1:3950407/93‑1 (MQ=255)

                                                                   ccATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAc                      >  1:1030385/1‑93 (MQ=255)

                                                                   ccATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAc                      >  1:4178630/1‑93 (MQ=255)

                                                                    cATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGTTTGGCGGCAATCAcc                     <  2:1985527/93‑1 (MQ=255)

                                                                    cATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAcc                     <  2:2655770/93‑1 (MQ=255)

                                                                    cATGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAcc                     <  2:494773/93‑1 (MQ=255)

                                                                    cATGTGGCGGACTCGGACTGATCGGCTAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCAcc                     <  2:1161533/93‑1 (MQ=255)

                                                                     gtgtgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGTCAATCACCt                    <  2:4004822/92‑1 (MQ=255)

                                                                     aTGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCt                    >  1:1345496/1‑93 (MQ=255)

                                                                     aTGTGGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCt                    <  1:1947543/93‑1 (MQ=255)

                                                                      tgtgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTg                   <  1:3699863/93‑1 (MQ=255)

                                                                       gtgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGt                  <  2:47645/93‑1 (MQ=255)

                                                                        tgGCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGtt                 >  1:2756573/1‑93 (MQ=255)

                                                                          gCGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTc               >  1:2667312/1‑93 (MQ=255)

                                                                          gCGGACTCGCCCTGATGGCCCAGCTGTTCAAAGGTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTc               <  1:4575747/93‑1 (MQ=255)

                                                                           cGGACTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTca              <  2:1313819/93‑1 (MQ=255)

                                                                           cGGACTCGGCCTGATGCCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTca              <  1:717534/93‑1 (MQ=255)

                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           <  1:1789591/93‑1 (MQ=255)

                                                                              aCTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           <  1:566144/93‑1 (MQ=255)

                                                                              aCTCGGCCTGATGCCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagca           <  2:1594684/93‑1 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          <  2:3946791/93‑1 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          <  1:1413691/93‑1 (MQ=255)

                                                                               cTCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTcagcag          >  1:659887/1‑93 (MQ=255)

                                                                                tCGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTCGTGACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGt         >  1:2905654/1‑93 (MQ=255)

                                                                                 cGGCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGGAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGtt        <  2:2027915/93‑1 (MQ=255)

                                                                                  ggCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTc       <  1:2193111/93‑1 (MQ=255)

                                                                                  ggCCTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTc       >  1:2382064/1‑93 (MQ=255)

                                                                                  ggCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTc       >  2:4562922/1‑93 (MQ=255)

                                                                                     cTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGc    <  1:1239904/93‑1 (MQ=255)

                                                                                     cTGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGc    >  1:4144998/1‑93 (MQ=255)

                                                                                      tGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGcc   <  2:3933777/93‑1 (MQ=255)

                                                                                      tGATGGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGcc   <  2:3570970/93‑1 (MQ=255)

TGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCAAAGTTATTCGCTTGCAGATTCTTCGCGACCTGGTCACCCGGGTTGGCGGCAATCACCTGTTTCAGCAGTTCTGCCG  >  NC_000913/4302562‑4302741

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