BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,507,420	T→C	T18A (ACG→GCG)	b1437 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,507,420	0	T	C	99.5%	606.9 / ‑6.8	190	T18A (ACG→GCG)	b1437	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (78/110); total (79/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 197 total aligned reads displayed.

TGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATG  >  NC_000913/1507334‑1507508

tGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCAcccc                                                                                    <  2:4141309/93‑1 (MQ=255)

 ggCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCt                                                                                   <  2:2912923/93‑1 (MQ=255)

 ggCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCt                                                                                   <  1:1597358/93‑1 (MQ=255)

  gCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTg                                                                                  <  1:939606/93‑1 (MQ=255)

   cTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCACTGc                                                                                 <  1:2003284/93‑1 (MQ=255)

     cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                               <  1:829974/93‑1 (MQ=255)

      gAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGAt                                                                              <  2:2604952/93‑1 (MQ=255)

        aTCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATga                                                                            <  2:1733207/93‑1 (MQ=255)

        aTCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATga                                                                            <  1:3700990/93‑1 (MQ=255)

         tctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgag                                                                           >  2:839692/1‑93 (MQ=255)

          ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                          <  1:1724832/93‑1 (MQ=255)

            cAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAtt                                                                        <  1:4302336/93‑1 (MQ=255)

            cAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAtt                                                                        <  1:2885817/93‑1 (MQ=255)

             aGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTa                                                                       <  1:3596776/93‑1 (MQ=255)

              gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaa                                                                      <  2:1409505/93‑1 (MQ=255)

              gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaa                                                                      <  1:778855/93‑1 (MQ=255)

               ggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaaa                                                                     <  1:226367/93‑1 (MQ=255)

                gCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCATGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAg                                                                    <  1:1450601/93‑1 (MQ=255)

                gCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAg                                                                    <  2:2129234/93‑1 (MQ=255)

                gCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAg                                                                    <  1:1248989/93‑1 (MQ=255)

                  gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                  >  1:239897/1‑93 (MQ=255)

                  gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                  >  1:235075/1‑93 (MQ=255)

                   tCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAAc                                                                 >  2:2322371/1‑93 (MQ=255)

                    cGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAAcc                                                                >  1:243391/1‑93 (MQ=255)

                     gATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCa                                                               <  2:4637165/93‑1 (MQ=255)

                      aTGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAt                                                              >  1:2899362/1‑93 (MQ=255)

                      aTGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAt                                                              <  1:2091436/93‑1 (MQ=255)

                      aTGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAt                                                              <  2:403385/93‑1 (MQ=255)

                       tctAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAtt                                                             <  1:3628657/91‑1 (MQ=255)

                       tGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAtt                                                             >  1:307236/1‑93 (MQ=255)

                        gTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTg                                                            >  2:377600/1‑93 (MQ=255)

                         tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                           <  2:1655627/93‑1 (MQ=255)

                         tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                           >  2:634784/1‑93 (MQ=255)

                           gCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGt                                                         >  1:2938787/1‑93 (MQ=255)

                            cGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTa                                                        <  2:1695308/93‑1 (MQ=255)

                             gAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTaa                                                       >  1:1422452/1‑93 (MQ=255)

                             gAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTaa                                                       <  1:392183/93‑1 (MQ=255)

                               aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                     >  2:4186350/1‑93 (MQ=255)

                                  gcagcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCTCGACACCCCTGCGATGAGATTAAAGAACCATTCCGTACAGCaa                                                  >  2:81316/1‑93 (MQ=255)

                                  gcagcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCaa                                                  <  1:1852969/93‑1 (MQ=255)

                                   cagcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAAt                                                 <  2:4492495/93‑1 (MQ=255)

                                    agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                <  2:1504198/93‑1 (MQ=255)

                                     gcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATcc                                               <  1:3731292/93‑1 (MQ=255)

                                     gcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATcc                                               >  2:1508139/1‑93 (MQ=255)

                                     gcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATcc                                               >  2:2078141/1‑93 (MQ=255)

                                       aaCCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTg                                             >  1:4538316/1‑93 (MQ=255)

                                       aaCCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTg                                             >  1:4661191/1‑93 (MQ=255)

                                          cATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaaa                                          >  2:3588016/1‑93 (MQ=255)

                                          cATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaaa                                          >  2:2021516/1‑93 (MQ=255)

                                           aTTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaac                                         >  1:3530267/1‑93 (MQ=255)

                                           aTTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaac                                         >  1:3821734/1‑93 (MQ=255)

                                            tttGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaaca                                        >  2:264491/1‑93 (MQ=255)

                                               gAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACt                                     <  1:296827/93‑1 (MQ=255)

                                                aagcTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTc                                    <  1:1828651/90‑1 (MQ=255)

                                                 aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                   <  2:1093292/93‑1 (MQ=255)

                                                 aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                   >  2:523195/1‑93 (MQ=255)

                                                 aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                   >  1:642022/1‑93 (MQ=255)

                                                 aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                   <  1:4572019/93‑1 (MQ=255)

                                                  aCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCgg                                  <  1:2157973/93‑1 (MQ=255)

                                                  aCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCgg                                  <  1:4745779/93‑1 (MQ=255)

                                                  aCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCgg                                  >  2:1264256/1‑93 (MQ=255)

                                                  aCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCgg                                  >  2:1983584/1‑93 (MQ=255)

                                                   cTCAGGTTTCATGGGAGGCGCAGTGTCAGGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                 <  1:74777/93‑1 (MQ=255)

                                                   cTCAGGGTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                 <  1:4248151/93‑1 (MQ=255)

                                                      agctttCATGGGAGGCGCAGTGCCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTg                              <  2:1795426/90‑1 (MQ=255)

                                                      aGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTg                              >  1:497978/1‑93 (MQ=255)

                                                        gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                            <  1:303764/93‑1 (MQ=255)

                                                         tttCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGt                           <  1:405328/93‑1 (MQ=255)

                                                            cATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTaa                        <  1:2687022/93‑1 (MQ=255)

                                                             aTGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAAt                       >  1:2056769/1‑93 (MQ=255)

                                                             aTGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAAt                       <  1:229290/93‑1 (MQ=255)

                                                              tGGTAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATc                      <  1:166294/93‑1 (MQ=255)

                                                              tGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATc                      >  1:1048121/1‑93 (MQ=255)

                                                                  aGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCa                  <  1:2565769/93‑1 (MQ=255)

                                                                  aGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCa                  >  2:4592899/1‑93 (MQ=255)

                                                                   ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                 <  1:4128610/93‑1 (MQ=255)

                                                                   ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                 >  2:3398467/1‑93 (MQ=255)

                                                                   ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                 <  2:4334809/93‑1 (MQ=255)

                                                                   ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                 <  1:3265955/93‑1 (MQ=255)

                                                                    gcgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGTAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAt                <  2:4294189/93‑1 (MQ=255)

                                                                      gcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTa              <  2:652609/93‑1 (MQ=255)

                                                                       cAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaa             >  2:4242384/1‑93 (MQ=255)

                                                                       cAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaa             >  2:718090/1‑93 (MQ=255)

                                                                        aGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaa            >  1:427046/1‑93 (MQ=255)

                                                                         gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaa           <  2:3732425/93‑1 (MQ=255)

                                                                         gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaa           <  1:2783432/93‑1 (MQ=255)

                                                                          tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa          <  2:4013446/93‑1 (MQ=255)

                                                                          tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa          <  2:3138758/93‑1 (MQ=255)

                                                                          tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa          <  1:4433775/93‑1 (MQ=255)

                                                                           gtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAg         >  1:3128062/1‑93 (MQ=255)

                                                                            tCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAgg        <  1:772224/93‑1 (MQ=255)

                                                                             cATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGt       <  1:1007497/93‑1 (MQ=255)

                                                                             cATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGt       <  2:1319720/93‑1 (MQ=255)

                                                                              aTGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGtt      >  1:3945960/1‑93 (MQ=255)

                                                                              aTGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGtt      >  2:2034572/1‑93 (MQ=255)

                                                                                gCCGCGCCACCCCTGCGCTGAGATTAAAGATCCATTGCGTAAAGCAATCCTGAAACAACTCGGTGTGAGTTAATCGCCAATTAAAAAGGTTaa    <  2:1354971/93‑1 (MQ=255)

                                                                                gCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTGTGAGGGAAGCGCCAATTAAAAAGGTTaa    >  1:4210428/1‑93 (MQ=255)

                                                                                  cgcgCCACCGCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATg  <  1:3387361/93‑1 (MQ=255)

                                                                                  cgcgCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATg  <  2:4483691/93‑1 (MQ=255)

                                                                                  cgcgCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATg  <  2:2757371/93‑1 (MQ=255)

TGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATG  >  NC_000913/1507334‑1507508

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