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 I175 R1	328	64.7	2582366	96.2%	2484236	120.9

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	918,237	T→A	N14K (AAT→AAA)	yljB →	protein YljB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	918,237	0	T	A	94.7%	135.0 / 0.0	38	N14K (AAT→AAA)	yljB	protein YljB
Reads supporting (aligned to +/- strand): ref base T (1/1); new base A (20/16); total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAAT  >  NZ_CP009273/918112‑918370

gTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCg                                                                                                                          >  2:140793/1‑139 (MQ=255)

           cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  1:1012619/1‑139 (MQ=255)

           cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGAc                                                                                                                                  <  1:633306/120‑1 (MQ=255)

           cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGAc                                                                                                                                  >  2:633306/1‑120 (MQ=255)

           cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  1:1158018/1‑139 (MQ=255)

           cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  2:716463/1‑139 (MQ=255)

            gCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCaa                                                                                                                <  2:1012619/137‑1 (MQ=255)

               aaCTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATa                                                                                                           <  1:987932/139‑1 (MQ=255)

               aaCTCTAGAAATACAATATAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATa                                                                                                           <  1:910126/139‑1 (MQ=255)

                      gAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAgatga                                                                                                      <  1:1124893/137‑1 (MQ=255)

                              aTTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatcta                                                                                            >  1:28506/1‑139 (MQ=255)

                                   gagaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCAt                                                                                        <  1:99274/138‑1 (MQ=255)

                                     gaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCTCTGATGCCAACGATAGATGATAGTTATCTATCAtgtg                                                                                     <  1:191715/139‑1 (MQ=255)

                                                aGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTgg                                                                          >  2:1112301/1‑139 (MQ=255)

                                                         gCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCaaa                                                                 <  2:322738/139‑1 (MQ=255)

                                                             ttGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGa                                                                              >  1:779643/1‑122 (MQ=255)

                                                             ttGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGa                                                                              <  2:779643/122‑1 (MQ=255)

                                                                           ggggCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGc                                               <  2:29778/139‑1 (MQ=255)

                                                                                aGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCa                                              >  1:258636/1‑135 (MQ=255)

                                                                                aGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCa                                              <  2:258636/135‑1 (MQ=255)

                                                                                aGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGc                                          >  2:972441/1‑139 (MQ=255)

                                                                                 ggTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCAt                                                                                        <  2:1104713/92‑1 (MQ=255)

                                                                                 ggTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCAt                                                                                        >  1:1104713/1‑92 (MQ=255)

                                                                                   tATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGtt                                       <  1:1112301/139‑1 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGtt                                                                                                 <  1:884146/77‑1 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGtt                                                                                                 >  2:884146/1‑77 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAAt                                                                >  2:875167/1‑110 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAAt                                                                <  1:875167/110‑1 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGc                                   >  2:593282/1‑139 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGc                                   >  1:135162/1‑139 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGc                                   >  2:718013/1‑139 (MQ=255)

                                                                                       aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGc                                   >  2:898862/1‑139 (MQ=255)

                                                                                            aaaaTTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGAt                               >  1:213749/1‑138 (MQ=255)

                                                                                                          ttAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAga                <  1:716463/139‑1 (MQ=255)

                                                                                                             aCGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAgatga             <  2:135162/139‑1 (MQ=255)

                                                                                                                   ttGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACt        >  2:956053/1‑138 (MQ=255)

                                                                                                                    tGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGa      >  1:331050/1‑139 (MQ=255)

                                                                                                                        aGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACTGAAGATGATAACTGACAAt  >  1:578814/1‑139 (MQ=255)

GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAAT  >  NZ_CP009273/918112‑918370

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

CCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAATGGGTAA  >  NZ_CP009273/918092‑918376

CCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGAC                                                                                                                                        <  SRR3722129.755476/151‑1 (MQ=60)

    aaatgtttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagaGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGAT                                                                                                                                    <  SRR3722129.279340/72‑1 (MQ=60)

           aagagacagGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTG                                                                                                                             <  SRR3722129.642719/142‑1 (MQ=60)

                    GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACG                                                                                                                    >  SRR3722129.1026917/1‑151 (MQ=60)

                    GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACG                                                                                                                    >  SRR3722129.1174168/1‑151 (MQ=60)

                                  CAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTA                                                                                                      <  SRR3722129.1001895/151‑1 (MQ=60)

                                  CAACTCTAGAAATACAATATAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTA                                                                                                      <  SRR3722129.923021/151‑1 (MQ=60)

                                       CTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTAT                                                                                                 >  SRR3722129.28871/1‑151 (MQ=60)

                                         AGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTAT                                                                                                 <  SRR3722129.1140571/149‑1 (MQ=60)

                                            catacgagattcgccttagtctcgtgggctcggagatgtgtataagagacAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGT                                                                                            <  SRR3722129.896744/101‑1 (MQ=60)

                                                      AGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGAT                                                                                   <  SRR3722129.100635/150‑1 (MQ=60)

                                                        AGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCTCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGG                                                                                <  SRR3722129.194399/151‑1 (MQ=60)

                                                                      CCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCActgtctc                                                                  >  SRR3722129.790947/1‑144 (MQ=60)

                                                                             agatgtgtataagagacAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTC                                                           <  SRR3722129.887647/134‑1 (MQ=60)

                                                                                         GGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCG                                               >  SRR3722129.262284/1‑151 (MQ=60)

                                                                                          GACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATctgtctcttatacacatctccgagcccacgagacta                                               >  SRR3722129.1120156/1‑114 (MQ=60)

                                                                                                GGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCC                                        >  SRR3722129.136986/1‑151 (MQ=60)

                                                                                                     GGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATA                                    >  SRR3722129.216729/1‑150 (MQ=60)

                                                                                                      GTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAAC                                  <  SRR3722129.1127828/151‑1 (MQ=60)

                                                                                                                             GTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGAC           >  SRR3722129.335821/1‑151 (MQ=60)

                                                                                                                             GTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGAC           <  SRR3722129.726974/151‑1 (MQ=60)

                                                                                                                                 ACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACTGAAGATGATAACTGACAATG       >  SRR3722129.587360/1‑151 (MQ=60)

                                                                                                                                      ATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtc  >  SRR3722129.1191296/1‑94 (MQ=60)

CCGGGGTCAGAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAATGGGTAA  >  NZ_CP009273/918092‑918376

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