BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	4,344,323	0	G	A	28.7%	256.9 / 18.7	145	intergenic (+54/‑102)	PP_3817 CDS (Polyamine ABC transporter_ ATP‑binding protein)/PP_3818 CDS (Periplamic phosphate‑binding protein)	–/–
Reads supporting (aligned to +/- strand): ref base G (97/5); new base A (41/0); total (140/5)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
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 146 total aligned reads displayed.

AGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACC  >  AE015451_phaM1/4344187‑4344454

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCTTCCCGACCCAACACCCCCCCCGCAGCACACAACC‑TGGcc                                                                                                                              >  2:2617471/1‑136 (MQ=255)

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCAGATCGAGCACCACCCCCACGGAACACCATC‑TTGCa                                                                                                                              >  1:12229443/1‑137 (MQ=255)

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCAGATCGAGAACCCCCCCCCCAGAACACAAAC‑TGGCa                                                                                                                              >  1:4793488/1‑137 (MQ=255)

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCCGATCGAGCACCACCCGCGCGGCACAAGATC‑TGGCa                                                                                                                              >  2:2860739/1‑137 (MQ=255)

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCAACAACCACGCGGAACAAGATC‑TGGCa                                                                                                                              >  1:2103621/1‑137 (MQ=255)

        tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCCGCCAGATCGAGCACCACCCGCGCGGCACAAGAGC‑TGGCa                                                                                                                              >  2:9390942/1‑137 (MQ=255)

             aGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGTCCAATCGAGAACCGCCCCCGCGGCACACGACC‑TTGCaagagc                                                                                                                         >  2:8583864/1‑132 (MQ=255)

             aGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GACTGCCAAATCGAACACCGCCCACACGGCACACAAGC‑TGGCAGgcgc                                                                                                                         >  2:8583863/1‑137 (MQ=255)

              gTGGAAGACGCTTGCGTTTATACCGAATGGGCCGAGAGCGACCGGAGCAAGTCGGCCGGGGCGCATGAACTCTCTTGTACGTTGCAGC‑GCCGGAAAGATCGAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCt                                                                                                                        <  1:7911947/137‑1 (MQ=255)

               tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCCAATCGAGCACCGCCCGCGCAGCACACGATC‑TGGAAAGAGCTg                                                                                                                       >  2:1371700/1‑137 (MQ=255)

                ggCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGT‑GCCCGCCAGAGCGAGAACCACCCGCAAGGAACACGAGC‑TGGGAGGCGCTG‑c                                                                                                                     >  2:17712853/1‑137 (MQ=255)

                ggCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGATGCCCGAGCGAGCACCACCCACACGGCACAAAACC‑TGGCAGGCGCTG‑c                                                                                                                     >  1:12080752/1‑137 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCTTCCAAATCGAGCACCCCCCCCACGGCACACGATC‑TGGCAGGCGCTG‑TAc                                                                                                                   >  1:18050650/1‑137 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCCGAGCGAGCACCACCCCCGCGGCACACGAGC‑CGGCAGGAGCTgtaa                                                                                                                    >  1:6675562/1‑134 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTTTCCGATCGAGCACCACCCGCGCGGCACACAACC‑CGGCAGGCGCTG‑ccc                                                                                                                   >  2:12460348/1‑135 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCTCGATCGAGCACCGCCCGCGAGGCACACGATC‑TGGCAGGCGGTG‑cga                                                                                                                   >  2:16048446/1‑135 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCTCGCACCACCCACACCGCCCACAACC‑TGGAAAGCGACG‑CAc                                                                                                                   >  1:5598062/1‑137 (MQ=255)

                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGACCGAACACCACCCAAACGGCAAACAATC‑TCGCCAGCGCTG‑CAc                                                                                                                   >  2:12460380/1‑137 (MQ=255)

                   aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCACCCACGTAGCACAAGATC‑TCGCCGGCGCTG‑AACt                                                                                                                  >  2:12233566/1‑137 (MQ=255)

                   aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGACCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GACTGCCAGATCGAGAACCACCCGCGTGGCAACCGATC‑TGGCAGGAGCTgtccc                                                                                                                   >  2:4166104/1‑133 (MQ=255)

                       gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCTTCCAGATCAAACACCGCCCGCTCTGCACACAATC‑TGGCAGGTGCTT‑CCCTACTg                                                                                                              >  1:7571838/1‑137 (MQ=255)

                       gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGTCCGAGCGAGCACCACCCCCGCGGAACACAATC‑TGGCAGGGGATG‑CAATACTg                                                                                                              >  1:6568811/1‑137 (MQ=255)

                       gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCCAACCTAACCCCACCCCCGCCGCACACAACC‑TTGCACGCGCCG‑CAcgcccg                                                                                                              >  1:12203579/1‑132 (MQ=255)

                            cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGTCAGGTCGAGGGCCGCCCGCGCGGTACACGATC‑TTGCAGGTGGTG‑TCCTCTTGGGGGc                                                                                                         >  2:6572750/1‑137 (MQ=255)

                              tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGTGCATTAACGATATTGTACGTTGCAGC‑GTCTGTCAGATCGAGCAACACCCACGCAGCACACAAGC‑TTGCAGGCACTG‑CACTATTTGAGGCgg                                                                                                       >  2:14832430/1‑136 (MQ=255)

                              tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTACCAGATCGAGCACCGACCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGt                                                                                                       >  1:17955382/1‑137 (MQ=255)

                              tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GACCTCCCAATCCAACACCGCCCGTACGGCACACGATC‑TTGCAGGCGCTG‑CACTATTTGAggagc                                                                                                       >  2:7335498/1‑134 (MQ=255)

                                    ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTTTCCGGGCGAGCACCGCCCNCGCGGAACACAAGC‑TTGCAAGGGCTT‑TACTCTTTGGGGGGTACGccc                                                                                                 >  1:15245362/1‑137 (MQ=255)

                                    ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATC‑TTGCAGGCGCTG‑CACTATTGGAGGCGTACGccc                                                                                                 <  2:12086867/137‑1 (MQ=255)

                                     cGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAAAACCGACCGCGCGGCACAAGATC‑TGGCAGGGGCTG‑CAATATTGGAGGCATACGcccc                                                                                                >  2:11446671/1‑137 (MQ=255)

                                        aTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCACGTGGCACCCGATC‑TTGCCGGTGCTG‑CACTATTTGCGGAGTACGCCCCGTc                                                                                             >  2:12800576/1‑137 (MQ=255)

                                         tGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCAGAGCGAGAACCACCCACGCGGCACAAAATC‑TGGAAAGAGCTG‑CACTATTTGAGGAGTACACCCCGTcc                                                                                            >  1:15831146/1‑136 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCCTCCCA‑AGCGAGAACCGCCCGCGCGGAAAACGATCTTGTCAGG‑GTTG‑CACTCTTTGAGGAGTCCTCCCCATcac                                                                                           >  1:9364770/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCTTTCCGATTTAGCACCCCCCCCGCTGCCCCCAATC‑CTGCCAGCGCTG‑TACTCTTTTAGGCGTCCACCCCGTccc                                                                                           >  2:15653472/1‑135 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCCTTCAGGACGAACACCCCCCACGCGGCACACAAAC‑TGGCAGGAACTGCCACCTTTTG‑GGCGTACCCCCCGTccc                                                                                           >  2:5264365/1‑135 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCCGTCCGATCGAGCACCACCCCCACGGCACAAAAGC‑TGGCAGGAGCTG‑TACTCTTGGGGGCGTACGCCCCATcac                                                                                           >  2:15629544/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GTCCGCCAGACCAAACACCGCCCCCGCGGCAAACAAGC‑TTGCCGGAGATG‑CACTCTTTGAGGCGTACGCCCCGTcaa                                                                                           >  2:3524877/1‑136 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCCGACCAAACACCACCCACACAGCACAAAACC‑CAGCAGGAGATG‑CACTACTTGAGGAGTACGCCCCCTcac                                                                                           >  1:12526561/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCTGCCAGAGCGAGAACCCCCCACCTAGCAAACAACC‑CTGCAAGAGCTG‑CACTCTTGGAGGCGGACGCCCCATcac                                                                                           >  1:10079404/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCCGCCCGATCGAGCACCGCCCCCGCGGCACACAACC‑TTGTAGGCGCTT‑CACCATTTGAGGAGTACGCCCCTTcac                                                                                           >  2:3411010/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCCGCCCAATCAAGCCACACCCACCCGGCACACAATC‑TGGCAGGCGATG‑TACTCCTTGGGGCCTACACCCCCTcac                                                                                           >  2:5622928/1‑137 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GGCCGCCAAAGCGAAAACCCCCCCCCCGGCACCCAACC‑CTGCAAGCGCTT‑CCCTCCTTGAGGAGTACACCCCCTccc                                                                                           >  1:4151860/1‑135 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTTCCAGATTGAGTACCGCCCACGCGGCACAAAATC‑TTGTAGGCGGTT‑CACCCTTGGAGGCGTGCTCCCCCTcaa                                                                                           >  1:12843228/1‑136 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCCGAACACGGCC‑CGGCCGGCGCTG‑TCATATTTGAGGCCTACACCCCATcca                                                                                           >  2:4678960/1‑135 (MQ=255)

                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GACTGCCCGATCGAGCACCGCCCGCGCGGAAAAAAATC‑TTGCAGGCGCTG‑CACCCTTGGAGGCGTACACCCCGTcac                                                                                           >  2:8998455/1‑137 (MQ=255)

                                            gCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCGTTCCAAATCTAGACCCACCCACGCGGCACAAAAAC‑TGGCAGGCGCTG‑GCCCAATGGAGGCGTAAACCCCATccccc                                                                                         >  2:16499967/1‑133 (MQ=255)

                                               gagaGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCCAACCGAACCCCCCCCCCCCCGCACCCCACC‑TCGCCAGCGCTT‑CACCCCTTGAGGCGTACACCCCCTcaccctcc                                                                                      >  2:1611365/1‑132 (MQ=255)

                                                         tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTTCCAGATCTAGCACCGCCCACGCGGCACACGATC‑TGGCAGGCGCTG‑ACCTCTTGGAGGCGTACACCCCGGCACATTACCATCccacca                                                                            >  2:8150211/1‑137 (MQ=255)

                                                         tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACAAGGTT‑TGGAAGGCGTTG‑TACTACTTGAGGCGTCCGACCCATCACATTGTTATCCCgacc                                                                            >  2:17528641/1‑134 (MQ=255)

                                                         tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAAATCAAACACCACCCACACAGCAAACAATC‑TGGCAGGCGCTG‑CACCACTGTAGGCGTACGCCCCGTCCCACTATCCTCCCGCCa                                                                            >  2:1755397/1‑137 (MQ=255)

                                                         tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCCGCCAAACCGAGCACCACCCCCGCGGCACACAATC‑CGGCCGGCGCTG‑CAACATTGGAGGCATAAGCCCCATCACATTACCATCAcaaca                                                                            >  2:8964918/1‑133 (MQ=255)

                                                              aaGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTTCCCGATCCAACACCCCCCGCGCAGCACCCAATC‑TGGCAAGCACTT‑TACCACTTGAAGCGTACACCCCATAAAAGTATCACCACGCCACTGTc                                                                       >  2:7707491/1‑137 (MQ=255)

                                                              aaGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATTGAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACCATTGGAGGAGTACGCCCCATCACACTCTCATCACGCCACCGTc                                                                       >  2:8380474/1‑137 (MQ=255)

                                                                   ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGACCGAACACCGCCCGCGCGGCACACGATC‑TGGCAGGAGCCG‑CACCACTAGAGGCGTACGCCCCGTCACCTTATCATCCCGCCACTGTCGtcac                                                                  >  2:13932717/1‑134 (MQ=255)

                                                                        gggCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTTCCAGATCGAGCACCGCCCGCGCGGAACACAATC‑TGGCAGGAGATG‑CAACATTGGAGGCGTACGCCCCGTCACATTATCATCCCACCACTTTCGTAACCCGcc                                                             >  2:17238295/1‑137 (MQ=255)

                                                                                     ctctTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATA‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                                <  2:14760570/137‑1 (MQ=255)

                                                                                                   aGC‑GGCTGCCCAATCGCGCACCGCCCACGCGGCACACAACC‑TGGCAGAAGACG‑CACCATTGGAGGCGTACGCCCCATCCCATTATCATCAAGCCACCGGCGTAATCTTCCGGGAAACTCTCTGCCACAGGTCACACg                                  >  2:13149265/1‑137 (MQ=255)

                                                                                                   aGC‑GCCTGCCCAATCTAGCCCCCCCCGCTCCGCACCCAACC‑TGGCAGGCGCTG‑CACTCTTTGAGGCGTCCGCCCCCCCACATTCTCATCCCGCCACCGTCGTAATCTGCCGGCCAACCCTATGCCACCAGTCAACCg                                  >  2:5682922/1‑137 (MQ=255)

                                                                                                   aGC‑GCCTGCCAGATCGAGCACCGTCCGCGCGGCACACGAGC‑TGGGCAGTGCTG‑CACTATTTGAGGCGTACGCCCCATCACACTGTCATCACACCCCTGTTGTCATCCTCCGGCCCACTCTCTGCCAAAGGTCACCCg                                  >  2:11366027/1‑137 (MQ=255)

                                                                                                   aGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGAGC‑TGGCAAGCGCCT‑CACCATTTGAGGCGTACGCCCCATCCAATTCTCCTCCCACCCCTGTCTTCACCTTTCCGCCCAATCTATGTCCCAAGTCAaccc                                  >  2:3211353/1‑136 (MQ=255)

                                                                                                   aGC‑GCCTGCCAGATCGAGCACCGCCCACGCGGAACACAAAC‑TGGAAAGCGCTG‑CCCTCTTTGGGGAGTACACCCCATCACACCCTCATCCCACCACCGTCGTAATCTTCCGGGCAAACCTCTGCCACAAGTCAccca                                  >  2:18746973/1‑136 (MQ=255)

                                                                                                   aGC‑GCCTGCCAGATCGAGCACCACCCGCACGGCACACAATC‑CTGCAAGCGCTG‑CACCACTGGAGGCGTCCCCCCCATCACACTCTCCTCCCACCACCTTCGTCATCTGCCGGCCCACTCTCTGCCACAGGTCAACCg                                  >  2:3211354/1‑137 (MQ=255)

                                                                                                     c‑gcccgccAAATCGAGCACCGCCCACGCGGAACACGATC‑TGGCAGGCGCTG‑CACCATTGGAGGCGTCCGCCCCATCAAATTATCATCAAACCACCGTCGTAACCTGCCGGACAAACCTATGCCAAAAGTCAACCgaa                                >  2:13157264/1‑135 (MQ=255)

                                                                                                     c‑gcccgccAAACCGAACACCACCCACGCGGTACACGACC‑CAGCAAGCGCTA‑TACCATTGGAGGGATACGTCCCGTCACATTATCATCACGCCACCGTTGTAATCTTCTGGTCAACTCTACGCCAAAGGTCGCCCGGc                                >  2:15477267/1‑137 (MQ=255)

                                                                                                     c‑gcCTTCCAGAGCGAGCAACACCCGCGCGGCACACGAGC‑TGGCAGGCGCTG‑CACTATTGGAGGGGTACGCCCCGGCCCATTCTTCTCCCGCCAATGTCGTCATCTTCCTGCCAACCCTCCTCCCCCGGTCCCCCGGc                                >  2:9231162/1‑137 (MQ=255)

                                                                                                                   gAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCNCATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAATCTATGCCACAGGTCAACCGGCATGCTCCGCCTGc                   >  1:17390492/1‑137 (MQ=255)

                                                                                                                   gAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                   >  2:3427307/1‑137 (MQ=255)

                                                                                                                       aCCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt               >  2:3397937/1‑137 (MQ=255)

                                                                                                                         cgcccgcGCGGCACACGATC‑TGGAAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt             >  1:15043768/1‑137 (MQ=255)

                                                                                                                            ccgcGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc          >  2:1882243/1‑137 (MQ=255)

                                                                                                                             cgcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         >  1:9702493/1‑137 (MQ=255)

                                                                                                                             cgcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         >  1:15426710/1‑137 (MQ=255)

                                                                                                                             cgcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         >  2:9190431/1‑137 (MQ=255)

                                                                                                                             cgcgcgACAAACGATT‑GGGCAGCCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         <  2:12806973/137‑1 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGc                                                          <  2:18739198/87‑1 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGc                                                          >  1:18739198/1‑87 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  1:1697855/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:15395666/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:9807972/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:6847558/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:8200901/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  1:7565640/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:3916225/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  1:9097127/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:10067431/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:18360330/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:17769277/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:15312862/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAACCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:18749287/1‑137 (MQ=255)

                                                                                                                              gcgcgGCACACGATC‑CGGCAGGCGCTG‑CACTCTTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        >  2:9947586/1‑137 (MQ=255)

                                                                                                                                  tgCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:2605828/2‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGTCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGACCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:974998/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  2:7038251/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:7892278/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:630446/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:4011872/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:1953028/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:16780144/1‑137 (MQ=255)

                                                                                                                                  ggCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  1:11446953/1‑137 (MQ=255)

                                                                                                                                    cacaCGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc  >  2:9181418/1‑137 (MQ=255)

AGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGC‑GCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATC‑TGGCAGGCGCTG‑CACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACC  >  AE015451_phaM1/4344187‑4344454

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