BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,561,521	C→T	S754N (AGC→AAC)	b1489 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,561,521	0	C	T	99.1%	373.7 / ‑5.2	109	S754N (AGC→AAC)	b1489	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (45/63); total (45/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 113 total aligned reads displayed.

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAT  >  NC_000913/1561431‑1561609

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                        >  2:1773131/1‑93 (MQ=255)

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                        <  2:2928061/93‑1 (MQ=255)

  tCCCTGAATAACGCGACAGTGGATCTTGCGTAGCCTCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTg                                                                                      >  1:2245413/1‑93 (MQ=255)

  tCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTg                                                                                      >  2:722021/1‑93 (MQ=255)

  tCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTg                                                                                      >  1:3830062/1‑93 (MQ=255)

    ccTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTgtc                                                                                    >  2:1058017/1‑91 (MQ=255)

    ccTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGcc                                                                                    >  1:3614194/1‑93 (MQ=255)

     cTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGccc                                                                                   <  2:3133459/93‑1 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                               >  1:2989524/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCGCGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                              >  1:2195877/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                              >  2:1441908/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                              <  2:4146598/93‑1 (MQ=255)

            aCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTg                                                                            <  2:3035635/93‑1 (MQ=255)

            aCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTg                                                                            <  1:2051865/93‑1 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                        >  2:4060176/1‑93 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                       <  1:1687619/93‑1 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                       >  1:2525405/1‑93 (MQ=255)

                  cAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATg                                                                      >  2:1292408/1‑93 (MQ=255)

                     tGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCt                                                                   <  2:77925/93‑1 (MQ=255)

                     tGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCt                                                                   <  2:2571921/93‑1 (MQ=255)

                     tGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCt                                                                   >  2:1804322/1‑93 (MQ=255)

                      ggATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCtt                                                                  <  1:1208224/93‑1 (MQ=255)

                      ggATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCtt                                                                  <  2:373801/93‑1 (MQ=255)

                      ggATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCtt                                                                  <  1:3468509/93‑1 (MQ=255)

                       gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTc                                                                 >  2:342468/1‑93 (MQ=255)

                        aTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCa                                                                >  2:4375477/1‑93 (MQ=255)

                        aTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCa                                                                >  1:347919/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                               >  2:699723/1‑93 (MQ=255)

                           ttGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGt                                                             <  1:4038024/93‑1 (MQ=255)

                            tGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTa                                                            <  1:267810/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                          <  2:1503441/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                          <  2:3270360/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                          <  2:1093455/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTGTGGTTTCGACGCCTTCCGCCACGACGTTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                          <  2:1858178/93‑1 (MQ=255)

                               gTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAgg                                                         >  2:627901/1‑93 (MQ=255)

                               gTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAgg                                                         >  1:962297/1‑93 (MQ=255)

                               gTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAgg                                                         >  2:3141150/1‑93 (MQ=255)

                                tAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGc                                                        >  1:924025/1‑93 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                      >  1:1329852/1‑93 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                      <  1:106232/93‑1 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                      >  2:583036/1‑93 (MQ=255)

                                   cATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAg                                                     <  1:3171141/93‑1 (MQ=255)

                                    aTCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAgg                                                    >  1:2123435/1‑93 (MQ=255)

                                     tctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGa                                                   <  2:2481079/93‑1 (MQ=255)

                                     tctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGa                                                   <  2:1864299/93‑1 (MQ=255)

                                     tctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGa                                                   >  2:4246735/1‑93 (MQ=255)

                                      ctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGAt                                                  <  1:2547203/93‑1 (MQ=255)

                                      ctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGAt                                                  <  1:195840/93‑1 (MQ=255)

                                      ctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGGAATGGCTTCAAGTAAGGCAAGGAt                                                  >  2:4597880/1‑93 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                 >  1:4206332/1‑93 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                 <  1:3241987/93‑1 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                 <  2:4222596/93‑1 (MQ=255)

                                          aaTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGt                                              <  2:3208861/93‑1 (MQ=255)

                                           aTTGCTCTTTTGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCggt                                             >  1:2862907/1‑91 (MQ=255)

                                           aTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtt                                             <  1:2742873/93‑1 (MQ=255)

                                           aTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGCCAAGGATGCGtt                                             <  2:3759355/93‑1 (MQ=255)

                                            ttGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttt                                            <  1:4221257/93‑1 (MQ=255)

                                             tGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtttt                                           <  2:433557/93‑1 (MQ=255)

                                              gCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttttt                                          >  1:1210865/1‑93 (MQ=255)

                                               ctctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTc                                         <  1:2816254/93‑1 (MQ=255)

                                                tctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCg                                        >  1:4632474/1‑93 (MQ=255)

                                                tctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCg                                        >  2:1392068/1‑93 (MQ=255)

                                                  tttGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGt                                      >  1:3399390/1‑93 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                  <  1:940967/93‑1 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                  <  1:3311001/93‑1 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                  <  1:910045/93‑1 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                  <  1:4532189/93‑1 (MQ=255)

                                                       tttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGAc                                 <  1:1621235/93‑1 (MQ=255)

                                                        ttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACa                                >  1:3710338/1‑93 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                              <  2:1606790/93‑1 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                              <  2:4316604/93‑1 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                              <  2:3055342/93‑1 (MQ=255)

                                                             cGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGAt                           <  1:1878558/93‑1 (MQ=255)

                                                              gCCTTCCGCCACGACGGTTACATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                          <  1:3128713/93‑1 (MQ=255)

                                                              gCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                          <  1:1587105/93‑1 (MQ=255)

                                                               ccTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCg                         >  1:1205180/1‑93 (MQ=255)

                                                                cTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                        >  2:1397444/1‑93 (MQ=255)

                                                                cTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                        >  1:2163357/1‑93 (MQ=255)

                                                                 ttCCGCCTCGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                       >  1:314744/1‑93 (MQ=255)

                                                                   ccGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACaa                     <  2:2833630/93‑1 (MQ=255)

                                                                       cACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACt                 <  1:1414781/93‑1 (MQ=255)

                                                                        acgacgGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtt                <  2:3248239/93‑1 (MQ=255)

                                                                         cgacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACttt               <  2:1721025/93‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt              <  2:2584908/93‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt              <  2:2306576/93‑1 (MQ=255)

                                                                           atggTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTg             <  1:805469/91‑1 (MQ=255)

                                                                           acgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTg             >  1:3887617/1‑93 (MQ=255)

                                                                           acgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTg             >  2:3710967/1‑93 (MQ=255)

                                                                           acgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTGCGGTCAGACAACGATCGACAAAACTTTTg             <  2:2888877/93‑1 (MQ=255)

                                                                             ggTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTc           <  1:2432346/93‑1 (MQ=255)

                                                                             ggTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTc           <  1:4486330/93‑1 (MQ=255)

                                                                               ttAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa         <  1:3579553/93‑1 (MQ=255)

                                                                               ttAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa         <  2:72173/93‑1 (MQ=255)

                                                                               gtAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa         <  1:3104122/92‑1 (MQ=255)

                                                                                tAAATTGAGGTTTTGCTCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAt        <  1:4183195/93‑1 (MQ=255)

                                                                                tAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAt        >  2:1593217/1‑93 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt       >  1:2754052/1‑93 (MQ=255)

                                                                                    ttGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTg    <  1:26849/93‑1 (MQ=255)

                                                                                    ttGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTg    >  2:4324755/1‑93 (MQ=255)

                                                                                     tGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGa   >  2:3823203/1‑93 (MQ=255)

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAT  >  NC_000913/1561431‑1561609

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