BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	98.1%	1577.1 / 17.0	474	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (7/2); new base C (252/213); total (259/215)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 518 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                           >  2:4953654/1‑93 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cc                                                                                        <  1:7378785/93‑1 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cc                                                                                        <  2:7443886/93‑1 (MQ=255)

   gcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccc                                                                                       <  1:5537956/93‑1 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                      >  1:4403099/1‑92 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                    >  2:2363034/1‑90 (MQ=255)

       cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                   >  2:3501215/1‑89 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑CAaccc                                                                                >  1:433099/1‑93 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑CAaccc                                                                                >  2:6492556/1‑93 (MQ=255)

            ctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCTCTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCAc                                                                              <  1:6857902/93‑1 (MQ=255)

            ctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCAc                                                                              <  2:6693945/93‑1 (MQ=255)

             tgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACg                                                                             >  1:1906545/1‑93 (MQ=255)

             tgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACg                                                                             >  1:2373278/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGcc                                                                           >  1:5771798/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGcc                                                                           >  1:6416836/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCaa                                                                         >  1:3887972/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCaa                                                                         >  1:1541873/1‑93 (MQ=255)

                  tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAAc                                                                        >  1:6621077/1‑93 (MQ=255)

                  tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAAc                                                                        <  2:4864009/93‑1 (MQ=255)

                   cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAAcc                                                                       >  2:4337020/1‑93 (MQ=255)

                    ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCa                                                                      >  2:794325/1‑93 (MQ=255)

                     cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAg                                                                     <  2:2965813/93‑1 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGc                                                                    >  1:3375427/1‑93 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGc                                                                    >  1:1013152/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGca                                                                   >  1:1680609/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGca                                                                   <  2:596971/93‑1 (MQ=255)

                        ccAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGcac                                                                  <  1:4725675/93‑1 (MQ=255)

                         cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGcaca                                                                 >  1:6158303/1‑93 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACaa                                                                <  1:4114877/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAata                                                              <  2:1890838/93‑1 (MQ=255)

                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATg                                                            <  1:36136/93‑1 (MQ=255)

                               aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGa                                                           <  1:1200256/93‑1 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGaa                                                          >  2:1693005/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGaa                                                          >  1:7279558/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGaa                                                          <  1:2949023/93‑1 (MQ=255)

                                   aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCg                                                       <  1:875823/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGc                                                      <  1:423357/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  2:2791609/1‑93 (MQ=255)

                                       cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGAGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAg                                                   <  2:4580911/93‑1 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  <  1:3453185/93‑1 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  2:7013984/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 >  2:2451904/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:2439839/93‑1 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  1:2175158/1‑93 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  2:4383912/93‑1 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  2:5833230/1‑93 (MQ=255)

                                                ccACTTCAACACTTAAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  2:6442895/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:6802808/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  2:7272552/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:7190918/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:927126/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  1:6201327/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  2:5410/1‑93 (MQ=255)

                                                    tctcacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:7114787/5‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  2:5723077/1‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     <  1:6107754/93‑1 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     <  1:4912344/93‑1 (MQ=255)

                                                      caacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                                    >  1:2801995/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:6871657/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:1057125/93‑1 (MQ=255)

                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcccgg                                  <  1:2646943/93‑3 (MQ=255)

                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                 >  2:4182339/1‑93 (MQ=255)

                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                 >  2:2626960/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                <  1:2139162/93‑1 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:2131945/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:2011847/1‑93 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             >  1:2567685/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  2:5328091/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  1:5021628/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:3175224/93‑1 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  1:4623551/93‑1 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:3388017/1‑93 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  1:1743032/93‑1 (MQ=255)

                                                                 aCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACg                         >  1:6558050/1‑93 (MQ=255)

                                                                 aCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACg                         >  2:4712819/1‑93 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        <  2:3982837/93‑1 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        >  1:7025406/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     <  2:6358699/93‑1 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     <  2:509628/93‑1 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:5850325/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  2:2029683/1‑93 (MQ=255)

                                                                         aCTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:1312789/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcggcgcc               >  2:1139792/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               >  2:2496803/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               <  2:1457557/93‑1 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  1:5251243/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:3116493/1‑93 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             <  2:3721046/93‑1 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             >  2:1396456/1‑93 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             <  1:5431711/93‑1 (MQ=255)

                                                                              gCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGGCGCCAccgcgccgcg            >  1:3742748/1‑93 (MQ=255)

                                                                              gCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  1:5666510/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  1:1260389/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  1:3282933/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  2:38424/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         <  2:2746561/93‑1 (MQ=255)

                                                                                  gAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  2:1238050/93‑1 (MQ=255)

                                                                                       tAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   >  2:5570449/1‑93 (MQ=255)

                                                                                        aGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  <  2:6085672/93‑1 (MQ=255)

                                                                                        aGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  <  1:454944/93‑1 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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