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 I176 R1	180	42.9	1710034	96.5%	1650182	119.3

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	443,081	C→T	L27L (CTG→CTA)	cyoE ←	heme o synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	443,081	0	C	T	100.0%	133.5 / NA	36	L27L (CTG→CTA)	cyoE	heme o synthase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (18/18); total (18/18)

CGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGT  >  NZ_CP009273/442949‑443201

cGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGa                                                                                                                       <  1:628192/137‑1 (MQ=255)

 gTCCTATCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAAt                                                                                                                     <  1:713557/138‑1 (MQ=255)

       tCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCaa                                                                                                              >  1:273000/1‑139 (MQ=255)

        ccATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCaa                                                                                                              <  1:432047/138‑1 (MQ=255)

                    cGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAg                                                                                                 <  2:263249/139‑1 (MQ=255)

                     gATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAgg                                                                                                <  2:693900/139‑1 (MQ=255)

                             tGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCa                                                                                          >  1:73828/1‑137 (MQ=255)

                                       gttgttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGcc                                                                               <  2:608207/138‑1 (MQ=255)

                                         tgttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTgg                                                                            >  2:735305/1‑139 (MQ=255)

                                           ttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                                                          >  2:573660/1‑139 (MQ=255)

                                           ttAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGtt                                                                          <  1:735305/139‑1 (MQ=255)

                                             aaaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                                                        >  1:86092/1‑139 (MQ=255)

                                             aaaCACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTc                                                                        >  2:661636/1‑139 (MQ=255)

                                                      aCCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑c                                                              <  2:82218/139‑1 (MQ=255)

                                                          gACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAgg                                                           >  2:633142/1‑138 (MQ=255)

                                                                 caACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCtt                                                   <  1:573660/139‑1 (MQ=255)

                                                                  aaCCATTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTa                                                  <  2:73828/139‑1 (MQ=255)

                                                                         gACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGcc                                                                                           >  2:23551/1‑92 (MQ=255)

                                                                         gACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGcc                                                                                           <  1:23551/92‑1 (MQ=255)

                                                                                cAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAAccgc                                    >  1:550184/1‑139 (MQ=255)

                                                                                  aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAg                                                                                                 <  2:598762/77‑1 (MQ=255)

                                                                                  aCCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAg                                                                                                 >  1:598762/1‑77 (MQ=255)

                                                                                      gCGTGTAGATAAACAGGGTATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTCTTGCTTAAACATCATAACCGCCGctct                              >  2:282259/1‑139 (MQ=255)

                                                                                      gCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGctct                              >  1:808646/1‑139 (MQ=255)

                                                                                           aaaaaaaaCGGGGAAAAATCATGGCTTCCCTTTG‑GGCCCTTTAGAACCCCCCAAACCCCCGAGTTGAGTTTCCAAAAAATGATCCCTGTTTTCTTTCCTTG‑CGGGTTTTGTTAAAACATCAAAACCGCCGCTCTTAGt                          <  2:332539/133‑1 (MQ=255)

                                                                                             gATAAACGGGGTATAATCAATGCTCCCCTTTG‑AGGCCAGTAGGAACCCCCCAATCCCCGAGATCAGGTTGCAAAAGATGATGCCTGTTTTCTTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGca                       <  1:633142/139‑1 (MQ=255)

                                                                                             cagaaaCAGGGAATAATAAATGCTGCCTTTGG‑AGGCCAGTAGGAACCCCCCAATCACCGAGATCAGGTTGCAAAAGATGATGCTTGTTTTCTTTACTTG‑CAGGTTTTGCTTAAACATCAAAACCGCCGCTCTTAGTGca                       <  1:87399/136‑1 (MQ=255)

                                                                                                aaaCGGGGGAAAATAAGTGCT‑CCCTTTGTAGGCCATTAGAATCCCCCCAAGCCCGGAGATCA‑GTTGCCAAAGACGCTTCCCGGTTTTGTTTTTTGTCGGGTTTTTCTAAAACTGCATAACCGCCGCTCTTAGTGcatca                    <  1:78745/139‑1 (MQ=255)

                                                                                                aaaCGGGGAAAAAACAATGCTGCCCTTTG‑AGGCAAGTAGGAACCCCCCAACCCCCGAGATCAGGTTGCAAAAGATGATGCCTGGTTCCTTTACTTG‑CAGTTATTGTTAAAACATCATAACCGCCGCTCTTAGTGcatca                    <  1:663375/139‑1 (MQ=255)

                                                                                                  aCAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGcatcat                   >  1:301044/1‑138 (MQ=255)

                                                                                                     ggggATAATCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtg               >  2:6941/1‑139 (MQ=255)

                                                                                                           aaTCAATGCTGCCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtgttgctg         >  2:87399/1‑136 (MQ=255)

                                                                                                                 tgctgcCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAgg   >  1:809912/1‑139 (MQ=255)

                                                                                                                  gctgcCCTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGt  >  1:5748/1‑139 (MQ=255)

                                                                                                                      cccTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTa                                                         >  1:16060/1‑80 (MQ=255)

                                                                                                                      cccTTTG‑AGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTa                                                         <  2:16060/80‑1 (MQ=255)

CGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTG‑AGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTG‑CAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGT  >  NZ_CP009273/442949‑443201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 24 ≤` `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

CAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACA  >  NZ_CP009273/442934‑443227

aagagacaGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAG                                                                                                                                                  <  SRR3722130.843646/142‑1 (MQ=60)

           TCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAAT                                                                                                                                      >  SRR3722130.277235/1‑151 (MQ=60)

              TCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATC                                                                                                                                     <  SRR3722130.637601/149‑1 (MQ=60)

               CGTCCTATCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCAC                                                                                                                                   <  SRR3722130.724101/150‑1 (MQ=60)

                      TCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATC                                                                                                                            <  SRR3722130.438856/150‑1 (MQ=60)

                                 GTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAA                                                                                                                  >  SRR3722130.74881/1‑149 (MQ=60)

                                         tagtctcgtgggctcggagatgtgtataagagacaGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGC                                                                                                         <  SRR3722130.23862/115‑1 (MQ=60)

                                                 ATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCG                                                                                                >  SRR3722130.87303/1‑151 (MQ=60)

                                                         GTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGC                                                                                        <  SRR3722130.746100/151‑1 (MQ=60)

                                                                               ACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAA                                                                  <  SRR3722130.582325/151‑1 (MQ=60)

                                                                                    CAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCC                                                             >  SRR3722130.558541/1‑151 (MQ=60)

                                                                                      GTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatg                                                           >  SRR3722130.607771/1‑99 (MQ=60)

                                                                                          CACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTT                                                       >  SRR3722130.820442/1‑151 (MQ=60)

                                                                                                      CGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATC                                            >  SRR3722130.305711/1‑150 (MQ=60)

                                                                                                           acagAAACAGGGAATAATAAATGCTGCCTTTGGAGGCCAGTAGGAACCCCCCAATCACCGAGATCAGGTTGCAAAAGATGATGCTTGTTTTCTTTACTTGCAGGTTTTGCTTAAACATCAAAACCGCCGCTCTTAGTGCATCATCATGTTG                                      <  SRR3722130.88630/147‑1 (MQ=60)

                                                                                                           AGATAAACGGGGTATAATCAATGCTCCCCTTTGAGGCCAGTAGGAACCCCCCAATCCCCGAGATCAGGTTGCAAAAGATGATGCCTGTTTTCTTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTG                                      <  SRR3722130.642618/151‑1 (MQ=60)

                                                                                                              aaaacggggaaaaaaCAATGCTGCCCTTTGAGGCAAGTAGGAACCCCCCAACCCCCGAGATCAGGTTGCAAAAGATGATGCCTGGTTCCTTTACTTGCAGTTATTGTTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAG                                   <  SRR3722130.673266/136‑1 (MQ=60)

                                                                                                                     GGGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGT                            >  SRR3722130.821717/1‑151 (MQ=60)

                                                                                                                      GGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCActgtctcttatacacatctccg                           >  SRR3722130.623796/1‑129 (MQ=60)

                                                                                                                      GGATAATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTT                           >  SRR3722130.5828/1‑151 (MQ=60)

                                                                                                                          AATCAATGCTGCCCTTTGAGGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgt                       >  SRR3722130.16274/1‑102 (MQ=60)

                                                                                                                            atgtgtataagagacagGGCCAGTAGGAACCCCCCAATCACCGAGATCAGTTTGCAAAAGATGAGGCCTGGTTTCTTTACTTGCGGGTTTTGCTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACAT                     <  SRR3722130.735192/134‑1 (MQ=60)

                                                                                                                             CAATGCTGCCCTTTGAGCCCAGTAGAATCCCCCCAACCACCGAGTTCAGGTTGCCAAATATGATCCTTGTTTTCTTTCTTTTCAGTTTTTGTTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACAT                     <  SRR3722130.427199/150‑1 (MQ=47)

                                                                                                                                             GGCCAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaagacaaataaaaatgaaatctaatgtat    >  SRR3722130.636959/1‑45 (MQ=60)

                                                                                                                                                CAGTAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACA  >  SRR3722130.458733/1‑150 (MQ=60)

CAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACA  >  NZ_CP009273/442934‑443227

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