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	753,637	A→G	N165D (AAT→GAT)	sdhB →	succinate dehydrogenase iron‑sulfur subunit SdhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	753,637	0	A	G	89.4%	125.3 / 7.0	47	N165D (AAT→GAT)	sdhB	succinate dehydrogenase iron‑sulfur subunit SdhB
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (28/14); total (30/17)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

TTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGC  >  NZ_CP009273/753508‑753768

ttGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGAt                                                                                                                             >  1:761526/1‑138 (MQ=255)

 tGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGAt                                                                                                                             >  1:621656/1‑137 (MQ=255)

               aaTCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCgg                                                                                                             <  1:1211501/139‑1 (MQ=255)

                 tCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTATTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCa                                                                                                           <  1:537057/139‑1 (MQ=255)

                                         gATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTc                                                                                   <  2:621656/139‑1 (MQ=255)

                                          aTGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTcc                                                                                  >  1:1261164/1‑139 (MQ=255)

                                            gCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTTCt                                                                                 >  2:1047927/1‑138 (MQ=255)

                                                    agcGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAg                                                                        <  2:1088907/139‑1 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCatat                                                                                         >  2:513582/1‑120 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCatat                                                                                         <  1:513582/120‑1 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGGTTGATAGcc                                                                      >  2:531771/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTTGTAGcc                                                                      >  2:26368/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTTATAGcc                                                                      >  1:1165899/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGcc                                                                      >  1:94488/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGcc                                                                      >  1:1134121/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGcc                                                                      >  1:758724/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGcc                                                                      >  1:685970/1‑139 (MQ=255)

                                                      cgcgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGcc                                                                      >  1:109147/1‑139 (MQ=255)

                                                        cgAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGt                                                                    <  1:902815/139‑1 (MQ=255)

                                                                          gTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACaac                                                  >  2:725216/1‑137 (MQ=255)

                                                                                 tGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCg                                           <  1:1005604/139‑1 (MQ=255)

                                                                                 tGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTTGCGGGATATTGTTTCCTGATTGATAGCCGTGATACCGGGGCTGACAGCCCCCCCg                                           >  2:469471/1‑139 (MQ=255)

                                                                                   tATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGGTTGTTTGCGGCCTTTTGTTTTCTGGTTTATAGCCGTGGTACCGGGACTGACAGCCGCCTCg                                           >  2:1141059/1‑137 (MQ=255)

                                                                                     ttCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGAc                                         >  2:1235895/1‑137 (MQ=255)

                                                                                     ttCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGAc                                         <  1:1235895/137‑1 (MQ=255)

                                                                                     ttCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACgg                                       >  2:123382/1‑139 (MQ=255)

                                                                                     ttCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACgg                                       >  1:286709/1‑139 (MQ=255)

                                                                                       ctctGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGt                                      <  2:109147/138‑1 (MQ=255)

                                                                                          tGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGttt                                    >  2:420135/1‑137 (MQ=255)

                                                                                          tGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGttt                                    >  1:944538/1‑137 (MQ=255)

                                                                                           gcgcATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAg                                 >  1:790291/1‑139 (MQ=255)

                                                                                               aTGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGGGACTGACAGCCGCCTCGACGGGTTGAGTGAt                             >  1:512865/1‑139 (MQ=255)

                                                                                                    gttCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCAt                         <  2:790291/138‑1 (MQ=255)

                                                                                                       cAACCTCTTGTCCGTCTTTCTGGTGGAAGCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCa                      <  1:95257/138‑1 (MQ=255)

                                                                                                          ccTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGGAGGGCTGTTCGCGGGCTCTCGGTTTCTTTTTGAGAGGCGCGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGAGGCATTTAGCg                   >  2:547588/1‑138 (MQ=255)

                                                                                                             ctTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTAtt               >  1:829564/1‑139 (MQ=255)

                                                                                                             ctTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGGTTGAGGGATGCATTCAGCGTAtt               >  1:407296/1‑139 (MQ=255)

                                                                                                              ttGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTAt                <  1:911867/137‑1 (MQ=255)

                                                                                                               tGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTTTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTcc             >  1:634188/1‑139 (MQ=255)

                                                                                                               tGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTcc             >  2:315144/1‑139 (MQ=255)

                                                                                                               tGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGGGATGCATTCAGCGTATTcc             >  1:1149814/1‑139 (MQ=255)

                                                                                                                    gTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGt        <  1:623925/139‑1 (MQ=255)

                                                                                                                    gTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGt        <  1:1265559/139‑1 (MQ=255)

                                                                                                                        ttCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTca      >  2:952170/1‑137 (MQ=255)

                                                                                                                         tCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTcaca    <  1:952170/138‑1 (MQ=255)

                                                                                                                         tCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTcaca    <  2:94488/138‑1 (MQ=255)

                                                                                                                             gtggGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGc  <  2:1261164/136‑1 (MQ=255)

TTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGC  >  NZ_CP009273/753508‑753768

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

AAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACT  >  NZ_CP009273/753491‑753778

AAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCC                                                                                                                                           <  SRR3722129.1101490/151‑1 (MQ=60)

      AGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATA                                                                                                                                      >  SRR3722129.772617/1‑150 (MQ=60)

       GCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATA                                                                                                                                      >  SRR3722129.630893/1‑149 (MQ=60)

                               AAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTA                                                                                                            <  SRR3722129.1228389/151‑1 (MQ=60)

                                 ATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTATTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGC                                                                                                          <  SRR3722129.544907/151‑1 (MQ=60)

                                                AGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGGATATCGTTTCCT                                                                                           >  SRR3722129.1278701/1‑151 (MQ=60)

                                                   aagagacagTGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGAT                                                                                        <  SRR3722129.521091/142‑1 (MQ=60)

                                                    agagacagTGCCAGAGCAGCGCGAACAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTGTCTGGTGGGATCCCGATAAGTATATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATT                                                                                       <  SRR3722129.1019626/143‑1 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgc                                                                               >  SRR3722129.788441/1‑98 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTTATAGCCG                                                                               >  SRR3722129.1182152/1‑151 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCG                                                                               >  SRR3722129.769777/1‑151 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCG                                                                               >  SRR3722129.1149912/1‑151 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCG                                                                               >  SRR3722129.696083/1‑151 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGGTTCCTGATTGATAGCCG                                                                               >  SRR3722129.95767/1‑151 (MQ=60)

                                                            TGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCG                                                                               >  SRR3722129.110631/1‑151 (MQ=60)

                                                                       gagatgtgtataagagacagGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGA                                                                    <  SRR3722129.438553/131‑1 (MQ=60)

                                                                        GCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGAC                                                                   <  SRR3722129.915619/151‑1 (MQ=60)

                                                                                           GTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGG                                                >  SRR3722129.290820/1‑151 (MQ=60)

                                                                                                AATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTG                                             >  SRR3722129.957899/1‑149 (MQ=60)

                                                                                                 ATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGT                                          <  SRR3722129.1019770/151‑1 (MQ=60)

                                                                                                 ATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGG                                          >  SRR3722129.801747/1‑151 (MQ=60)

                                                                                                   GTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGA                                        <  SRR3722129.1253060/151‑1 (MQ=60)

                                                                                                     ATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGGGACTGACAGCCGCCTCGACGGGTTGAGTGATG                                      >  SRR3722129.520369/1‑151 (MQ=60)

                                                                                                                   TTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCtgtctcttatacacatctccgagcccacgagactaaggcgaatatcgtatgccgtcttctgcttgaaaaaaaaaaagatagcgcgcagtcttacccggtactgaaactgcc                           >  SRR3722129.62626/1‑37 (MQ=60)

                                                                                                                   TTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTC                        >  SRR3722129.841488/1‑151 (MQ=60)

                                                                                                                   TTGTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGGTTGAGGGATGCATTCAGCGTATTC                        >  SRR3722129.413203/1‑151 (MQ=60)

                                                                                                                     GTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTTTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCG                      >  SRR3722129.643619/1‑151 (MQ=60)

                                                                                                                     GTTCAACCTCTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGGTTGAGGGATGCATTCAGCGGATTCCG                      >  SRR3722129.1165790/1‑151 (MQ=60)

                                                                                                                       TCAACCTCTTGTCCGTCTTTCTGGTGGAAGCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGC                     <  SRR3722129.96546/150‑1 (MQ=60)

                                                                                                                              CTTGTCCGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCAC               <  SRR3722129.924798/149‑1 (MQ=60)

                                                                                                                                    CGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCAT       <  SRR3722129.633198/151‑1 (MQ=60)

                                                                                                                                    CGTCTTTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCAT       <  SRR3722129.1283149/151‑1 (MQ=60)

                                                                                                                                         TTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAAC   <  SRR3722129.965625/150‑1 (MQ=60)

                                                                                                                                         TTCTGGTGGGATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGGTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACC  >  SRR3722129.497588/1‑151 (MQ=60)

AAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACT  >  NZ_CP009273/753491‑753778

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