BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AP012306	2,073,848	Δ2 bp	intergenic (+106/+161)	lysV → / ← xapR	tRNA‑Lys/DNA‑binding transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP012306	2,073,847	0	G	.	99.0%	451.9 / ‑3.1	97	intergenic (+105/+163)	lysV/xapR	tRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (28/68); total (29/68)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AP012306	2,073,848	0	T	.	99.0%	452.1 / ‑1.3	97	intergenic (+106/+162)	lysV/xapR	tRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand): ref base T (0/0); major base . (28/68); minor base C (1/0); total (29/68)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 112 total aligned reads displayed.

CTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073698‑2073987

||

cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              >  1:448573/1‑150 (MQ=32)

cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              >  1:962791/1‑150 (MQ=32)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             <  1:3092908/150‑2 (MQ=14)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             <  1:4121725/150‑2 (MQ=14)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             >  1:2387398/1‑149 (MQ=14)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             <  2:3682114/150‑2 (MQ=14)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             >  2:432263/1‑149 (MQ=14)

 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             >  1:3807086/1‑149 (MQ=14)

  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAGAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcg                                                                                                                                            <  2:3379073/150‑3 (MQ=11)

  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcg                                                                                                                                            <  2:3443200/150‑3 (MQ=11)

  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGCTTGAATTGAGCGCACCGAAACAACGGAGCCGCTCgcg                                                                                                                                            >  2:1585072/1‑148 (MQ=2)

  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑CGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                                                           >  2:1272323/1‑147 (MQ=0)

  tttAATCAATTGGTCGCAGGCTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGCCTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcg                                                                                                                                            >  1:3437431/1‑148 (MQ=2)

   ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                                                           >  2:611463/1‑147 (MQ=1)

   ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                                                           >  2:513390/1‑147 (MQ=1)

                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:4015876/150‑1 (MQ=2)

                                                                                                      aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     <  1:1694981/150‑1 (MQ=255)

                                                                                                      aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     >  2:3205515/1‑150 (MQ=255)

                                                                                                      aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     <  1:1075733/150‑1 (MQ=255)

                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  2:3400028/150‑1 (MQ=255)

                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  2:2175407/150‑1 (MQ=255)

                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  1:1119281/150‑1 (MQ=255)

                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  1:862784/150‑1 (MQ=255)

                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   >  1:3157798/1‑150 (MQ=255)

                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:3872826/150‑1 (MQ=255)

                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:373887/150‑1 (MQ=255)

                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  1:1575202/150‑1 (MQ=255)

                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:586845/150‑1 (MQ=255)

                                                                                                           cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                >  2:1427897/1‑150 (MQ=255)

                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  1:4183770/1‑150 (MQ=255)

                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               <  1:3217814/150‑1 (MQ=255)

                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  2:3402198/1‑150 (MQ=255)

                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCG‑‑CGCGC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  1:4045720/1‑150 (MQ=11)

                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             >  1:1051668/1‑150 (MQ=255)

                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             >  1:2860038/1‑150 (MQ=255)

                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  1:3098999/150‑1 (MQ=255)

                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  2:3019307/150‑1 (MQ=255)

                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  1:537660/150‑1 (MQ=255)

                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  1:2628/150‑1 (MQ=255)

                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  2:3161988/150‑1 (MQ=255)

                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  1:4336811/150‑1 (MQ=255)

                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  1:3723292/150‑1 (MQ=255)

                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGATCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  1:3703764/1‑150 (MQ=14)

                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  2:733697/1‑150 (MQ=255)

                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  1:891768/150‑1 (MQ=255)

                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTCTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  2:3762134/150‑1 (MQ=14)

                                                                                                                  tttGAGTTGAGCGCAGCGAACCAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAGGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  2:776823/150‑1 (MQ=14)

                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  2:2971004/150‑1 (MQ=255)

                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  1:1000432/150‑1 (MQ=255)

                                                                                                                    tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCt                       <  2:2498077/150‑1 (MQ=255)

                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  1:1833600/150‑1 (MQ=255)

                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      >  1:872423/1‑150 (MQ=255)

                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  1:189718/150‑1 (MQ=255)

                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  2:160957/150‑1 (MQ=255)

                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTGTGCTATCTGCGATTTGCGAAATTGCCTGATgcgcgt                      >  2:3521705/1‑148 (MQ=14)

                                                                                                                      agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     <  1:1651686/150‑1 (MQ=255)

                                                                                                                      agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     >  1:2344528/1‑150 (MQ=255)

                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:2612442/150‑1 (MQ=255)

                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:2799873/150‑1 (MQ=255)

                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    >  2:2182354/1‑150 (MQ=255)

                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  2:855998/150‑1 (MQ=255)

                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:760028/150‑1 (MQ=255)

                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:376377/150‑1 (MQ=255)

                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:693702/150‑1 (MQ=255)

                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  >  1:778307/1‑150 (MQ=255)

                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 >  1:1181528/1‑150 (MQ=255)

                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  1:3218002/150‑1 (MQ=255)

                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  1:4054102/150‑1 (MQ=255)

                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                >  1:3463916/1‑150 (MQ=255)

                                                                                                                            gcgcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGTACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:3809508/150‑1 (MQ=17)

                                                                                                                            gcgcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:676173/150‑1 (MQ=255)

                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:816252/150‑1 (MQ=255)

                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  2:4083384/150‑1 (MQ=255)

                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             >  2:3179027/1‑150 (MQ=255)

                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             >  2:3676229/1‑150 (MQ=255)

                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  2:1749129/150‑1 (MQ=255)

                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  1:2385825/150‑1 (MQ=255)

                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            >  1:1142357/1‑150 (MQ=255)

                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGTAGa           >  2:3364043/1‑150 (MQ=17)

                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  2:5878/150‑1 (MQ=255)

                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  1:1400566/150‑1 (MQ=255)

                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  1:2457136/150‑1 (MQ=255)

                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  1:2068169/150‑1 (MQ=255)

                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          >  2:1724900/1‑150 (MQ=255)

                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  1:2373800/150‑1 (MQ=255)

                                                                                                                                  cGAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGact         <  2:3751044/150‑1 (MQ=255)

                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:3286698/150‑1 (MQ=255)

                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        >  1:1281787/1‑150 (MQ=255)

                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        >  1:49753/1‑150 (MQ=255)

                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:3797214/150‑1 (MQ=255)

                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  2:2301957/150‑1 (MQ=255)

                                                                                                                                      acaacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  1:1667249/150‑1 (MQ=255)

                                                                                                                                      acaacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:1207967/150‑1 (MQ=255)

                                                                                                                                       caacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  1:1297519/150‑1 (MQ=255)

                                                                                                                                       caacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  2:105564/150‑1 (MQ=255)

                                                                                                                                        aacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:2631576/150‑1 (MQ=255)

                                                                                                                                        aacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   >  1:3877935/1‑150 (MQ=255)

                                                                                                                                        aacGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   >  1:1084919/1‑150 (MQ=255)

                                                                                                                                         acGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:231072/150‑1 (MQ=255)

                                                                                                                                         acGGAGCCGCTC‑‑GC‑GGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:3539677/150‑1 (MQ=255)

||

CTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073698‑2073987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`