Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.

ALE, Flask, Isolate	Predicted Mutations	Mean Coverage	Total Reads	Percent Mapped	Mapped Reads	Average Read Length
A1 F1 I1 R1	1548	369.5	20095359	96.9%	19472402	93.0

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,519	T→C	T122A (ACG→GCG)	yddE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,519	0	T	C	99.6%	795.1 / ‑5.1	249	T122A (ACG→GCG)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (130/118); total (131/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 258 total aligned reads displayed.

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608

aGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAc                                                                                         >  2:3879428/1‑93 (MQ=255)

  ccTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCt                                                                                       <  2:5172532/93‑1 (MQ=255)

   cTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCtt                                                                                      >  1:4893861/1‑93 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    >  1:4522130/1‑93 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    <  2:4944633/93‑1 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    <  2:8458222/93‑1 (MQ=255)

      gCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGtt                                                                                   >  1:179136/1‑93 (MQ=255)

        aGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTcc                                                                                 <  1:6277272/93‑1 (MQ=255)

         gAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCa                                                                                >  2:7157806/1‑93 (MQ=255)

          aaTATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAg                                                                               <  1:5609027/93‑1 (MQ=255)

          aaTATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAg                                                                               <  2:8317780/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCTCCTTCCAGCGGTGGCTCAAATCCCGGCGCACCTTGTTCCAGc                                                                              <  2:2533018/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGc                                                                              <  1:923109/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGc                                                                              >  1:2510511/1‑93 (MQ=255)

            tatCGTCCTCAGTTAGATGATGCTCGTTGATAATCTCCGCTCGTGTTTCACCTTCCAGTGGTGGCTCAAAGCCCGTCGCACCTTGTTCCATCg                                                                             >  1:8940645/1‑93 (MQ=255)

            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                             >  1:3701996/1‑93 (MQ=255)

               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                          <  2:904770/93‑1 (MQ=255)

                gTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAAt                                                                         >  1:8667279/1‑93 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                        <  2:4457307/93‑1 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                        <  1:9657714/93‑1 (MQ=255)

                  ccTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAc                                                                       <  2:2516948/93‑1 (MQ=255)

                   cTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACg                                                                      >  1:6550177/1‑93 (MQ=255)

                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     <  2:4152026/93‑1 (MQ=255)

                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGAt                                                                    <  2:1723221/93‑1 (MQ=255)

                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGCGAAATACGAt                                                                    <  1:3650947/93‑1 (MQ=255)

                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:6633193/1‑93 (MQ=255)

                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:6695418/1‑93 (MQ=255)

                       gTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                  <  2:4385253/93‑1 (MQ=255)

                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 <  2:1188644/93‑1 (MQ=255)

                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 >  2:7818327/1‑93 (MQ=255)

                         gagaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACTATAatc                                                                >  1:7504659/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               >  2:4262074/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               <  1:2564877/93‑1 (MQ=255)

                           gaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcat                                                              >  2:5373350/1‑93 (MQ=255)

                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            <  2:9363140/93‑1 (MQ=255)

                              taaGCGCGTTGATAATCTCCGCCTGTGTTTCACCGTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCGGCGAAATACGATAATCATCGt                                                           <  1:3474488/92‑1 (MQ=255)

                              gAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGt                                                           <  2:9902022/93‑1 (MQ=255)

                               aaGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGtt                                                          <  2:7233092/93‑1 (MQ=255)

                               aaGCGCGTTGATAATCGCCCCACGGGCTTCACCTTCCAGCGGGGGCTCAAAGCCGGGGGCACCTTGTTCACGCGAAATACGTTATTCATAGtt                                                          >  2:848418/1‑93 (MQ=255)

                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         <  1:7656166/93‑1 (MQ=255)

                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         >  1:3842309/1‑93 (MQ=255)

                                aGCGCGTTAATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         <  2:6696697/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                      <  2:5304499/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                      <  1:7083039/93‑1 (MQ=255)

                                    cgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCt                                                     <  2:6555757/93‑1 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  1:2056802/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  2:5751460/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  2:7363229/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   <  1:5928081/93‑1 (MQ=255)

                                        gATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                 >  1:7882945/1‑93 (MQ=255)

                                        gATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                 <  2:7636255/93‑1 (MQ=255)

                                          tAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCTCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCgt                                               >  1:5274776/1‑92 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGCTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                              >  1:6867062/1‑93 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:6987802/93‑1 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:7454279/93‑1 (MQ=255)

                                               cccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                          <  1:6840248/92‑1 (MQ=255)

                                                 cgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAc                                        >  2:2692894/1‑93 (MQ=255)

                                                    aCGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgcg                                     >  1:727894/1‑93 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgcca                                    >  2:1939619/1‑91 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:1605234/93‑1 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:8912745/93‑1 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    >  2:4731064/1‑93 (MQ=255)

                                                      gtgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAt                                   >  2:5412496/1‑93 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  >  1:6397984/1‑93 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  <  2:7957714/93‑1 (MQ=255)

                                                       tgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATg                                  >  1:8782451/1‑93 (MQ=255)

                                                        gtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGt                                 >  2:1214351/1‑93 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                <  1:9798985/93‑1 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                <  2:3123535/93‑1 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                                >  1:4406163/1‑93 (MQ=255)

                                                              ccTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCt                           <  1:2308866/93‑1 (MQ=255)

                                                               cTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTg                          >  2:9692743/1‑93 (MQ=255)

                                                                ttCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGc                         >  2:3690960/1‑93 (MQ=255)

                                                                 tCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                        >  1:2798815/1‑93 (MQ=255)

                                                                 tCCAGCGGGGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                        <  1:2695225/93‑1 (MQ=255)

                                                                   cAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAg                      >  1:3211488/1‑93 (MQ=255)

                                                                     gCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                    >  2:664221/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   >  1:7155764/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   >  1:9112649/1‑93 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                   <  2:5882691/93‑1 (MQ=255)

                                                                       ggTGGCTCAAAGCCCTGCGCGCCTTGTTCCAGCAAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAt                  <  2:8188586/93‑1 (MQ=255)

                                                                       ggTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAt                  >  1:1647252/1‑93 (MQ=255)

                                                                        gTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATg                 >  1:6774947/1‑93 (MQ=255)

                                                                        gTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATg                 >  1:3370984/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  2:2795712/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  2:1442592/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              >  1:467850/1‑93 (MQ=255)

                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGgtt              >  2:2405970/1‑93 (MQ=255)

                                                                              cAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGcc           >  1:9563859/1‑93 (MQ=255)

                                                                              cAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGcc           >  2:3620577/1‑93 (MQ=255)

                                                                               aaaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCa          <  2:344021/93‑1 (MQ=255)

                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg         >  1:9486087/1‑93 (MQ=255)

                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg         <  1:829393/93‑1 (MQ=255)

                                                                                 aGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGa        <  1:3878966/93‑1 (MQ=255)

                                                                                    ccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCg     >  1:4407808/1‑93 (MQ=255)

                                                                                    ccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCg     <  1:3663094/93‑1 (MQ=255)

                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt    >  1:839197/1‑93 (MQ=255)

                                                                                      ggCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg   <  1:5107383/93‑1 (MQ=255)

                                                                                      cgcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg   <  2:3195599/92‑1 (MQ=255)

                                                                                       gcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGc  <  1:7912434/93‑1 (MQ=255)

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

AAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATT  >  NC_000913/1533428‑1533612

AAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAC                                                                                             >  SRR4302540.7777081/1‑94 (MQ=60)

AAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAC                                                                                             >  SRR4302540.8260277/1‑94 (MQ=60)

 AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACC                                                                                            <  SRR4302540.3774950/94‑1 (MQ=60)

  aaaaaaagttcctgccggaagggaaggggcggcggtttgaaaatTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCT                                                                                           <  SRR4302540.8175493/50‑1 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTT                                                                                          >  SRR4302540.961267/1‑94 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTT                                                                                          >  SRR4302540.363523/1‑94 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTT                                                                                          <  SRR4302540.3898619/94‑1 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTT                                                                                          >  SRR4302540.4946491/1‑94 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTT                                                                                          >  SRR4302540.6693366/1‑94 (MQ=60)

   CCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAACCCCGGCGCACCTT                                                                                          >  SRR4302540.2092692/1‑94 (MQ=60)

    CTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTG                                                                                         >  SRR4302540.8647837/1‑94 (MQ=60)

     TGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGT                                                                                        >  SRR4302540.4571104/1‑94 (MQ=60)

     TGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGT                                                                                        >  SRR4302540.7579671/1‑94 (MQ=60)

     TGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGT                                                                                        <  SRR4302540.3837121/94‑1 (MQ=60)

     TGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGT                                                                                        <  SRR4302540.8860112/94‑1 (MQ=60)

      GGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTT                                                                                       <  SRR4302540.7575037/94‑1 (MQ=60)

      GGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTT                                                                                       <  SRR4302540.1749507/94‑1 (MQ=60)

      GGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTT                                                                                       >  SRR4302540.180911/1‑94 (MQ=60)

        CAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCC                                                                                     <  SRR4302540.7068649/94‑1 (MQ=60)

        CAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCC                                                                                     <  SRR4302540.5378480/94‑1 (MQ=60)

         AGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCA                                                                                    >  SRR4302540.985440/1‑94 (MQ=60)

         AGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCA                                                                                    >  SRR4302540.6775317/1‑94 (MQ=60)

         AGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCA                                                                                    <  SRR4302540.6347105/94‑1 (MQ=60)

           AATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGC                                                                                  <  SRR4302540.5670732/94‑1 (MQ=60)

           AATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGC                                                                                  >  SRR4302540.2537625/1‑94 (MQ=60)

            ATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCG                                                                                 <  SRR4302540.932842/94‑1 (MQ=60)

            ATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCG                                                                                 >  SRR4302540.3742407/1‑94 (MQ=60)

             TATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGA                                                                                >  SRR4302540.3253900/1‑94 (MQ=60)

             TATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGA                                                                                >  SRR4302540.4649849/1‑94 (MQ=60)

                CGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAAT                                                                             >  SRR4302540.5146207/1‑94 (MQ=60)

                CGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAAT                                                                             >  SRR4302540.8767958/1‑94 (MQ=60)

                 GTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCGCGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATA                                                                            <  SRR4302540.8961501/94‑1 (MQ=60)

                 GTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATA                                                                            >  SRR4302540.2722574/1‑94 (MQ=60)

                  TCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAC                                                                           <  SRR4302540.4640769/94‑1 (MQ=60)

                  TCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAC                                                                           <  SRR4302540.9771404/94‑1 (MQ=60)

                   CCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACG                                                                          >  SRR4302540.6622831/1‑94 (MQ=60)

                   CCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGGGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTTCAGCGAAATACG                                                                          >  SRR4302540.7201945/1‑94 (MQ=60)

                    CTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGA                                                                         >  SRR4302540.257251/1‑94 (MQ=60)

                    CTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGA                                                                         >  SRR4302540.8612186/1‑94 (MQ=60)

                    CTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGA                                                                         >  SRR4302540.347513/1‑94 (MQ=60)

                      CAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATA                                                                       >  SRR4302540.6706608/1‑94 (MQ=60)

                      CAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATA                                                                       >  SRR4302540.6659455/1‑94 (MQ=60)

                      CAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATA                                                                       >  SRR4302540.6769334/1‑94 (MQ=60)

                      CAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCAAGCGAAATACGATA                                                                       <  SRR4302540.3690939/94‑1 (MQ=60)

                         TGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACTATAATC                                                                    >  SRR4302540.7588534/1‑94 (MQ=60)

                         TGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATC                                                                    <  SRR4302540.1985666/94‑1 (MQ=60)

                         TGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATA                                                                    >  SRR4302540.4150903/1‑94 (MQ=60)

                           AGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCAT                                                                  <  SRR4302540.2592527/94‑1 (MQ=60)

                             ATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCG                                                                >  SRR4302540.245139/1‑94 (MQ=60)

                               TAAGCGCGTTGATAATCTCCGCCTGTGTTTCACCGTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCGGCGAAATACGATAATCATCGTT                                                              <  SRR4302540.3512322/94‑1 (MQ=60)

                               GAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTT                                                              >  SRR4302540.6420761/1‑94 (MQ=60)

                               GAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTT                                                              >  SRR4302540.4539486/1‑94 (MQ=60)

                                AAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTG                                                             >  SRR4302540.3884095/1‑94 (MQ=60)

                                 AGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGT                                                            <  SRR4302540.5425649/94‑1 (MQ=60)

                                 AGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGT                                                            <  SRR4302540.7742368/94‑1 (MQ=60)

                                 AGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGT                                                            <  SRR4302540.9810542/94‑1 (MQ=60)

                                  GCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTG                                                           <  SRR4302540.5198473/94‑1 (MQ=60)

                                    GCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCT                                                         <  SRR4302540.7161871/94‑1 (MQ=60)

                                     CGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTT                                                        <  SRR4302540.2651396/94‑1 (MQ=60)

                                      GTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGACATACGATACTCATCGTTGTGCTTT                                                       >  SRR4302540.6726107/1‑94 (MQ=60)

                                      GTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTT                                                       >  SRR4302540.2078943/1‑94 (MQ=60)

                                      GTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTT                                                       >  SRR4302540.2374228/1‑94 (MQ=60)

                                      GTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTT                                                       >  SRR4302540.9290738/1‑94 (MQ=60)

                                       TTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTT                                                      <  SRR4302540.5993582/94‑1 (MQ=60)

                                       TTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATAGTTGTGCTTTT                                                      <  SRR4302540.9736812/94‑1 (MQ=60)

                                        TGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTC                                                     >  SRR4302540.7974988/1‑94 (MQ=60)

                                        TGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTC                                                     >  SRR4302540.7974989/1‑94 (MQ=60)

                                         GATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCG                                                    <  SRR4302540.3263094/94‑1 (MQ=60)

                                          ATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCTCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGT                                                   >  SRR4302540.5332366/1‑94 (MQ=60)

                                           TAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAT                                                  >  SRR4302540.6943184/1‑94 (MQ=60)

                                           TAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAT                                                  >  SRR4302540.9131287/1‑94 (MQ=60)

                                           TAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGCTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAT                                                  >  SRR4302540.6943183/1‑94 (MQ=60)

                                             ATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAG                                                <  SRR4302540.7537385/94‑1 (MQ=60)

                                             ATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAG                                                <  SRR4302540.7537384/94‑1 (MQ=60)

                                             ATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAG                                                <  SRR4302540.7065381/94‑1 (MQ=60)

                                                GCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCA                                             >  SRR4302540.7028327/1‑94 (MQ=60)

                                                GCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCA                                             <  SRR4302540.3672131/94‑1 (MQ=60)

                                                GCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCA                                             <  SRR4302540.6916010/94‑1 (MQ=60)

                                                CCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCA                                             <  SRR4302540.6916011/94‑1 (MQ=60)

                                                 CCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAC                                            <  SRR4302540.1898987/94‑1 (MQ=60)

                                                 CCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAC                                            <  SRR4302540.4008737/94‑1 (MQ=60)

                                                  CGCACGTGTTTCACCTTCCAGCGGTGGCTCCAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATAATCGTTGTGGTTTTCGATAGTCACG                                           >  SRR4302540.8963709/1‑94 (MQ=60)

                                                    CACGTGTTTCCCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCG                                         <  SRR4302540.48138/94‑1 (MQ=60)

                                                    CACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCG                                         >  SRR4302540.735691/1‑94 (MQ=60)

                                                     ACGTGTTTCACCTTCGAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGA                                        >  SRR4302540.885262/1‑94 (MQ=60)

                                                      CGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAT                                       <  SRR4302540.6920404/94‑1 (MQ=60)

                                                      CGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAT                                       <  SRR4302540.6920405/94‑1 (MQ=60)

                                                       GTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATG                                      >  SRR4302540.8884120/1‑94 (MQ=60)

                                                       GTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATG                                      >  SRR4302540.6469211/1‑94 (MQ=60)

                                                         GTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTT                                    >  SRR4302540.4453931/1‑94 (MQ=60)

                                                          TTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTT                                   <  SRR4302540.9913729/94‑1 (MQ=60)

                                                           TTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTT                                  >  SRR4302540.3714634/1‑94 (MQ=60)

                                                           TTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTT                                  >  SRR4302540.3714632/1‑94 (MQ=60)

                                                           TTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTT                                  >  SRR4302540.7693771/1‑94 (MQ=60)

                                                             CACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCC                                <  SRR4302540.5837221/94‑1 (MQ=60)

                                                             CACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCC                                >  SRR4302540.622754/1‑94 (MQ=60)

                                                               CCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTG                              <  SRR4302540.2333720/94‑1 (MQ=60)

                                                               CCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTG                              <  SRR4302540.7503213/94‑1 (MQ=60)

                                                                 TTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCC                            >  SRR4302540.2829037/1‑94 (MQ=60)

                                                                 TTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCC                            >  SRR4302540.2829038/1‑94 (MQ=60)

                                                                  TCCAGCGGGGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCA                           <  SRR4302540.2724325/94‑1 (MQ=60)

                                                                   CCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAG                          >  SRR4302540.3246351/1‑94 (MQ=60)

                                                                    CAGCGGTGTCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGC                         <  SRR4302540.6212049/94‑1 (MQ=60)

                                                                     AGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCG                        <  SRR4302540.8487893/94‑1 (MQ=60)

                                                                      GCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGCAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAcgga                       >  SRR4302540.7235517/1‑90 (MQ=60)

                                                                      GCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGA                       >  SRR4302540.9222169/1‑94 (MQ=60)

                                                                      GCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGA                       >  SRR4302540.7235516/1‑94 (MQ=60)

                                                                       CGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAT                      <  SRR4302540.2812667/94‑1 (MQ=60)

                                                                       CGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAT                      >  SRR4302540.1665227/1‑94 (MQ=60)

                                                                       CGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGAT                      <  SRR4302540.7843909/94‑1 (MQ=60)

                                                                        GGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATT                     >  SRR4302540.6527335/1‑94 (MQ=60)

                                                                        GGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATG                     >  SRR4302540.6849820/1‑94 (MQ=60)

                                                                        GGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATG                     >  SRR4302540.3407835/1‑94 (MQ=60)

                                                                        GGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATG                     >  SRR4302540.7203178/1‑94 (MQ=60)

                                                                         GTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGT                    >  SRR4302540.9621581/1‑94 (MQ=60)

                                                                          TGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGCTAGTCACGCGCTGTTTTCCTGCCAGCGATGTT                   >  SRR4302540.6804501/1‑94 (MQ=60)

                                                                           GGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTT                  >  SRR4302540.472854/1‑94 (MQ=60)

                                                                              TCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCC               >  SRR4302540.9676906/1‑94 (MQ=60)

                                                                               CAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCA              <  SRR4302540.3146742/94‑1 (MQ=60)

                                                                                AAAGCGCGGCGCACCTTGTTCCAGCGAAATCCGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG             <  SRR4302540.8269078/94‑1 (MQ=60)

                                                                                AAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG             <  SRR4302540.3694161/94‑1 (MQ=60)

                                                                                AAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAG             >  SRR4302540.9598556/1‑94 (MQ=60)

                                                                                 AAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGA            <  SRR4302540.838210/94‑1 (MQ=60)

                                                                                 AAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGA            <  SRR4302540.5152141/94‑1 (MQ=60)

                                                                                  AGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGAT           <  SRR4302540.3921120/94‑1 (MQ=60)

                                                                                    CCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCG         >  SRR4302540.4455597/1‑94 (MQ=60)

                                                                                    CCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCG         >  SRR4302540.6349558/1‑94 (MQ=60)

                                                                                    CCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCG         >  SRR4302540.4455596/1‑94 (MQ=60)

                                                                                     CCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGT        <  SRR4302540.9136753/94‑1 (MQ=60)

                                                                                     CCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGT        <  SRR4302540.3703185/94‑1 (MQ=60)

                                                                                     CCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGT        <  SRR4302540.4531676/94‑1 (MQ=60)

                                                                                     CCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGT        >  SRR4302540.848098/1‑94 (MQ=60)

                                                                                       GGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC      >  SRR4302540.9399728/1‑94 (MQ=60)

                                                                                       GGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC      <  SRR4302540.5162670/94‑1 (MQ=60)

                                                                                        GCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCA     <  SRR4302540.8004689/94‑1 (MQ=60)

                                                                                         CGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAT    >  SRR4302540.2073430/1‑94 (MQ=60)

                                                                                          GCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAAT   <  SRR4302540.6360469/94‑1 (MQ=60)

                                                                                          GCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACCCGATGTTTTCCTGCCACCGATGTTTGCCATATCGTCCAAT   >  SRR4302540.2649481/1‑94 (MQ=60)

                                                                                           CACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATT  <  SRR4302540.2530486/94‑1 (MQ=60)

AAGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATT  >  NC_000913/1533428‑1533612

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