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
A2 F5 I184 R1	278	29.6	1166432	97.3%	1134938	124.1

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,568,721	A→G	I145T (ATC→ACC)	mcrB ←	5‑methylcytosine‑specific restriction endonuclease subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,568,721	0	A	G	100.0%	139.4 / NA	37	I145T (ATC→ACC)	mcrB	5‑methylcytosine‑specific restriction endonuclease subunit McrB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (19/18); total (19/18)

CGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTT  >  NZ_CP009273/4568597‑4568856

cGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtg                                                                                                                                    >  2:405427/1‑130 (MQ=255)

cGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtg                                                                                                                                    <  1:405427/130‑1 (MQ=255)

             gTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATc                                                                                                              >  2:171348/1‑139 (MQ=255)

                tCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCggg                                                                                                           >  2:35316/1‑139 (MQ=255)

                tCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCggg                                                                                                           >  1:497304/1‑139 (MQ=255)

                tCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCgg                                                                                                            <  2:197740/138‑1 (MQ=255)

                  aGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgt                                                                                                          <  2:396537/138‑1 (MQ=255)

                  aGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgt                                                                                                          <  1:452550/138‑1 (MQ=255)

                   gggATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                                                                        <  2:562643/139‑1 (MQ=255)

                   gggATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgtgt                                                                                                        <  2:242637/139‑1 (MQ=255)

                                                cAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGaa                                                                           <  1:258371/139‑1 (MQ=255)

                                                    cAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                                                                       <  2:458148/139‑1 (MQ=255)

                                                    cAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                                                                       <  2:200378/139‑1 (MQ=255)

                                                     aGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGAc                                                                                     >  1:228499/1‑124 (MQ=255)

                                                     aGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGAc                                                                                     <  2:228499/124‑1 (MQ=255)

                                                      gTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  2:383452/1‑139 (MQ=255)

                                                       tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATa                                                                                                <  1:277078/111‑1 (MQ=255)

                                                       tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATa                                                                                                >  2:277078/1‑111 (MQ=255)

                                                       tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  1:248223/1‑138 (MQ=255)

                                                       tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  2:36882/1‑138 (MQ=255)

                                                       tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  2:438560/1‑138 (MQ=255)

                                                            ttCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAg                                                               <  1:171348/139‑1 (MQ=255)

                                                               aGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGAc                                                            >  2:393937/1‑139 (MQ=255)

                                                                 tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAAt                                                                                                     >  1:456440/1‑96 (MQ=255)

                                                                 tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAAt                                                                                                     <  2:456440/96‑1 (MQ=255)

                                                                 tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTg                                                          >  2:485401/1‑139 (MQ=255)

                                                                           aGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAAt                                                                                                               >  2:158055/1‑76 (MQ=255)

                                                                           aGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAAt                                                                                                               <  1:158055/76‑1 (MQ=255)

                                                                                  aTAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGa                                          >  1:199925/1‑138 (MQ=255)

                                                                                           ttGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACAccc                                 >  2:197822/1‑138 (MQ=255)

                                                                                                 gAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTg                          <  2:140134/139‑1 (MQ=255)

                                                                                                                   taGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGa        <  1:197822/139‑1 (MQ=255)

                                                                                                                     gTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGAt       <  2:199925/138‑1 (MQ=255)

                                                                                                                        gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGAt       <  1:558442/135‑1 (MQ=255)

                                                                                                                        gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGAt       >  2:558442/1‑135 (MQ=255)

                                                                                                                         ttGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGttt  >  1:340326/1‑139 (MQ=255)

                                                                                                                         ttGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGttt  >  2:500770/1‑139 (MQ=255)

CGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTT  >  NZ_CP009273/4568597‑4568856

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

TTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  >  NZ_CP009273/4568573‑4568871

TTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGT                                                                                                                                                       <  SRR3722139.279180/150‑1 (MQ=60)

TTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGT                                                                                                                                                       <  SRR3722139.4239/150‑1 (MQ=60)

       GTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGT                                                                                                                                               <  SRR3722139.9115/151‑1 (MQ=60)

              GTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCAT                                                                                                                                        <  SRR3722139.410220/151‑1 (MQ=60)

                             CTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGT                                                                                                                         >  SRR3722139.503055/1‑151 (MQ=60)

                                   gcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATC                                                                                                                   <  SRR3722139.159985/98‑1 (MQ=60)

                                         CAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATAC                                                                                                              <  SRR3722139.457843/150‑1 (MQ=60)

                                                  gatgtgtataagagacagCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGAC                                                                                                    <  SRR3722139.280574/133‑1 (MQ=60)

                                                                  ATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACctgtc                                                                                    >  SRR3722139.231327/1‑146 (MQ=60)

                                                                    CTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGT                                                                                   >  SRR3722139.251368/1‑150 (MQ=60)

                                                                       CCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATA                                                                               <  SRR3722139.261629/151‑1 (MQ=60)

                                                                              GTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTctgtctcttatacacatctccgagcccacg                                                                           >  SRR3722139.461795/1‑118 (MQ=60)

                                                                                   ATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTA                                                                   <  SRR3722139.173489/151‑1 (MQ=60)

                                                                                               ACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGAT                                                        >  SRR3722139.202390/1‑150 (MQ=60)

                                                                                                                                      ATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTT                >  SRR3722139.344463/1‑151 (MQ=60)

                                                                                                                                          ATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGT            <  SRR3722139.200260/151‑1 (MQ=60)

                                                                                                                                            AGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTA           <  SRR3722139.564835/150‑1 (MQ=60)

                                                                                                                                             GTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgt         >  SRR3722139.375501/1‑95 (MQ=60)

                                                                                                                                                  TGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTG     >  SRR3722139.291403/1‑150 (MQ=60)

                                                                                                                                                   GGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGA    >  SRR3722139.381564/1‑150 (MQ=60)

                                                                                                                                                    GTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  <  SRR3722139.35698/151‑1 (MQ=60)

                                                                                                                                                    GTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  <  SRR3722139.218340/151‑1 (MQ=60)

                                                                                                                                                    GTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  <  SRR3722139.5877/151‑1 (MQ=60)

                                                                                                                                                    GTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  <  SRR3722139.162762/151‑1 (MQ=60)

TTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAG  >  NZ_CP009273/4568573‑4568871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 0 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`