BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,354,373	G→A	G254G (GGC→GGT)	sapA ←	SapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,354,373	0	G	A	99.5%	702.2 / NA	205	G254G (GGC→GGT)	sapA	SapA
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (100/104); minor base C (0/1); total (100/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGC  >  NC_000913/1354283‑1354462

gtcgtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccacc                                                                                         >  1:5168895/1‑93 (MQ=255)

  cgtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCgg                                                                                       >  2:506451/1‑93 (MQ=255)

   gtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGa                                                                                      <  2:5393761/93‑1 (MQ=255)

     cACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGc                                                                                    <  2:6841399/93‑1 (MQ=255)

      aCGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGcc                                                                                   <  2:5368622/93‑1 (MQ=255)

            aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             >  2:3659634/1‑93 (MQ=255)

            aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             <  1:5754205/93‑1 (MQ=255)

              tGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATcca                                                                           >  2:1549706/1‑93 (MQ=255)

                gATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATccacc                                                                         >  1:3096778/1‑93 (MQ=255)

                 aTAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATccacca                                                                        <  2:6121847/93‑1 (MQ=255)

                  tAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCcaccac                                                                       <  2:4559473/93‑1 (MQ=255)

                   aGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACt                                                                      >  1:6504198/1‑93 (MQ=255)

                    gctggctgGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTa                                                                     <  1:829054/93‑1 (MQ=255)

                      tggctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGGCACCGGAGCCTAAATCCACCACTAcc                                                                   >  2:5583080/1‑93 (MQ=255)

                        gctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTg                                                                 <  2:2372203/93‑1 (MQ=255)

                        gctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTg                                                                 <  2:7766961/93‑1 (MQ=255)

                          tggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCg                                                               >  1:1088153/1‑93 (MQ=255)

                          tggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCg                                                               <  2:1148862/93‑1 (MQ=255)

                           ggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCgg                                                              <  2:217450/93‑1 (MQ=255)

                           ggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCgg                                                              <  2:4183245/93‑1 (MQ=255)

                            gcagcaGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGc                                                             <  2:342601/93‑1 (MQ=255)

                             cagcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCa                                                            >  2:322438/1‑93 (MQ=255)

                              agcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAt                                                           >  2:1458720/1‑93 (MQ=255)

                              agcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAt                                                           <  2:543032/93‑1 (MQ=255)

                               gcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAtt                                                          >  1:3824899/1‑93 (MQ=255)

                                caggccaggcCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTa                                                         <  1:1270363/93‑1 (MQ=255)

                                 aggccaggccAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTaa                                                        <  2:251915/93‑1 (MQ=255)

                                  ggccaggccaGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAAc                                                       <  2:2283915/93‑1 (MQ=255)

                                  ggccaggccaGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAAc                                                       >  1:5667877/1‑93 (MQ=255)

                                   gccaggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACg                                                      <  2:3890590/93‑1 (MQ=255)

                                   gccaggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACg                                                      <  2:4405041/93‑1 (MQ=255)

                                      aggcctgAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGCCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGtt                                                   >  1:4643763/1‑93 (MQ=255)

                                      aggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGtt                                                   >  1:7798832/1‑93 (MQ=255)

                                      aggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGtt                                                   >  2:5256688/1‑93 (MQ=255)

                                       ggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGttt                                                  <  1:5861213/93‑1 (MQ=255)

                                        gccagAACGTCGCATTCCCCGGTCAGGAGGTTCGACAGACGGCCGGTGCCCCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTa                                                 <  1:6007720/93‑1 (MQ=255)

                                         ccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAc                                                <  1:6905294/93‑1 (MQ=255)

                                          cagAAGGTCGCATTCCACGGGCAGGAGGTTCGACAGACGTCCGGTGCTACCGGAGCCTAAATCTACCACTACCTGCGGCATTAACGGTTTAcc                                               <  2:4562698/93‑1 (MQ=255)

                                          cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               >  2:6911082/1‑93 (MQ=255)

                                          cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               >  2:7394306/1‑93 (MQ=255)

                                          cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               >  2:5742868/1‑93 (MQ=255)

                                            gAACGTCGCATTCCCCGGTCAGGAGTTTGGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  1:7137131/1‑93 (MQ=255)

                                            gAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  2:6540231/1‑93 (MQ=255)

                                            gAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  2:4075643/1‑93 (MQ=255)

                                            gAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  2:7919631/1‑93 (MQ=255)

                                             aaCGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgcg                                            >  2:4327782/1‑93 (MQ=255)

                                              aCGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCgcgc                                           >  1:1571184/1‑93 (MQ=255)

                                                gTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCa                                         >  1:4365378/1‑93 (MQ=255)

                                                  cGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGa                                       <  2:3177516/93‑1 (MQ=255)

                                                     aTTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGt                                    <  1:801164/93‑1 (MQ=255)

                                                      ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTCACGGTTTACCGCTCCAGAAGtc                                   >  1:318387/1‑93 (MQ=255)

                                                      ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  1:4514446/1‑93 (MQ=255)

                                                      ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  2:1299749/1‑93 (MQ=255)

                                                      ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  1:3921433/1‑93 (MQ=255)

                                                      ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  1:3469593/1‑93 (MQ=255)

                                                       tCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtca                                  <  2:3108239/93‑1 (MQ=255)

                                                       tCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtca                                  <  2:5834267/93‑1 (MQ=255)

                                                        ccccGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtcat                                 <  2:1676863/93‑1 (MQ=255)

                                                         cccGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGGCAGAAGtcata                                >  1:2013771/1‑92 (MQ=255)

                                                           cGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTcatcat                              <  2:4632205/93‑1 (MQ=255)

                                                           cGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTcatcat                              <  1:4587459/93‑1 (MQ=255)

                                                            ggTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAAGCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATg                             <  1:2954898/93‑1 (MQ=255)

                                                             gTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGa                            <  1:7173775/93‑1 (MQ=255)

                                                             gTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGa                            <  1:2681718/93‑1 (MQ=255)

                                                              tCAGGATTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           <  2:2965170/93‑1 (MQ=255)

                                                              tCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           <  1:5654927/93‑1 (MQ=255)

                                                              tCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTNANCGCGCCAGAANTCATCATGAc                           >  1:7275089/1‑93 (MQ=255)

                                                              tCAGGAGTTTCGACAGACGCCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           >  2:6713019/1‑93 (MQ=255)

                                                              gcAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           <  1:4941358/92‑1 (MQ=255)

                                                                aGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGt                         <  1:2769357/93‑1 (MQ=255)

                                                                aGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGt                         <  1:7626534/93‑1 (MQ=255)

                                                                 ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTACCGGTTTACCGCGCCAGAAGTCATCATGACGtt                        <  2:5593514/93‑1 (MQ=255)

                                                                 ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        >  2:185922/1‑93 (MQ=255)

                                                                 ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        <  1:4024979/93‑1 (MQ=255)

                                                                 ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        >  1:349298/1‑93 (MQ=255)

                                                                  gAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTg                       <  1:2133591/93‑1 (MQ=255)

                                                                   aGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGt                      <  1:6726537/93‑1 (MQ=255)

                                                                    gTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTa                     <  2:73008/93‑1 (MQ=255)

                                                                     tttCGACAGACGTCCGGTGCCACCGGGGCCTAAATACACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAg                    >  2:73041/1‑93 (MQ=255)

                                                                     tttCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAg                    <  2:4551583/93‑1 (MQ=255)

                                                                      ttCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAgg                   >  1:7581329/1‑93 (MQ=255)

                                                                        cGACAGACGTCCCGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGGAGGCg                 <  1:3721317/93‑1 (MQ=255)

                                                                         gACAGACGCCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGa                >  2:3779572/1‑93 (MQ=255)

                                                                          aCAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGaa               >  2:2099465/1‑93 (MQ=255)

                                                                          aCAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGaa               >  2:1176539/1‑93 (MQ=255)

                                                                           cAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAAt              <  1:5483857/93‑1 (MQ=255)

                                                                            aGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATa             >  2:4014580/1‑93 (MQ=255)

                                                                             gACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATaa            >  2:4729856/1‑93 (MQ=255)

                                                                              aCGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATaaa           >  2:4882714/1‑93 (MQ=255)

                                                                              aCGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATaaa           <  1:7567630/93‑1 (MQ=255)

                                                                              aCGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATaaa           <  1:1214852/93‑1 (MQ=255)

                                                                                gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         >  1:4841097/1‑93 (MQ=255)

                                                                                gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         <  1:7442393/93‑1 (MQ=255)

                                                                                gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         <  1:506451/93‑1 (MQ=255)

                                                                                gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         <  2:4315770/93‑1 (MQ=255)

                                                                                  ccGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGc       <  1:6601577/93‑1 (MQ=255)

                                                                                   cGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGcc      <  2:5884653/93‑1 (MQ=255)

                                                                                   cGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGcc      <  1:4398261/93‑1 (MQ=255)

                                                                                     gTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCg    >  2:2801656/1‑93 (MQ=255)

                                                                                      tGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCgg   <  2:4745868/93‑1 (MQ=255)

GTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGC  >  NC_000913/1354283‑1354462

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