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 I2 R1	884	257.2	12958910	96.9%	12557183	93.0

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,886,608	G→A	A84T (GCC→ACC)	b3708 →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,886,608	0	G	A	99.6%	901.2 / ‑5.6	235	A84T (GCC→ACC)	b3708	tryptophanase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (90/144); total (90/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 243 total aligned reads displayed.

AGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACT  >  NC_000913/3886518‑3886696

aGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAAcc                                                                                        <  1:875676/93‑1 (MQ=255)

 gCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCg                                                                                       >  2:1440923/1‑93 (MQ=255)

    gCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAgg                                                                                    >  1:522860/1‑92 (MQ=255)

    gCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGt                                                                                    <  1:3422421/93‑1 (MQ=255)

    gCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGt                                                                                    <  1:5331077/93‑1 (MQ=255)

    gCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGt                                                                                    >  2:5787882/1‑93 (MQ=255)

     cACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTc                                                                                   >  2:4266311/1‑93 (MQ=255)

     cACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTc                                                                                   <  2:6379146/93‑1 (MQ=255)

      aCCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCa                                                                                  >  1:2684816/1‑93 (MQ=255)

        cGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGt                                                                                <  1:3923313/93‑1 (MQ=255)

        cGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGt                                                                                <  1:5788148/93‑1 (MQ=255)

        cGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGt                                                                                <  2:1563717/93‑1 (MQ=255)

         ggggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTg                                                                               <  2:4628655/93‑1 (MQ=255)

          gggCGATGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGa                                                                              >  2:232997/1‑93 (MQ=255)

              ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGaaaaa                                                                          <  2:2714417/93‑1 (MQ=255)

              ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGaaaaa                                                                          <  1:1355241/93‑1 (MQ=255)

               gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAat                                                                         <  2:2600937/93‑1 (MQ=255)

               gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAat                                                                         >  1:6134144/1‑93 (MQ=255)

               gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAat                                                                         >  1:3794382/1‑93 (MQ=255)

               gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAat                                                                         >  1:5061978/1‑93 (MQ=255)

                tGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAata                                                                        <  2:346664/93‑1 (MQ=255)

                tGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAata                                                                        <  1:3042965/93‑1 (MQ=255)

                tGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAata                                                                        <  2:2010813/93‑1 (MQ=255)

                   cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCt                                                                     <  1:4530459/93‑1 (MQ=255)

                   cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCt                                                                     >  2:2206511/1‑93 (MQ=255)

                    gCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCtt                                                                    <  1:2872955/93‑1 (MQ=255)

                     cAGAGAATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCttt                                                                   <  2:6476583/93‑1 (MQ=255)

                       gagCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTgg                                                                 <  2:3883323/93‑1 (MQ=255)

                         gCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGtt                                                               <  1:2372969/93‑1 (MQ=255)

                          cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTa                                                              <  1:5174560/93‑1 (MQ=255)

                          cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTa                                                              <  2:4953642/93‑1 (MQ=255)

                           aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTAt                                                             <  1:2616635/93‑1 (MQ=255)

                               aGGCTACGATGATGCGCGGCAACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAAt                                                         <  1:1239201/93‑1 (MQ=255)

                                ggCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATa                                                        >  2:2560504/1‑93 (MQ=255)

                                ggCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATa                                                        <  2:5648622/93‑1 (MQ=255)

                                 gCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATac                                                       <  2:403667/93‑1 (MQ=255)

                                 gCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATac                                                       >  1:6233861/1‑93 (MQ=255)

                                  cTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGATAAATATCTTTGGTTATCAATaca                                                      <  2:4376310/93‑1 (MQ=255)

                                    gCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACAcc                                                    <  1:5291973/93‑1 (MQ=255)

                                     cGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCa                                                   <  1:1013974/93‑1 (MQ=255)

                                      gatgatGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAt                                                  <  1:416231/93‑1 (MQ=255)

                                      gatgatGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAt                                                  >  1:1607618/1‑93 (MQ=255)

                                       atgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAtt                                                 <  2:4869545/93‑1 (MQ=255)

                                       atgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAtt                                                 <  1:4857385/93‑1 (MQ=255)

                                        tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                <  1:2891511/93‑1 (MQ=255)

                                        tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                <  2:165476/93‑1 (MQ=255)

                                        tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                <  1:3627394/93‑1 (MQ=255)

                                        tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                <  1:2407866/93‑1 (MQ=255)

                                        tgatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTc                                                >  2:3082862/1‑93 (MQ=255)

                                         gatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAttac                                               <  1:425033/93‑3 (MQ=255)

                                         gatgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTcc                                               <  1:5352278/93‑1 (MQ=255)

                                          atgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCg                                              <  2:6031645/93‑1 (MQ=255)

                                            gcgcGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGAc                                            <  2:1333817/93‑1 (MQ=255)

                                            gcgcGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGAc                                            <  2:3617302/93‑1 (MQ=255)

                                             cgcgGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACt                                           >  2:1253201/1‑93 (MQ=255)

                                              gcggcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTc                                          >  1:333932/1‑93 (MQ=255)

                                                 gcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCAcc                                       <  1:1550003/93‑1 (MQ=255)

                                                 gcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCAcc                                       >  2:196753/1‑93 (MQ=255)

                                                  cgacgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCa                                      <  2:735583/93‑1 (MQ=255)

                                                    acgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                    <  1:4427074/93‑1 (MQ=255)

                                                    acgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                    >  2:6328199/1‑93 (MQ=255)

                                                    acgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                    >  2:4729772/1‑93 (MQ=255)

                                                    acgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                    <  2:2808546/93‑1 (MQ=255)

                                                     cgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAggg                                   >  1:5661432/1‑93 (MQ=255)

                                                     cgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAggg                                   <  1:232667/93‑1 (MQ=255)

                                                      gaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGc                                  >  2:4660228/1‑93 (MQ=255)

                                                      gaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGc                                  >  2:5390660/1‑93 (MQ=255)

                                                       aaGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGcc                                 <  1:1656158/93‑1 (MQ=255)

                                                       aaGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGcc                                 >  2:3160759/1‑93 (MQ=255)

                                                       aaGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGcc                                 <  2:5346015/93‑1 (MQ=255)

                                                        aGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCg                                >  2:2156144/1‑93 (MQ=255)

                                                        aGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCg                                <  2:5906736/93‑1 (MQ=255)

                                                         gCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGt                               <  2:3651585/93‑1 (MQ=255)

                                                         gCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGt                               <  2:4288231/93‑1 (MQ=255)

                                                            tACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGgc                            <  2:3414042/93‑1 (MQ=255)

                                                            tACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGgc                            <  1:2190180/93‑1 (MQ=255)

                                                            tACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGgc                            <  1:3949098/93‑1 (MQ=255)

                                                              cAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGgcgc                          >  2:5166639/1‑93 (MQ=255)

                                                              cAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGgcgc                          <  1:2029220/93‑1 (MQ=255)

                                                                gCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCACGGCCGTGGCGCag                        >  1:6362178/1‑93 (MQ=255)

                                                                 cGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGTTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCaga                       <  1:1043976/93‑1 (MQ=255)

                                                                 cGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAACATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCaga                       <  2:3504153/93‑1 (MQ=255)

                                                                    cAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCa                    <  1:3166124/93‑1 (MQ=255)

                                                                      gTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                  <  2:5125765/93‑1 (MQ=255)

                                                                      gTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                  <  2:5576208/93‑1 (MQ=255)

                                                                          tAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCta              <  1:4556100/93‑1 (MQ=255)

                                                                           aGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGCCCGTGGCGCAGAGCAAATCtat             >  2:1607229/1‑93 (MQ=255)

                                                                            gCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCtata            >  2:443277/1‑93 (MQ=255)

                                                                              tactaTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATAtt          <  1:4687849/93‑1 (MQ=255)

                                                                               actaTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTc         >  2:2054066/1‑93 (MQ=255)

                                                                               actaTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTc         <  2:2823685/93‑1 (MQ=255)

                                                                                ctaTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTcc        >  1:4945517/1‑93 (MQ=255)

                                                                                ctaTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTcc        >  1:3764891/1‑93 (MQ=255)

                                                                                  aTGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCgg      <  1:678802/93‑1 (MQ=255)

                                                                                   tGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt     <  1:6452296/93‑1 (MQ=255)

                                                                                   tGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt     <  1:2979240/93‑1 (MQ=255)

                                                                                      gTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACt  >  1:800464/1‑93 (MQ=255)

                                                                                      gTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACt  <  2:1786589/93‑1 (MQ=255)

                                                                                      gTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACt  <  2:2240762/93‑1 (MQ=255)

                                                                                      gTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACt  <  2:4203592/93‑1 (MQ=255)

AGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACT  >  NC_000913/3886518‑3886696

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  >  NC_000913/3886515‑3886701

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAA                                                                                               >  SRR4302534.2857287/1‑94 (MQ=60)

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAA                                                                                               <  SRR4302534.1091936/94‑1 (MQ=60)

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAA                                                                                               >  SRR4302534.2783641/1‑94 (MQ=60)

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAA                                                                                               <  SRR4302534.560936/94‑1 (MQ=60)

  CAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACC                                                                                             <  SRR4302534.3292336/94‑1 (MQ=60)

  CAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACC                                                                                             <  SRR4302534.6308931/94‑1 (MQ=60)

   AGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCG                                                                                            <  SRR4302534.897827/94‑1 (MQ=60)

     CGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAG                                                                                          >  SRR4302534.657178/1‑94 (MQ=60)

      GGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGT                                                                                         >  SRR4302534.4364474/1‑94 (MQ=60)

      GGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGG                                                                                         >  SRR4302534.535857/1‑94 (MQ=60)

       GCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTC                                                                                        <  SRR4302534.3514380/94‑1 (MQ=60)

       GCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTC                                                                                        <  SRR4302534.5474859/94‑1 (MQ=60)

       GCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTC                                                                                        <  SRR4302534.4568017/94‑1 (MQ=60)

        CACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCA                                                                                       >  SRR4302534.2813571/1‑94 (MQ=60)

        CACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCA                                                                                       >  SRR4302534.2755777/1‑94 (MQ=60)

           CGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTG                                                                                    <  SRR4302534.3376743/94‑1 (MQ=60)

           CGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTG                                                                                    <  SRR4302534.5945162/94‑1 (MQ=60)

           CGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTG                                                                                    <  SRR4302534.4028431/94‑1 (MQ=60)

           CGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTG                                                                                    <  SRR4302534.6592980/94‑1 (MQ=60)

            GGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGA                                                                                   >  SRR4302534.2396373/1‑94 (MQ=60)

            GGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGA                                                                                   >  SRR4302534.3109418/1‑94 (MQ=60)

            GGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGA                                                                                   >  SRR4302534.3670563/1‑94 (MQ=60)

            GGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGA                                                                                   >  SRR4302534.5373378/1‑94 (MQ=60)

              GGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAA                                                                                 >  SRR4302534.666566/1‑94 (MQ=60)

                CGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAA                                                                               >  SRR4302534.1692481/1‑94 (MQ=60)

                 GGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTGTGGGTTTACCGAGTCAGTGAAAAAT                                                                              >  SRR4302534.2396117/1‑94 (MQ=60)

                 GGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAAT                                                                              <  SRR4302534.1390958/94‑1 (MQ=60)

                 GGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAAT                                                                              >  SRR4302534.6301001/1‑94 (MQ=60)

                 GGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAAT                                                                              >  SRR4302534.3896259/1‑94 (MQ=60)

                 GGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAAT                                                                              >  SRR4302534.5198671/1‑94 (MQ=60)

                  GTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTAAACCGAGTCAGTGAAAAATA                                                                             <  SRR4302534.6473402/94‑1 (MQ=60)

                   TGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATAT                                                                            <  SRR4302534.3122778/94‑1 (MQ=60)

                   TGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATAT                                                                            <  SRR4302534.2948428/94‑1 (MQ=60)

                   TGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATAT                                                                            >  SRR4302534.4191111/1‑94 (MQ=60)

                      CGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTT                                                                         <  SRR4302534.5184809/94‑1 (MQ=60)

                      CGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTT                                                                         <  SRR4302534.4653224/94‑1 (MQ=60)

                       GCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTT                                                                        <  SRR4302534.2948397/94‑1 (MQ=60)

                        CAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTG                                                                       <  SRR4302534.2121074/94‑1 (MQ=60)

                        CAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTG                                                                       <  SRR4302534.2537367/94‑1 (MQ=60)

                          GAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGT                                                                     >  SRR4302534.1646219/1‑94 (MQ=60)

                            GCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTA                                                                   <  SRR4302534.2435979/94‑1 (MQ=60)

                             CATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTAT                                                                  <  SRR4302534.5314006/94‑1 (MQ=60)

                              ATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATC                                                                 <  SRR4302534.2685879/94‑1 (MQ=60)

                                 CAGGCTGCGATGATGCGCGGCGACGACGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAAT                                                              <  SRR4302534.5342812/94‑1 (MQ=60)

                                  AGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATA                                                             >  SRR4302534.4780478/1‑94 (MQ=60)

                                  AGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATA                                                             <  SRR4302534.4807015/94‑1 (MQ=60)

                                  AGGCTACGATGATGCGCGGCAACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATA                                                             <  SRR4302534.1271860/94‑1 (MQ=60)

                                   GGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATAC                                                            >  SRR4302534.6403073/1‑94 (MQ=60)

                                    GCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACA                                                           <  SRR4302534.3994559/94‑1 (MQ=60)

                                    GCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACA                                                           <  SRR4302534.1209910/94‑1 (MQ=60)

                                       GCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCA                                                        <  SRR4302534.86060/94‑1 (MQ=60)

                                       GCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCA                                                        <  SRR4302534.5434796/94‑1 (MQ=60)

                                       GCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCA                                                        <  SRR4302534.6326145/94‑1 (MQ=60)

                                        CGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAT                                                       >  SRR4302534.1649600/1‑94 (MQ=60)

                                        CGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCAT                                                       <  SRR4302534.1040069/94‑1 (MQ=60)

                                         GATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATT                                                      <  SRR4302534.426872/94‑1 (MQ=60)

                                          ATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTC                                                     <  SRR4302534.4654826/94‑1 (MQ=60)

                                          ATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTC                                                     <  SRR4302534.4989169/94‑1 (MQ=60)

                                          ATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTC                                                     <  SRR4302534.542737/94‑1 (MQ=60)

                                           TGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCC                                                    <  SRR4302534.3077864/94‑1 (MQ=60)

                                           TGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCC                                                    <  SRR4302534.2471863/94‑1 (MQ=60)

                                           TGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCC                                                    <  SRR4302534.3724903/94‑1 (MQ=60)

                                           TGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCC                                                    <  SRR4302534.2967383/94‑1 (MQ=60)

                                            TATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCG                                                   <  SRR4302534.1320706/94‑1 (MQ=60)

                                            GATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCGGTGAAAAATATCTTTGGTTATCAATACACCATTCCG                                                   <  SRR4302534.3536184/94‑1 (MQ=60)

                                            GATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCG                                                   <  SRR4302534.5496583/94‑1 (MQ=60)

                                            GATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTACG                                                   <  SRR4302534.435893/94‑1 (MQ=60)

                                                CGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTC                                               >  SRR4302534.342500/1‑94 (MQ=60)

                                                CGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTC                                               >  SRR4302534.5750656/1‑94 (MQ=60)

                                                CGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTC                                               <  SRR4302534.4509909/94‑1 (MQ=60)

                                                 GCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTTCGACTCC                                              >  SRR4302534.5954027/1‑94 (MQ=60)

                                                   GGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACC                                            <  SRR4302534.5074619/94‑1 (MQ=60)

                                                    GCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCA                                           <  SRR4302534.777578/94‑1 (MQ=60)

                                                    GCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCA                                           <  SRR4302534.6163482/94‑1 (MQ=60)

                                                    GCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCA                                           <  SRR4302534.1590616/94‑1 (MQ=60)

                                                      GACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGG                                         <  SRR4302534.1925907/94‑1 (MQ=60)

                                                      GACGAAGCCGACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCAACAGG                                         >  SRR4302534.1042700/1‑94 (MQ=60)

                                                       ACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGG                                        <  SRR4302534.4546611/94‑1 (MQ=60)

                                                       ACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGG                                        >  SRR4302534.5814710/1‑94 (MQ=60)

                                                        CGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGC                                       <  SRR4302534.238756/94‑1 (MQ=60)

                                                          AAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCG                                     <  SRR4302534.1699323/94‑1 (MQ=60)

                                                          AAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCG                                     >  SRR4302534.3838872/1‑94 (MQ=60)

                                                            GCCTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTG                                   <  SRR4302534.5782875/94‑1 (MQ=60)

                                                              CTACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGC                                 >  SRR4302534.393476/1‑94 (MQ=60)

                                                               TACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCG                                <  SRR4302534.4054818/94‑1 (MQ=60)

                                                               TACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCG                                <  SRR4302534.2248090/94‑1 (MQ=60)

                                                                ACAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGC                               <  SRR4302534.526208/94‑1 (MQ=60)

                                                                 CAGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCA                              <  SRR4302534.2082286/94‑1 (MQ=60)

                                                                  AGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAG                             <  SRR4302534.4862977/94‑1 (MQ=60)

                                                                  AGCGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCACGGCCGTGGCGCAG                             >  SRR4302534.6533978/1‑94 (MQ=60)

                                                                    CGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGTTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAG                           <  SRR4302534.1070939/94‑1 (MQ=60)

                                                                    CGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAG                           >  SRR4302534.226294/1‑94 (MQ=60)

                                                                    CGGCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAG                           <  SRR4302534.388893/94‑1 (MQ=60)

                                                                      GCAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCA                         <  SRR4302534.5440984/94‑1 (MQ=60)

                                                                       CAGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAA                        <  SRR4302534.3249618/94‑1 (MQ=60)

                                                                        AGTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAA                       >  SRR4302534.2468747/1‑94 (MQ=60)

                                                                         GTCGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAT                      >  SRR4302534.1608816/1‑94 (MQ=60)

                                                                           CGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGCGCAAATCT                    >  SRR4302534.4076054/1‑94 (MQ=60)

                                                                           CGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCT                    <  SRR4302534.1811698/94‑1 (MQ=60)

                                                                           CGTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCT                    <  SRR4302534.3199052/94‑1 (MQ=60)

                                                                            GTAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTA                   >  SRR4302534.5255861/1‑94 (MQ=60)

                                                                            GTAGCGACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTA                   <  SRR4302534.546912/94‑1 (MQ=60)

                                                                             TAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTAT                  <  SRR4302534.4679625/94‑1 (MQ=60)

                                                                             TAGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTAT                  <  SRR4302534.4047369/94‑1 (MQ=60)

                                                                              AGCTACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATA                 <  SRR4302534.2994082/94‑1 (MQ=60)

                                                                                 TACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTC              <  SRR4302534.4815094/94‑1 (MQ=60)

                                                                                 TACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTC              <  SRR4302534.5931916/94‑1 (MQ=60)

                                                                                  ACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCCATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCC             >  SRR4302534.713857/1‑94 (MQ=60)

                                                                                  ACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCC             >  SRR4302534.5079501/1‑94 (MQ=60)

                                                                                  ACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCC             <  SRR4302534.1370679/94‑1 (MQ=60)

                                                                                  ACTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCC             >  SRR4302534.3865961/1‑94 (MQ=60)

                                                                                   CTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCG            <  SRR4302534.5011550/94‑1 (MQ=60)

                                                                                   CTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCG            >  SRR4302534.3401183/1‑94 (MQ=60)

                                                                                   CTATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCG            <  SRR4302534.2752635/94‑1 (MQ=60)

                                                                                    TATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGG           >  SRR4302534.2908669/1‑94 (MQ=60)

                                                                                     ATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGT          <  SRR4302534.695378/94‑1 (MQ=60)

                                                                                     ATGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGT          >  SRR4302534.4398254/1‑94 (MQ=60)

                                                                                     ATGCGTTAACCGAGACAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGCGCCGTGGCGCAGAGCAAATCTATATTCCGGT          >  SRR4302534.3318129/1‑94 (MQ=60)

                                                                                      TGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTA         <  SRR4302534.3057209/94‑1 (MQ=60)

                                                                                      TGCGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTA         <  SRR4302534.6626162/94‑1 (MQ=60)

                                                                                        CGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACT       >  SRR4302534.1417613/1‑94 (MQ=60)

                                                                                        CGTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACT       >  SRR4302534.820569/1‑94 (MQ=60)

                                                                                         GTTAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTG      <  SRR4302534.3451962/94‑1 (MQ=60)

                                                                                           TAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAT    >  SRR4302534.3803332/1‑94 (MQ=60)

                                                                                           TAACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAT    >  SRR4302534.3554433/1‑94 (MQ=60)

                                                                                            AACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATT   <  SRR4302534.1298679/94‑1 (MQ=60)

                                                                                            AACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATT   <  SRR4302534.5052167/94‑1 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  <  SRR4302534.4609520/94‑1 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  <  SRR4302534.395947/94‑1 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  >  SRR4302534.1719401/1‑94 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  <  SRR4302534.6461818/94‑1 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  >  SRR4302534.5667835/1‑94 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  >  SRR4302534.208447/1‑94 (MQ=60)

                                                                                             ACCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCAATATTCCGGTACTGATTA  <  SRR4302534.5553000/94‑1 (MQ=60)

GACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTA  >  NC_000913/3886515‑3886701

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