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	98.8%	613.3 / ‑5.1	173	P130S (CCA→TCA)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (80/91); total (80/93)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 174 total aligned reads displayed.

AGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGGTCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGA  >  NC_000913/1723179‑1723358

aGTTCGATACGCAAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATc                                                                                         <  2:6935944/93‑1 (MQ=255)

 gTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCg                                                                                        >  2:3239463/1‑93 (MQ=255)

 gTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCg                                                                                        <  2:2850172/93‑1 (MQ=255)

 gTTCGATACGCCAAAATGACGCACTTTCCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCg                                                                                        <  1:2197124/93‑1 (MQ=255)

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

    cGATACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTg                                                                                     <  2:7837634/93‑1 (MQ=255)

      aTACGCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGa                                                                                   >  1:1841566/1‑93 (MQ=255)

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

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

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

          gCCAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTaa                                                                               >  1:7103320/1‑93 (MQ=255)

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

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

           ccAAAATGACGCTCTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAAc                                                                              <  2:6294121/93‑1 (MQ=255)

            cAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACa                                                                             >  2:726810/1‑93 (MQ=255)

            cAAAATGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACa                                                                             <  2:6683047/93‑1 (MQ=255)

             aaaaTGACGCACTTTGCCGTTCTGATGCAGATTTTTGAACGCGTTCGCCACTTCATCGGCATCCTTTAACGGGGCTTATCGGTGGATTAACAg                                                                            >  1:7004552/1‑93 (MQ=255)

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

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

             aaaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAg                                                                            <  2:260155/93‑1 (MQ=255)

             aaaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAg                                                                            <  2:3840631/93‑1 (MQ=255)

              aaaTGGCGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGc                                                                           <  2:5362277/93‑1 (MQ=255)

              aaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGc                                                                           >  1:5465167/1‑93 (MQ=255)

              aaaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGc                                                                           <  1:3588738/93‑1 (MQ=255)

               aaTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCa                                                                          >  1:5594441/1‑93 (MQ=255)

                aTGACGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCaa                                                                         >  1:7810490/1‑93 (MQ=255)

                   aCGCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGCGGCTTAACAGCAAATc                                                                      >  1:1060980/1‑93 (MQ=255)

                     gCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCa                                                                    <  1:7347490/93‑1 (MQ=255)

                     gCACTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCa                                                                    <  2:4580353/93‑1 (MQ=255)

                       ccTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  <  1:2116736/92‑1 (MQ=255)

                       aCTTTGCCTCTGTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  <  1:2661379/93‑1 (MQ=255)

                       aCTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  >  2:1443524/1‑93 (MQ=255)

                       aCTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  <  1:1434436/93‑1 (MQ=255)

                       aCTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  <  2:3170804/93‑1 (MQ=255)

                       aCTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAga                                                                  >  2:732046/1‑93 (MQ=255)

                        cTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAgat                                                                 <  2:1749969/93‑1 (MQ=255)

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

                        cTTTGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCAAgat                                                                 >  1:1585390/1‑93 (MQ=255)

                          ttGCCGCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAgatga                                                               >  1:2911124/1‑93 (MQ=255)

                               gCTCTGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGt                                                          >  1:1321473/1‑93 (MQ=255)

                                ctctGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCTGATGATCGGTc                                                         <  1:1681709/93‑1 (MQ=255)

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

                                 tctGATGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTcg                                                        <  1:7685525/93‑1 (MQ=255)

                                      tGCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCgaga                                                   >  1:3011117/1‑93 (MQ=255)

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

                                       gCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGat                                                  >  1:3309094/1‑93 (MQ=255)

                                       gCAGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGat                                                  >  2:4921219/1‑93 (MQ=255)

                                         aGATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGatta                                                <  1:4972060/93‑1 (MQ=255)

                                          gATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaa                                               <  1:6466803/93‑1 (MQ=255)

                                          gATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaa                                               >  1:2344115/1‑93 (MQ=255)

                                          gATGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaa                                               <  1:3533289/93‑1 (MQ=255)

                                           aTGTTTGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaat                                              <  2:653422/93‑1 (MQ=255)

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

                                             ttttGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGattaatta                                            <  2:5593639/92‑1 (MQ=255)

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

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

                                              tttGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAg                                           <  2:5650989/93‑1 (MQ=255)

                                              tttGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAg                                           >  2:2186821/1‑93 (MQ=255)

                                               ttGAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGc                                          <  1:6933446/93‑1 (MQ=255)

                                                 gAACGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGa                                        <  1:3053581/93‑1 (MQ=255)

                                                  aaCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGAc                                       <  2:1262019/93‑1 (MQ=255)

                                                   aCGCGTCCGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACt                                      >  2:3985463/1‑93 (MQ=255)

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

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

                                                          cGCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgc                               >  2:5435752/1‑93 (MQ=255)

                                                           gCCACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGACGATCGGTCGCGAGCTTAATTAGCGACTGTTCGgcg                              <  1:1400346/93‑1 (MQ=255)

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

                                                            ccACTTCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGgcgc                             <  2:3897986/93‑1 (MQ=255)

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

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

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

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

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

                                                                ttCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCtt                         <  2:4260180/93‑1 (MQ=255)

                                                                ttCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCtt                         <  1:1706448/93‑1 (MQ=255)

                                                                ttCATCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCtt                         >  2:5480258/1‑93 (MQ=255)

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

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

                                                                    tCGGCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatg                     >  2:5766527/1‑93 (MQ=255)

                                                                       gCATCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAatgatg                  <  2:7379444/93‑1 (MQ=255)

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

                                                                        cATCCATTAACGGGTCTGGTCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGAtgt                 <  2:2649213/93‑1 (MQ=255)

                                                                          tCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGa               >  2:4380112/1‑93 (MQ=255)

                                                                          tCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGa               >  2:914492/1‑93 (MQ=255)

                                                                          tCCATTAACGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGa               >  2:2196807/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                                aaCGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCgg         <  2:6958137/93‑1 (MQ=255)

                                                                                  cGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTCATTAGCGACTGTTCGGCGCTCTTAATGATGTGATcgcgcga       >  1:1891489/1‑90 (MQ=255)

                                                                                  cGGGTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTc       >  2:5206468/1‑93 (MQ=255)

                                                                                    ggTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAg     <  2:5170508/93‑1 (MQ=255)

                                                                                     gTCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTCATGATGTGATCGCGGTCAGt    >  2:5833398/1‑93 (MQ=255)

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

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

                                                                                      tCTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTg   <  2:134715/93‑1 (MQ=255)

                                                                                       cTGATCGGTGGATTAACAGCAAATCCAGATGATCGGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGa  <  1:4416277/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` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`