BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,723,269	G→A	P130S (CCA→TCA)	b1647 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,723,269	0	G	A	99.1%	808.5 / ‑4.3	228	P130S (CCA→TCA)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (101/124); total (104/124)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 235 total aligned reads displayed.

AGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGGTCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGA  >  NC_000913/1723179‑1723358

aGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATc                                                                                         <  1:967815/93‑1 (MQ=255)

aGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATc                                                                                         <  2:245150/93‑1 (MQ=255)

  ttCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCgg                                                                                       >  1:7950207/1‑93 (MQ=255)

  ttCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCgg                                                                                       <  1:1429162/93‑1 (MQ=255)

  ttCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCgg                                                                                       <  2:8805205/93‑1 (MQ=255)

  tgcGATGCGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGCCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCgg                                                                                       <  1:1246265/91‑1 (MQ=255)

   tCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGt                                                                                      <  1:9142391/93‑1 (MQ=255)

   tCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGt                                                                                      <  1:2162209/93‑1 (MQ=255)

     gATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGttg                                                                                    >  2:3452128/1‑91 (MQ=255)

       tACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGGTCGGTGGAt                                                                                  >  2:1299265/1‑93 (MQ=255)

        aCGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGAtt                                                                                 <  1:1951933/93‑1 (MQ=255)

        aCGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGAtt                                                                                 <  1:6949435/93‑1 (MQ=255)

        aCGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGAtt                                                                                 <  2:2638947/93‑1 (MQ=255)

         cGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTa                                                                                <  1:3820483/93‑1 (MQ=255)

          gCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTaa                                                                               <  1:1895871/93‑1 (MQ=255)

             aaaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAg                                                                            >  2:6424702/1‑93 (MQ=255)

             aaaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGAGTAACac                                                                            >  1:734303/1‑92 (MQ=255)

               aaTGACGCCCTTTGCCGCTCTGATGCAGATTTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCa                                                                          <  1:4466425/93‑1 (MQ=255)

                  gACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAAt                                                                       >  2:9302466/1‑93 (MQ=255)

                  gACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAAt                                                                       >  2:9383546/1‑93 (MQ=255)

                   aCGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATc                                                                      <  1:7551813/93‑1 (MQ=255)

                   aCGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATc                                                                      <  2:9541206/93‑1 (MQ=255)

                   aCGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATc                                                                      <  2:3901583/93‑1 (MQ=255)

                    cGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATcc                                                                     <  1:4521935/93‑1 (MQ=255)

                     gCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCa                                                                    >  1:3154254/1‑93 (MQ=255)

                      cACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAg                                                                   <  2:4882614/93‑1 (MQ=255)

                      cACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAg                                                                   >  2:5769847/1‑93 (MQ=255)

                      cACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAg                                                                   >  1:9939322/1‑93 (MQ=255)

                        cTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAgat                                                                 >  2:5310229/1‑93 (MQ=255)

                         tttGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAgatg                                                                <  1:7094444/93‑1 (MQ=255)

                          ttGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAgatga                                                               <  2:6059210/93‑1 (MQ=255)

                            gCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATc                                                             >  2:5025629/1‑93 (MQ=255)

                            gCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATc                                                             <  2:6150198/93‑1 (MQ=255)

                            gCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATc                                                             >  2:7753261/1‑93 (MQ=255)

                               gCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGt                                                          <  2:7421310/93‑1 (MQ=255)

                                ctctGATGCAGATGTTTGAACGCATCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATACACATGATCgggc                                                         >  2:5612465/1‑91 (MQ=255)

                                 tctGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTcg                                                        >  1:2735200/1‑93 (MQ=255)

                                   tGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTcgcg                                                      <  2:1634456/93‑1 (MQ=255)

                                    gATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCga                                                     <  1:6270577/93‑1 (MQ=255)

                                    gATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCga                                                     <  2:1186858/93‑1 (MQ=255)

                                     aTGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCgag                                                    >  2:9510575/1‑93 (MQ=255)

                                      tGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCgaga                                                   <  1:4292780/93‑1 (MQ=255)

                                         aGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGatta                                                >  2:5716236/1‑93 (MQ=255)

                                           aTGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaat                                              >  1:5588287/1‑93 (MQ=255)

                                            tGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatt                                             >  1:7855403/1‑93 (MQ=255)

                                            tGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatt                                             >  2:323975/1‑93 (MQ=255)

                                             gTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatta                                            <  1:8811741/93‑1 (MQ=255)

                                             gTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatta                                            <  1:1010707/93‑1 (MQ=255)

                                             gTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatta                                            <  1:174428/93‑1 (MQ=255)

                                              tttGACCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAg                                           <  2:9846887/93‑1 (MQ=255)

                                              tgtaaaCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCTCGAGATTAATTAg                                           <  2:5149666/89‑1 (MQ=255)

                                               ttaaaCGCGTCCGCCACTTCATCGGCATCCATTATCGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGc                                          <  1:309875/90‑1 (MQ=255)

                                                  aaCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGAc                                       >  1:8495471/1‑93 (MQ=255)

                                                   aCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGCATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACt                                      <  1:5045478/93‑1 (MQ=255)

                                                    cgcgTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTg                                     <  2:2986875/93‑1 (MQ=255)

                                                      cgTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGtt                                   <  2:7526315/93‑1 (MQ=255)

                                                       gcccGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTc                                  <  2:9008296/91‑1 (MQ=255)

                                                       gTCCCCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTc                                  <  2:2662633/93‑1 (MQ=255)

                                                         ccGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCgg                                >  2:6659835/1‑93 (MQ=255)

                                                          cGCCACTTCATCGGCATCCATTTACGGGTCTGTTCGGTGGTTTAACAGCAAATCCAGATGATTGTTTGCGAGATTAATTAGCGACCGTTCGgc                               >  2:9264081/1‑93 (MQ=255)

                                                          cGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgc                               <  1:7454159/93‑1 (MQ=255)

                                                          cGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgc                               <  1:3657721/93‑1 (MQ=255)

                                                           gCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGGGAGATTAATTAGCGACTGTTCGgcg                              >  2:5867856/1‑93 (MQ=255)

                                                           gCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgcg                              >  2:9512242/1‑93 (MQ=255)

                                                            ccACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgcgc                             >  1:7676832/1‑93 (MQ=255)

                                                              ccTTTATCGGCAACCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctc                           <  1:8184095/92‑1 (MQ=255)

                                                              aCTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctt                           >  2:2801862/1‑92 (MQ=255)

                                                              aCTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctc                           >  2:2747512/1‑93 (MQ=255)

                                                              aCTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctc                           >  1:9865426/1‑93 (MQ=255)

                                                               cTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctct                          <  1:338031/93‑1 (MQ=255)

                                                               cTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctct                          >  1:2597342/1‑93 (MQ=255)

                                                               cTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGctct                          >  2:5692223/1‑93 (MQ=255)

                                                                 tCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTa                        <  2:1129788/93‑1 (MQ=255)

                                                                  cATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTaa                       >  1:29042/1‑93 (MQ=255)

                                                                  cATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTaa                       <  1:5364310/93‑1 (MQ=255)

                                                                  cATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTaa                       >  2:4667223/1‑93 (MQ=255)

                                                                   aTCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAat                      <  1:4689612/93‑1 (MQ=255)

                                                                   aTCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAat                      <  2:7245471/93‑1 (MQ=255)

                                                                     cGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatga                    <  1:2846665/93‑1 (MQ=255)

                                                                     cGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatga                    <  2:2321786/93‑1 (MQ=255)

                                                                      ggCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatgat                   >  2:8204114/1‑93 (MQ=255)

                                                                       gCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatgatg                  >  2:543751/1‑93 (MQ=255)

                                                                        cATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGAtgt                 >  2:7340727/1‑93 (MQ=255)

                                                                        cATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGAtgt                 <  1:4664243/93‑1 (MQ=255)

                                                                           ccATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGAt              <  2:2342701/93‑1 (MQ=255)

                                                                            cATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATc             >  2:9810475/1‑93 (MQ=255)

                                                                            cATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATc             <  2:1638338/93‑1 (MQ=255)

                                                                             aTTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcg            >  1:9347853/1‑93 (MQ=255)

                                                                             aTTAACGGGTCTGACCGGTGGATTAACAACAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcg            <  2:3302106/93‑1 (MQ=255)

                                                                              ttAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcgc           <  2:2077208/93‑1 (MQ=255)

                                                                              ttAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcgc           >  1:4987048/1‑93 (MQ=255)

                                                                               tAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcgcg          <  2:1893503/93‑1 (MQ=255)

                                                                                 aCGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGt        <  1:659625/93‑1 (MQ=255)

                                                                                   gggTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCACGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCa      >  1:4033785/1‑93 (MQ=255)

                                                                                     gTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGt    <  1:6369427/93‑1 (MQ=255)

                                                                                      tCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTg   >  2:4276020/1‑93 (MQ=255)

                                                                                      tCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTg   >  2:798457/1‑93 (MQ=255)

                                                                                      tCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTg   >  2:6838855/1‑93 (MQ=255)

                                                                                      tCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTg   >  2:5843427/1‑93 (MQ=255)

                                                                                       cTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGa  <  1:71945/93‑1 (MQ=255)

AGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGGTCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGA  >  NC_000913/1723179‑1723358

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