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	579,413	(T)_6→5	coding (277/750 nt)	appY →	DNA‑binding transcriptional activator AppY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	579,408	0	T	.	96.1%	207.2 / 2.8	51	coding (272/750 nt)	appY	DNA‑binding transcriptional activator AppY
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (20/29); total (21/30)
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.25e-01

CTCTAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACA  >  NZ_CP009273/579271‑579532

ctctAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGtt                                                                                                                             >  1:1023805/1‑139 (MQ=255)

ctctAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGtt                                                                                                                             >  2:1223387/1‑139 (MQ=255)

     aaCTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGtttttct                                                                                                                        <  2:926369/139‑3 (MQ=255)

     aaCTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGtttttct                                                                                                                        <  1:1223387/139‑3 (MQ=255)

          aGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGtt                                                                                                                  >  2:177923/1‑139 (MQ=255)

          aGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGtt                                                                                                                  >  2:484528/1‑139 (MQ=255)

              ttgttgATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGaa                                                                                                              <  1:514459/139‑1 (MQ=255)

              ttgttgATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGaa                                                                                                              <  1:718050/139‑1 (MQ=255)

               tgttgATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAg                                                                                                             >  1:770758/1‑139 (MQ=255)

                           tAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAg                                                                                                 >  1:678220/1‑139 (MQ=255)

                             aGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCa                                                                                               <  1:452695/139‑1 (MQ=255)

                             aGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCa                                                                                               <  2:1023805/139‑1 (MQ=255)

                                 aaCATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACt                                                                                           <  2:1119545/139‑1 (MQ=255)

                                 aaCATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACt                                                                                           >  2:572305/1‑139 (MQ=255)

                                    aTGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcac                                                                                        <  1:133304/139‑1 (MQ=255)

                                     tGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcaca                                                                                       <  2:340217/139‑1 (MQ=255)

                                     tGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcaca                                                                                       <  2:281791/139‑1 (MQ=255)

                                     tGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcaca                                                                                       <  2:20965/139‑1 (MQ=255)

                                     tGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcaca                                                                                       <  1:1087072/139‑1 (MQ=255)

                                     tGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTcaca                                                                                       <  2:841916/139‑1 (MQ=255)

                                      gTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCacac                                                                                      >  2:1230672/1‑139 (MQ=255)

                                         tATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTc                                                                                    >  2:463399/1‑139 (MQ=255)

                                                ttttATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGt                                                                            >  1:589474/1‑139 (MQ=255)

                                                   tATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTa                                                                          >  1:1246558/1‑138 (MQ=255)

                                                           tGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 <  2:591620/139‑1 (MQ=255)

                                                           tGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 <  2:770758/139‑1 (MQ=255)

                                                           tGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 <  1:177923/139‑1 (MQ=255)

                                                           tGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 <  1:905106/139‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  1:463399/139‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  1:572305/138‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  2:589474/138‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  1:657627/138‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  1:660248/138‑1 (MQ=255)

                                                             aataatGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACt                                                               <  2:678220/139‑1 (MQ=255)

                                                                  tGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGta                                                          >  2:944797/1‑139 (MQ=255)

                                                                      cTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAaaaga                                                                     >  2:1089424/1‑124 (MQ=255)

                                                                      cTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAaaaga                                                                     <  1:1089424/124‑1 (MQ=255)

                                                                      cTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATc                                                      <  2:895511/139‑1 (MQ=255)

                                                                         ccTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 <  1:1086764/125‑1 (MQ=255)

                                                                         ccTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAaagaaaga                                                                 >  2:1086764/1‑125 (MQ=255)

                                                                         ccTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCt                                                    >  2:359955/1‑138 (MQ=255)

                                                                                       tctTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTc                                      >  2:150322/1‑138 (MQ=255)

                                                                                                     gCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAgc                       >  1:321376/1‑138 (MQ=255)

                                                                                                        atatTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  1:73812/95‑1 (MQ=255)

                                                                                                        atatTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                >  2:73812/1‑95 (MQ=255)

                                                                                                               tGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCAt             >  2:765315/1‑139 (MQ=255)

                                                                                                                   gagGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTta          >  2:809164/1‑138 (MQ=255)

                                                                                                                    agGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTtata        <  1:809164/139‑1 (MQ=255)

                                                                                                                    agGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTtata        <  2:1289406/139‑1 (MQ=255)

                                                                                                                      gAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATAc       >  2:681307/1‑138 (MQ=255)

                                                                                                                        aCCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTaa    >  2:600588/1‑139 (MQ=255)

                                                                                                                          cTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACa  <  1:681307/139‑1 (MQ=255)

                                                                                                                          cTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACa  <  1:150322/139‑1 (MQ=255)

                                                                                                                           tCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                >  1:1133417/1‑76 (MQ=255)

                                                                                                                           tCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGAc                                                                <  2:1133417/76‑1 (MQ=255)

CTCTAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGACGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACA  >  NZ_CP009273/579271‑579532

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

ATGCCCCTCTAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATATAAATACA  >  NZ_CP009273/579265‑579544

ATGCCCCTCTAAACTAAGAATTGTTGATATAGATAAGAAAACACGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTC‑TGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTC                                                                                                                                  <  SRR3722129.805061/150‑1 (MQ=60)

          AAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGG                                                                                                                        <  SRR3722129.1240414/151‑1 (MQ=60)

          AAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGG                                                                                                                        >  SRR3722129.781951/1‑151 (MQ=60)

                   ATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAAT                                                                                                               <  SRR3722129.728581/151‑1 (MQ=60)

                   ATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAAT                                                                                                               <  SRR3722129.521978/151‑1 (MQ=60)

                      GTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGC                                                                                                            >  SRR3722129.688245/1‑151 (MQ=60)

                                  AAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTC                                                                                                <  SRR3722129.459244/151‑1 (MQ=60)

                                         CATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGG                                                                                         <  SRR3722129.135101/151‑1 (MQ=60)

                                          ATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGT                                                                                        <  SRR3722129.1102287/151‑1 (MQ=60)

                                           TGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTT                                                                                       >  SRR3722129.598193/1‑151 (MQ=60)

                                              TTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAA                                                                                     >  SRR3722129.1263931/1‑150 (MQ=60)

                                                             GACAGTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTA                                                                      <  SRR3722129.1104668/151‑1 (MQ=60)

                                                                 acagATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATC                                                                  <  SRR3722129.1101973/147‑1 (MQ=60)

                                                                 GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATC                                                                  <  SRR3722129.917933/151‑1 (MQ=60)

                                                                 GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATC                                                                  <  SRR3722129.180364/151‑1 (MQ=60)

                                                                  gtctcgtgggctcggagatgtgtataagagacagTAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCC                                                                 <  SRR3722129.74773/117‑1 (MQ=60)

                                                                   GAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCC                                                                 <  SRR3722129.470138/151‑1 (MQ=60)

                                                                   GAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCC                                                                 <  SRR3722129.580773/150‑1 (MQ=60)

                                                                   GAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCC                                                                 <  SRR3722129.667404/150‑1 (MQ=60)

                                                                   GAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCC                                                                 <  SRR3722129.670050/150‑1 (MQ=60)

                                                                                                 TGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGCG                                  >  SRR3722129.326008/1‑151 (MQ=60)

                                                                                                                       GAAGAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgc            >  SRR3722129.1149199/1‑98 (MQ=60)

                                                                                                                          GAGGAACCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATAT         <  SRR3722129.820808/151‑1 (MQ=60)

                                                                                                                                CCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATATAAATAC   <  SRR3722129.691371/151‑1 (MQ=60)

                                                                                                                                CCTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATATAAATAC   <  SRR3722129.152417/151‑1 (MQ=60)

                                                                                                                                 gtgtataagagacaggTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATATAAATA    <  SRR3722129.32943/134‑1 (MQ=60)

                                                                                                                                 CTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACctgtctcttatacacat  >  SRR3722129.640461/1‑134 (MQ=60)

                                                                                                                                 CTCCATTATCGCTTG‑TTTTTCTCTGTTTGAAGGGATTAAAAT                                                                                                               >  SRR3722129.576132/1‑42 (MQ=60)

ATGCCCCTCTAAACTAAGAATTGTTGATATAGATAAGAAAACATGTTTATCCTTTTTTATCGAC‑GTGAATAATGAGCTGCCTGGCAAATTTACTCTTGATAAGAATGGCTATATTGCTGAAGAGGAACCTCCATTATCGCTTGTTTTTTCTCTGTTTGAAGGGATTAAAATAGCAGACTCACACTCCCTTTGGTTAAAAGAAAGACTATGTATATCCTTACTTGCCATGTTCAAAAAACGCGAAAGTGTAAATTCATTTATACTAACAAATATAAATACA  >  NZ_CP009273/579265‑579544

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: `‑`