BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	5,671,671	0	C	T	27.1%	519.8 / 33.4	315	intergenic (+90/‑80)	PP_4978 CDS (conserved protein of unknown function)/PP_4979 CDS (putative Periplasmic binding protein)	–/–
Reads supporting (aligned to +/- strand): ref base C (88/138); new base T (0/84); total (88/227)
Fisher's exact test for biased strand distribution p-value = 3.63e-15
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 324 total aligned reads displayed.

GGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTG  >  AE015451_phaM1/5671535‑5671807

ggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGcc                                                                                                                                          <  2:17677706/137‑1 (MQ=255)

    cAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGgtg                                                                                                                                      <  1:2950468/137‑1 (MQ=255)

             cGAGGCGCATCTGCTGCGATCCACCCAATCCTGATGTGCCCGCAAAACTCATTTTCTTGCGCGGTCTTGTGGGAGCGGGCGTTCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGTGGGc                                                                                                                             <  2:3561091/137‑1 (MQ=255)

               aGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGGGTGCCCAAAAAAAGAAAAACCCGATCGTTGGGACCGGCTTCGTCGGTGGTCGggggtaa                                                                                                                           >  2:8856265/1‑134 (MQ=255)

                     aTCTGCTGGGATCCACCCAATCCGGGTGTGCCTGCAAAACCCATTTTTTTACGCGGTCTTTGGGGAGTGGGCGTGCCCGCGAAGAGAGAACCCCAGTCGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGc                                                                                                                     <  2:13234347/137‑1 (MQ=21)

                        tgctgcGATCCACCCAAGCCTGATGTGCCCGCAAAACCCATTGCCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                  <  2:12269998/137‑1 (MQ=255)

                        tgctgcGATCCACCAAAGCTTGATGTGCCCGAAAAACTCATTTTCTTTCGCGTTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                  <  2:4269002/137‑1 (MQ=37)

                         tcttcGATCCACCCAAGCCTGGTGTGCCCGCAAACCTTATTGTCTTGCTCGTTCTTGTGGGATCGTTCGTCCCCGCGAAGAAAGCAACCCGATCTTTGGCACCGGCTTCGCCGGTGTGCGCGGGCACGCCCGCTccc                                                                                                                 <  2:13951274/133‑1 (MQ=25)

                                    cccAATCCTGCTTTTCCCGAAAACCCCATTTTCTTGCTCGTTCTTGTGGGAGCGGGCGTGCCCGCGAAGAAAACAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGTGGACACGCCCGCTCCCACATCGATAat                                                                                                      <  2:8214305/137‑1 (MQ=21)

                                    cccAAGCTTGATGTCCCCGCAAAACTCGTGTTTTTGCGCGGTCTGGTGGGAGCGGGCGTGCCTGCGAAAAGAGAAACCCGATCGGTGGCACCGGCTTCGCTGGTGTTCGAGGGCACGCCCGCTCCCACATCGATAat                                                                                                      <  2:751654/137‑1 (MQ=25)

                                     ccAACCCTGATTTGCCCGCAAAACCCATTGTCTTATTCGGCCTTGTTGGAGCGAGCGTGCCCGCGAAGAGAGCAACCCGGTCGTTGGCACCGGCTTCTCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatg                                                                                                     <  2:18641663/137‑1 (MQ=25)

                                       aaGCTTGGTGTGCCCGCAAAACTCATTGTCTTGCGTGGTTTTGTGGGAGTGGGTGTGCCCGAGAAGAGAGCAACCCGATTGTTGGCACCGGTTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgca                                                                                                   <  2:14907126/137‑1 (MQ=25)

                                         tccgggtgtgCCCGAAAAATTCATTGTTTTGTGTGGTTTTGTGGGAGTGGGTGTGCCCGTGAAGAGAGCAACCCGATCGTTGGCACTGGCTTCGCTGGTGTTCGTGGGCACTCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:6155501/131‑1 (MQ=255)

                                         tccTGATTTGCCCGAAAAACTCATTTTCTTGCGCGGTCTTGTGGGAGTGGGCGTGCCCGCGAAGAGAACAACCCGATCGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:12786295/136‑1 (MQ=37)

                                         tccTGATGTGCCCGCAAAACTCGTTGTCTTGCGTGGTCTGGTGGGAGCGGGTGTGCCTGCGAAGAGAGCAACCCGATCGTGGGCACCGGCTTCGCCGGTGTTTGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:6573659/136‑1 (MQ=37)

                                         tccTGATGTGCCCACAAACCCCATTGTTTTGCGTGTTTTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGGTGGCACCGGCTTCGCTGGTGTTCGCGGACACACCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:5914567/136‑1 (MQ=21)

                                         tccTGATGTACCCGCAAACCCCATTTTCTTCCGCGGTCTTGTGGGAGCGGGCGTACCCGCGAAGAGAGAAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:8060944/136‑1 (MQ=37)

                                         gctgggtgtgCCCGCAAAACTCGTTGTCTTGCGCGGTTTTGTGGGAGCGGGCGTGCCCGCGAAGAAAGCACCCAGATTGTTGGCACCGGTTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:7836913/131‑1 (MQ=25)

                                         gctgggtgtgCCAGAAAAACTCGTTGTCTTGCGCGTTTTTGGGGGAGTGGGTGTGCCCGCGAAGAGAGAAACCCGGTTGTTGGCACTGGTTTCGCTGGTGTTTGTGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:10279137/131‑1 (MQ=255)

                                         gctggATGTACCCAAAAAACTCGTTGTTTTGCGTGGTTTGGGGGGAGTGGGTGTGCCCGCGAAGAGAGCAACCCGATTGTTGACACTGTTTTCGCCGGTTTTCGCGGGCACCCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:15437849/133‑1 (MQ=255)

                                         gCCTGGTGTGCCCGCAAAACTCGTTGTCTTGCGTGGTTTTGTGGGAGCGGGTGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:15871322/137‑1 (MQ=37)

                                         gCCTGGTGTGCCCGCAAAACTCATTTTTTTGTGCGGTTTTGTGGGGGCGGGCGTTCCAGAGAAAAGAGCAACCCGATCGTTGACACCGGTTTCGCCGGTGTTCGCGGACACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:8782433/137‑1 (MQ=255)

                                         gCCTGGTGTGCCCGAAAAACTCATTGTTTTGCGTGGTCGTGTGGGGGCGGGCGTGCCCGCGAAGAGGGCACCCCGATTGTTGGCGTTGGCTTCGCTGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:9999365/137‑1 (MQ=25)

                                         gCCTGATGTGCCCGCAAAACTCATTGTCTTGCGTGGTCTGGGGGGAGCGGGCGTGCCTGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGGGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:13404038/137‑1 (MQ=38)

                                         gCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCACGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:1239637/137‑1 (MQ=255)

                                         gCCTGATGTCCCCAAAAAACTCATTGTCTTGCGCGGTTTTGTGGGAGCGGGCGTACCCGCGAAGAGAACAACCCAATCGTTGGCACCGGCTTCCCCGGTGTTCGCGGGCACGCCCACTCCCACATCGATAatgcatg                                                                                                 <  2:10184608/137‑1 (MQ=37)

                                         cccTGATGTTCCGGCAAAACTCGTTGTCTTTTTCGGTTTGGTGGGGGTGGGTTTGCCCGCGAAGAGAGCACCCCGGTCTTTGGCACCGTCTTCGCCGGTGTTCTCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:7641682/136‑1 (MQ=18)

                                         cccTGATGTTCCCGCAAACCCCATTGCCTTTCGCGTTCTTGTGGGAGCGGGCTTCCCCGCGAAGAAAGCAACCCGATCGTTGGCACCGGCTTCGCCGTTGTTCGCGGACACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:8734073/136‑1 (MQ=255)

                                         cccTGATGTACCAACAAAACTCATTTTCTTGCTCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGTTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:14557850/136‑1 (MQ=37)

                                         acttGATGTGCCCGAAAAACTCATTTTTTTGCGTGGTTTGGTGGGAGTGGGCGTGCCAGAGAAGAGAGCAACCCGATTGTGGGCACCGGTTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  1:4145983/134‑1 (MQ=25)

                                         accggATGTGCCCGCAAACCTCATTTTTTTTCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                                                 <  2:2835786/133‑1 (MQ=37)

                                          ccgggtgtgCCCAAAAAACTCATTTTTTTGTGTGTTTTTGGGGGAGTGGGTGTGCCCGAGAAGAGAGCAACCCGGTTGTGGACACCGGTTTCGCTGGTGTTCGCGGGCACCCCCGCTCCCACATCGATAatgcatgc                                                                                                <  1:16939946/132‑1 (MQ=255)

                                               tgtgCCCGCAAAACTCATTGTCTGTCGCGGTCTTGTGGGAGCGGACGTGCCCGCGAGGAGAGCAACCCGATCGTTGGCACTGGCTTCGCCGGGGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTa                                                                                           <  2:2630366/137‑1 (MQ=37)

                                                                 tGTCTTGCGCGGTCTTGTGGGAGCGGGAGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgcg                                                                         <  2:9663566/137‑1 (MQ=255)

                                                                          cgGTCTTGTGGGAGCGGTTGTGCCCGCGAAGAGAGCAACCTGATTGTGGGCACTGGCTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAAtctc                                                                <  2:6036913/137‑1 (MQ=37)

                                                                            tttttGGGGGGGGGGGTGTTCCCGCGAGGAGAGCAACCAGATCGTGGGCACCGGCTTCGCTGGTGTTCGCGGACACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  1:13279441/134‑1 (MQ=255)

                                                                            tttttGGGGGAGGGGGGGTGCAAGCGAAGAGAGCAACCTGATTGTGGACACTGTCTTCGCTGTTGTTCTTGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  1:3452432/134‑1 (MQ=255)

                                                                            gttttGGGGGAGTGGGTGTACCGGAGAAGAAAACAACCCGATTGTGGACACCGGCTTCGCTGGTGTTCGTGGGCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  1:559696/134‑1 (MQ=255)

                                                                            gTCTTGTGGGAGCGGGCGTGCCCGCGAGAAGAGCAACCCGATCGTTGGCACAGGCTTCGCCGGTGTTTACGGGCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  2:12568801/137‑1 (MQ=38)

                                                                            gTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAACAACCCGATCTTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  2:9901934/137‑1 (MQ=255)

                                                                            gTCTTGTGGGAGCGGGCGTCCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                                              <  2:1613872/137‑1 (MQ=255)

                                                                                tgtgGGAGTGGGCGTGCCCGCAAAGAGAACAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  1:5638078/137‑1 (MQ=255)

                                                                                tgtgGGAGGGGGTGTGCCCGCGAAGAGAACAACCCGCTTGTTGGCACCGGCTTCGCCGGTGTTTGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  1:12016691/137‑1 (MQ=38)

                                                                                tgtgGGAGCGGGTGTACCAGCGAAGAGAGCAACCCGGTTGTGGGCACTGTTTTCGCTGGTGTTCGTGGGAACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  1:14924840/137‑1 (MQ=255)

                                                                                tgtgGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTTTCCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  2:2247513/137‑1 (MQ=255)

                                                                                tgtgGGAGCGGGCGTGCCCGAGAAGAGAAAAACCCGATCGTTGGCACCGGTTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  2:12870510/137‑1 (MQ=255)

                                                                                tgggggAGGGGGCGTGCCCGTGAAGAGAGAAACCCGGTTGTGGGCACTGGTTTTGCTGGTTTTTGCGGCCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  1:17388451/134‑1 (MQ=255)

                                                                                ggggggAGGGGGCGTGCCAGAGAAGGGAGAAACCTGATTGTGGGCACTGGTTTCGCTGTTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                                          <  1:8221446/134‑1 (MQ=37)

                                                                                       gtgggCGTGCCCGCGAAGAGAGCAACCCGATTGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCcaca                                                   <  2:722853/135‑1 (MQ=255)

                                                                                         gttcGTTCCAGCGAAGAGGGCACCCTGATTGTGGGCACCGGCTTCGCCGGTGTTCGCGGGCATACCCGCTCCCACATCGATAATGCATGCCTTTACCTTACACCGGCCGAGCGCGTAATCTCCTGCCATTCCACata                                                 <  2:11896249/134‑1 (MQ=25)

                                                                                         gggCGCGCCCGCGAGGAAAGCAACCCGATCGTTGGCACCGTCTTCGCCGGTGTTTGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACata                                                 <  1:10001798/137‑1 (MQ=38)

                                                                                          ggCATGCCCGCGAAGAGAGCAACCCGACCGTTGGCACCAGCTTCACCGGTGTTCGCGGGGAAGCACGCTCCCACATCGATGATGCAGGCCTTTCCCGTACACCAGCCCAACGCGTAATCCCCAGCCACGCCACacat                                                >  2:18359465/1‑137 (MQ=255)

                                                                                              tGCCCGCGAAGAGAACAACCCGATTGTTGGCACTGGCTCCGCTGGTGTTCGCGGCCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:903337/137‑1 (MQ=255)

                                                                                              tGCCCGAGAAGAGAGAAACCCGATTGTTGGCACTGGTTTCGCCGGTGTTCGCGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:14353924/137‑1 (MQ=38)

                                                                                              tGCCCGAGAAGAGAGAAACCCGATTGTGGGAACTGGTTTCGCTGGTGTTCGCGGACACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:18639768/137‑1 (MQ=255)

                                                                                              tGCCAGGGAAGGGAGCAACCTGATTGTTGGCACTGGTTTCGCTGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:1352219/137‑1 (MQ=38)

                                                                                              tGCCAGCGAAGAGAGAAACCCGATTGTGGGCACTGGTTTTGCTGGTGTTCGGGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  2:11139/137‑1 (MQ=255)

                                                                                              ggCCTGCGAAGAGAGCAACCCGGTTGTGGGCACCGGTTTCGCCGGTGTTCGCGGACACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:4575409/136‑1 (MQ=255)

                                                                                              ggCCCGCGAAGAGAGCAACCCGATCGTTGGCACTGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  2:13820101/136‑1 (MQ=255)

                                                                                              cacccGCGAAGAGGGAAACCGACTTTTTGGCACTGGTTTTGCTGTTGTTTGTGGGCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                                            <  1:16739856/135‑1 (MQ=255)

                                                                                                  cgcaAAGAGAACAACCCGATCGTTGACACTGGTTTCGCCGTTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTc                                        <  2:16469357/137‑1 (MQ=39)

                                                                                                   gcgAAGAGAGCAACCAGGTCGTTGGCACCGGCTTCGCCGGTTTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGACCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCt                                       <  2:880499/137‑1 (MQ=255)

                                                                                                     gAAGAGAGAAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGt                                     <  2:4505922/137‑1 (MQ=255)

                                                                                                       agagagCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTc                                   >  1:17940022/1‑137 (MQ=255)

                                                                                                        gagagCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCt                                  >  1:10224653/1‑137 (MQ=255)

                                                                                                        gagagCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCt                                  >  1:7343602/1‑137 (MQ=255)

                                                                                                             aaaCCTGATCTTTGGCACCGTCTTCCCAGGGGTTCGCGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATc                             <  1:10778393/136‑1 (MQ=255)

                                                                                                              aaCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATcc                            >  1:2406978/1‑137 (MQ=255)

                                                                                                              aaCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATcc                            >  1:3837148/1‑137 (MQ=255)

                                                                                                              aaCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATcc                            >  2:11772733/1‑137 (MQ=255)

                                                                                                              aaCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATcc                            >  1:14536954/1‑137 (MQ=255)

                                                                                                              aaCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGACCGTCTGTTCTGGATcc                            >  1:678759/1‑137 (MQ=255)

                                                                                                               aCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGAACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCTGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCct                           <  2:3381388/137‑1 (MQ=255)

                                                                                                                cccGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCctg                          >  2:15737224/1‑137 (MQ=255)

                                                                                                                   gATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCctgctg                       >  2:4862487/1‑137 (MQ=255)

                                                                                                                   gATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCctgctg                       >  2:3772511/1‑137 (MQ=255)

                                                                                                                    aTCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCtgctgc                      >  1:10457499/1‑137 (MQ=255)

                                                                                                                     tagTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTgctgct                     <  2:14108867/135‑1 (MQ=255)

                                                                                                                     tCGTTGGCACCTTCTTCGCTGTTGTTTGCGGCCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTgctgct                     <  1:13964732/137‑1 (MQ=255)

                                                                                                                        ttGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGt                  >  2:9377595/1‑137 (MQ=255)

                                                                                                                        ttGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGt                  >  2:15307653/1‑137 (MQ=255)

                                                                                                                        ttGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGt                  >  1:7039563/1‑137 (MQ=255)

                                                                                                                         gggCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGtt                 <  1:11680818/136‑1 (MQ=255)

                                                                                                                             aCCGGCTTCGCCGTTGTTCGCGGGCACGCCCGCCCCCACATCGATAATGCATGCCTTCACCGTACACCGGCCGAGCGCGTAGTCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCCGCTCGTTTGcc             >  2:10960605/1‑137 (MQ=39)

                                                                                                                             aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGcc             >  2:10617208/1‑137 (MQ=255)

                                                                                                                             aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGcc             >  2:11492658/1‑137 (MQ=255)

                                                                                                                             aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGcc             >  1:3024014/1‑137 (MQ=255)

                                                                                                                             aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGcc             >  2:12337082/1‑137 (MQ=255)

                                                                                                                              ccGGTTTCGCGGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCt            <  2:685586/137‑1 (MQ=255)

                                                                                                                              ccGGCTTCGCTGGTGTTCGCGGGCACACCCGCTCCCACATCGATAATNCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCt            <  1:15075111/137‑1 (MQ=255)

                                                                                                                              ccGGCTTCGCCGGTGTTCGGGGACACGCCTGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCt            <  2:1750778/137‑1 (MQ=255)

                                                                                                                              ccGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCt            <  2:13141330/137‑1 (MQ=255)

                                                                                                                                  cTTCTCTGGTTTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGc        <  1:10383794/137‑1 (MQ=255)

                                                                                                                                  cTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGc        <  2:16690153/137‑1 (MQ=255)

                                                                                                                                  cTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTATGTTCTGGATCCTGCTGCTCGTTTGCCTGAGc        >  2:1720047/1‑137 (MQ=255)

                                                                                                                                    tCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGccc      >  2:11528706/1‑137 (MQ=255)

                                                                                                                                    tCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGccc      >  2:16334363/1‑137 (MQ=255)

                                                                                                                                     cGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTACGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCg     >  2:3237011/1‑137 (MQ=255)

                                                                                                                                        cGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGANCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg  >  1:17184839/1‑137 (MQ=255)

                                                                                                                                        cGGTGTTCGCGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTg  <  1:13752104/137‑1 (MQ=255)

GGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGGCGTGCCCGCGAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTCGTTTGCCTGAGCCCGTTG  >  AE015451_phaM1/5671535‑5671807

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