BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	4,675,336	0	T	G	29.1%	229.3 / 17.5	186	intergenic (+35/+69)	PP_4137 CDS (putative Outer membrane siderophore receptor)/chrR CDS (Chromate reductase)	–/–
Reads supporting (aligned to +/- strand): ref base T (44/85); new base G (0/53); total (46/140)
Fisher's exact test for biased strand distribution p-value = 2.21e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

ACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCG  >  AE015451_phaM1/4675200‑4675472

     gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAtt                                                                                                                                          <  1:4324179/137‑1 (MQ=255)

      aCGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                         <  1:6232484/137‑1 (MQ=255)

         ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGC‑TGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCg                                                                                                                                     <  2:2820561/137‑1 (MQ=255)

               tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCa                                                                                                                                <  2:803369/137‑1 (MQ=255)

               tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCa                                                                                                                                <  1:584390/137‑1 (MQ=255)

               tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCa                                                                                                                                <  1:10334124/137‑1 (MQ=255)

               tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCa                                                                                                                                <  1:15740305/137‑1 (MQ=255)

               tACAACGACCCCGGCCAGGACAATGCCCATTCGACACCGGGCAAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCa                                                                                                                                <  2:16017646/137‑1 (MQ=255)

                accacgacgCCGGCCCGGACAATGCCCGTTCGACACCTGGCGAAGGATTTGGCGTGTATGCCGGGGTTTACTCCAGCTTCCACTGACGGGTGGGCGCGGGGTTCCCCGAGCATTTGTCAGTGAAATTTCCAACGCaa                                                                                                                               >  2:11492662/4‑136 (MQ=255)

                aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAt                                                                                                                               <  1:1674055/137‑1 (MQ=255)

                         ccGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTTCCCCCGACCACGTCCCTGAATTTTCCGCCACACTCGCGGGTa                                                                                                                      >  2:6287900/1‑137 (MQ=255)

                          tgtccAGGACAGTGCCCGTTCGACACTGGGCGAAGGGTTTGGTGGGTGTGCTGGGGTTTCCGAAAGCTCCCTCTAATGGGTGGGAGTGGGTTTACCCGCGAATGCGTCATGGAATTTGCCAACGCATTCGCGGGGaa                                                                                                                     <  2:9485346/134‑1 (MQ=255)

                          cGGCCAGGACAATGCTCGGTCGACCCCGGCCGACGGATTTGGTGGGTGTTCCGGGGTTTCCTACAGCTTCCGTGGACGGGTGGGCGCGGGTTTACCCGGGGGTACGTCGTTGAATTTGCCGACGCGTTCGCGGGGaa                                                                                                                     <  2:6518070/137‑1 (MQ=255)

                          cGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCTGGGTTTCCTACAGCCCCCGCGGGTGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTaa                                                                                                                     <  2:1380556/137‑1 (MQ=38)

                          cGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCTCCCCGAACGAGTCCCTTAATTTTCCTACGCCCTcgcccggta                                                                                                                     >  2:308371/1‑131 (MQ=255)

                          cGGCCAGGACAATGCCCGTTCGACACCGGGCAAAGGATTGGGCGGGTATGCCTGGGTTCCCTACAGCTTCCTCTGACGGGTGGGCGCGGGTTTACCTGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGGaa                                                                                                                     <  2:11163110/137‑1 (MQ=38)

                          cGGCCAGGACAATCCCCGTTCGACACCGGGCGAAGGGTTTGGCGGGTATGCAGGGGTTTCCTACAGCTTCCGCGGATGGGTGGGCGTGGGTTTACCCGCGAGTGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGaa                                                                                                                     <  1:17499542/137‑1 (MQ=255)

                          cGGCCAGAACAATGCCCGTTCGACACTGGGCGAAGGGTTTTGCGTGTATGCCGGGGTTTCCTACAGCTTCATCTGATGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTaa                                                                                                                     <  2:13106311/137‑1 (MQ=37)

                          cGGCAAGGACAGTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTCCCTACAGCTTCCGCTGACGGGTGGGCGTGGGTTTACCCGCGAATGTGTCATTGAATTTGCCGACGCATTCGCGGGGaa                                                                                                                     <  2:4451661/137‑1 (MQ=39)

                          cGGCAAGGACAGTGCCCGTTCGAACCCGGGTGAAGGATTGGGCGTGTATGCCGGGGTTCCCTACAGCTTCCTCGGATGGGTGGGCGCGGGTTTACCCGAGAATCCGTCGTGGACTTTGCCGATGCATTCGTGGGGaa                                                                                                                     <  2:2994186/137‑1 (MQ=255)

                          cGGCAAGGACAATGCCTGTTCGACACGGGGCGAAGGATTTGTTGTGTGTGCCGGGGTTTCCTACAGCTTCCTCTGACGGGTGGGCGCGGTTTTACCCGCGAAGGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGaa                                                                                                                     <  2:12380303/137‑1 (MQ=255)

                            gCCAGGACAATTCCCGTTCGACACTGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTTCCGACACATTCGCGGGGAAAc                                                                                                                   <  2:13013934/137‑1 (MQ=255)

                               agaaCAATTCCCGTTCGACACCGGGCAAAGGCTTTGGTGTGTGTGCCGGGGTTTCCTACAGCTTCCGCGGACGGGGGGGCACGGTTTTACCCGCGAGTTTGTCATTGAGTTTGCCGACGCATTCGCGGGGAAACCCg                                                                                                                <  2:5762309/134‑1 (MQ=255)

                               aGGACAATACCTGTTCGACACTGGGCGAAGGATTTGGCGTGTATGCCGGGGTTCCCCACAGCTCCCGCTGACGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGCCACGTTCGCGGGGAAACCCg                                                                                                                <  2:10765581/137‑1 (MQ=37)

                                       tcccGTTCGACACCGGGCGAAGGGTTTGGCGTGTATGCCGGGTTTTCCTACAGCTCCCGCTGACGGGTGGGCGCGGTTTTACCCGCGAAGGTGTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTaca                                                                                                        <  2:6458137/136‑1 (MQ=38)

                                         ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCGCGAAGGTGTCAGTTAATTCGTCGACGCAATCTCAGGTAAACCCACCCCtcccca                                                                                                      >  2:7709355/1‑132 (MQ=255)

                                               gACACCGGGAGAAGGATTTGACGTGTATGACGGGGTTCCCTACAGCTCCCGCTGACGGGCGGGCGCGGGTTCACCCGCGAATGCGTCAGTGAATTTGTCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCg                                                                                                <  2:16492210/137‑1 (MQ=37)

                                               aacacCGTGCGAAGTTTTTGGCTTGCATTCCTGGTTTTCTTCCACCTTCCTCTGACGGTTGGTCTCGGTTTTATCCGCGAGTGTGTCGTTAACTTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCg                                                                                                <  2:15892036/136‑1 (MQ=255)

                                               aacacCGGGCGAAGGATTGGGCGTGTATGCGGGGGTTTCCGACAGCTTCCTTTGATGGGTGGGCGCGGGTTTACCCACGAATGCGTCGTGGAATTTGCCGACGAGTTCGCGGGGAAACCCGCTCCTACACAGGGGCg                                                                                                <  2:8723034/136‑1 (MQ=255)

                                                 cacaGGGCAAAGGATTTGGCGTGTATTCTGGGGTTTCCGACATCTCCCTCGGACGGGTGGGCGTGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACACGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGtt                                                                                              <  2:9331645/137‑1 (MQ=25)

                                                 cacCGGGCGAAGGATTTGGCGTTTATGCCGGGTTTTCCTAAATCTTCCGCTGACGGGTGGGCACGGGTTTACCCGCGAGTACGTCGTGGAATTTGCCAACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGtt                                                                                              <  2:7100514/137‑1 (MQ=255)

                                                 cacCGGGAAAAGGATTTGACGTGTATGTCGGGGTTCCCTCCAACCTCCGCGGATGGGTGGGCGAGGGTTTACCTGCGAGTGCGTCATGGAGTTTGCCGACACGTTCGCGGGCAAACCCGCTCCTACACAGGGGCGtt                                                                                              <  2:5144104/137‑1 (MQ=255)

                                                   ccGGGCAAAGGATTTGGCGTGCATGCTGGGTTTTCCTACGGCTTCCGCCGATGGTGGGGCGTGGGTTTACCCACGAATTCGTCATTGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGc                                                                                            <  1:11367035/137‑1 (MQ=25)

                                                    cGGGCGAATGGTTTGGCGTGTATGCCGGGGTTTCCTACATCTTCCTCTGACGGTTGGGCGCGGGTTTACCCGCGAATTCGTCATTGAATTTGCCCACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                           <  1:18067060/137‑1 (MQ=37)

                                                     gggAGAAGGGTTTGGCGTGTGTGCTGGGGTTCCCGACAGTTTCCGCGGACGGGGGGGCGTGGGTTTACCCGCGAATGTGTCATGGAATTTTCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                          <  1:1904418/137‑1 (MQ=255)

                                                     gggAGAAGGGTTGGGGGTGGGTGCTGGGGTTTCCGACAGTTCCCGCGGATGGGGGGGCGCGGTTTTACCCGCGAATGCTTCGTTGACTTTGCTAACGCATTCGCGGGGAACCCCGCTCCTACACAGGGGCGTTGcac                                                                                          <  1:15371017/137‑1 (MQ=255)

                                                            gggtttGGTGTGTGTGCTGGGGTTTCCTACATCTTCCTTGGATGGGTGGGTGCGGTTTTACCCGCGAGTGTGTCGTGGAATTTCCCGACGCGTTCGGGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  1:2674685/134‑1 (MQ=255)

                                                            gggtttGGCGTGTATGCCGGGGTTTCCTACATCTTCCGCTGATGGGTGGGCGTGGGTTTACCCGCGAGTGCGTCATTGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  1:2003642/134‑1 (MQ=255)

                                                            ggATTTGGGGTGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  2:18406191/137‑1 (MQ=255)

                                                            ggATTTGGCGTGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  2:1023159/137‑1 (MQ=255)

                                                            ggATTTGGCGTGGATGCCGGGGTTTCCTACAGCTTCCGCTGACGGTTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  2:4044113/137‑1 (MQ=255)

                                                            ggATTTGGCGGGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGTTGAATTGGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                              <  2:15384635/137‑1 (MQ=38)

                                                             gATTTGGCGTGTATGCCGGGTTTTCCTACAGCTTCCGCTGACGGGGGGTCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GAt                                                                             <  2:8228663/137‑1 (MQ=255)

                                                                 tGGCGGGGATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAgg                                                                         <  1:11662398/137‑1 (MQ=38)

                                                                    cGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCCCAAACCCACCCTCGAATTTTCCAACACCCTCTCGGGTAAACCCACTCCCCCACAAAGGCGTTTTACCA‑CAA‑‑‑‑GATTCggggg                                                                      >  2:18712150/1‑136 (MQ=255)

                                                                      tgtAGGCAGGGGTTCCCTACACTTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAGTGCGTCGTGGAGTTTGCCGATGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTc                                                                    <  1:5070809/137‑1 (MQ=255)

                                                                         gtgtcGGGGTTTCCTACATTTCCCGCGGATGTGTGGGCGCGGGTTTACCCGTGAGTGTGTCATGGAGTTTTCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGa                                                                 <  1:7316993/133‑1 (MQ=255)

                                                                         gggctgggtTTTCCTAAAGTTTCTGTTGAGGGGTGGGTGTGGGTTTACCTGGGAGTGTGTCGTGGAGTTGGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGa                                                                 <  1:5390998/132‑1 (MQ=255)

                                                                           gCCGGGTTTTCCTACAGCTTCCGCTGATGGGGGGGCGTGGGTTTACCCGCGAATGCGTCGTGGAGTTTGCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                               <  1:3922992/137‑1 (MQ=37)

                                                                           gCCGGGGTTTCTTACATCTTCCTCTGACGGTTGGGCGCGGGTTTACCTGCGAGTGCGTCGTGGACTTTGCCGACGCGTTTGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                               <  1:11135919/137‑1 (MQ=255)

                                                                           gCCGGGGTTTCCTACAGCTTCCGTTGACGGGTGGGCGTGGGTTTACCCGCGAGTGTGTCGTGGAATTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                               <  1:4595407/137‑1 (MQ=255)

                                                                           gCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCTCGAAGGCATCCTTTAACTTTCCCAAACATTCTCCGGTCACCCCCCCTCCCCCCCGTGGCGGTTCACAGACAG‑‑‑‑GATCAGGGTTTCaacc                                                                >  2:9342467/1‑134 (MQ=255)

                                                                            ccGGGGTTTCCTACAGCTTCCGCTAAAGGGGGGTTGCGGGTTTCCCCGCGAATGTGTCGGTGAATTGGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCg                                                              <  1:13146583/137‑1 (MQ=255)

                                                                               gtgTTTCCTACGTCTTCTGCTGACGGGGGGGTGGGGGTTTACCCGAGAATGCGTCATGAAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGccc                                                           <  1:5659926/135‑1 (MQ=255)

                                                                               gggTTTCCTACGGCTTCCGCGGAGGGGGGGGCGGGGGTTGACCCGAGAATGCGTAGTGGAGTTTTCAAACGCATTCGTGGGAAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGccc                                                           <  1:848119/137‑1 (MQ=255)

                                                                                ggTTTCCTACGGCCTCCGCTGACGGGTTGTCGTGGTTTTACCCGCGAATGCGTCCTTGACTTTACCGCCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  2:17368107/137‑1 (MQ=37)

                                                                                ggTTTCCTACATCTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAATGTTTCATTGAATTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  2:6934058/137‑1 (MQ=38)

                                                                                ggTTTCCTACAGTTTCCGCTGACGGGTGGGCGCGGGTTTACCCGAGAATGCGTCATTGAGTTTTCCGACGCATCCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  1:12295457/137‑1 (MQ=255)

                                                                                ggTTTCCTACAGCTTCCTCTGATGGGTGGATGCGGGTTGACCCGCGAAGGCGTCATTGAGTTGGCCGATGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  2:1327540/137‑1 (MQ=255)

                                                                                ggTTTCCTACACCTTCCTCTGACGGGTGGGCGGGGGTTTACCAGCGAATGTTTCGTTGAGTTTTCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  2:2119720/137‑1 (MQ=255)

                                                                                ggTTCCCAACAGCTCCCGCGGACGGGTGGGCGCGGTTCTACCCGCAAAGGCGTCGTTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                          <  2:14291994/137‑1 (MQ=37)

                                                                                            ctccctcTGACGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCGTTCGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                              <  2:3383991/131‑1 (MQ=38)

                                                                                            cTTCCGCTGACGGGTGGGCGTGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGAATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                              <  2:6595131/137‑1 (MQ=255)

                                                                                                   tGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                       <  2:7418738/137‑1 (MQ=255)

                                                                                                   tGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAGTTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTCGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                       <  1:15851409/137‑1 (MQ=38)

                                                                                                   tGACGGGGGGGCGTGGGTTTACCCGCGAATATGTCATTGAATTTTCCGACGCATTCGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                       <  1:14784103/137‑1 (MQ=255)

                                                                                                      tggggggGCGTGGGTTTACCCGCGAATGTGTCGTGGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                    <  1:11964961/136‑1 (MQ=38)

                                                                                                      cGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCTGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                    <  1:10479509/137‑1 (MQ=255)

                                                                                                      cGGGTGGGCACGGTTTTACCCGCGAATTCGTCGTTGAATTTGCCGACGCATTCGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                    <  2:13319782/137‑1 (MQ=38)

                                                                                                        ggtgggtgGGGGTTTACCCGCGAGTGTGTCGTGGAGTTTGCCGACGCGTTTGCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                  <  1:3296426/137‑1 (MQ=255)

                                                                                                          tGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAcc                                <  2:1795709/137‑1 (MQ=255)

                                                                                                                 gggTTTACCCGCGAGTGCGTCGTGGAGTTTGCAGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGAc                         <  1:4194827/137‑1 (MQ=255)

                                                                                                                  ggTTTACCCGCGAATGCTCCATTGAATTTGCCGACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACa                        <  2:6154905/137‑1 (MQ=255)

                                                                                                                         ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTCCACAAGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                 >  1:6525214/1‑137 (MQ=255)

                                                                                                                          cgcgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                >  1:12950778/1‑137 (MQ=255)

                                                                                                                           gagAATGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGc               <  1:10698588/135‑1 (MQ=255)

                                                                                                                             aagtGCGCCATTGAATTTTCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGTGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCaa             <  2:8968193/134‑1 (MQ=255)

                                                                                                                               aTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGa           >  1:5934262/1‑137 (MQ=255)

                                                                                                                               aTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCCGTCGACAAACGCCTGCAAGa           >  2:6651636/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAgca          >  1:4811959/1‑135 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  1:8173921/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  1:668151/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  1:16897297/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  1:8717339/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  1:9351370/1‑137 (MQ=255)

                                                                                                                                tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCACTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa          >  2:3812360/1‑137 (MQ=255)

                                                                                                                                 gtgTCGTTGACTTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc         <  1:6579640/135‑1 (MQ=255)

                                                                                                                                 gtgTCGTTGACTTTGCCGACGCATTCGAGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc         <  1:7680368/135‑1 (MQ=255)

                                                                                                                                 gCGTCGTTGAGTTTGCAGAAGCGTTAGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc         <  1:4942020/137‑1 (MQ=255)

                                                                                                                                 gCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc         >  1:18071369/1‑137 (MQ=255)

                                                                                                                                 gCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc         >  1:1238650/1‑137 (MQ=255)

                                                                                                                                     cgttGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc     <  1:15134824/135‑1 (MQ=255)

                                                                                                                                     cATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc     <  1:1715925/137‑1 (MQ=255)

                                                                                                                                     cATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc     <  2:13195995/137‑1 (MQ=255)

                                                                                                                                       ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc   <  2:17314145/137‑1 (MQ=255)

                                                                                                                                       ttGAATTTGCCGACGCATTCGCGGGGNAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc   <  1:18626527/137‑1 (MQ=255)

                                                                                                                                       ttGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc   <  1:3894894/137‑1 (MQ=255)

                                                                                                                                        tGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAGGATTGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGt      >  2:13289899/1‑137 (MQ=38)

                                                                                                                                        tGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg  <  2:13957107/137‑1 (MQ=255)

                                                                                                                                        tGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg  <  2:17199821/137‑1 (MQ=255)

ACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCG  >  AE015451_phaM1/4675200‑4675472

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