BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	5,671,653	0	C	T	30.9%	406.9 / 34.5	278	intergenic (+72/‑98)	PP_4978 CDS (conserved protein of unknown function)/PP_4979 CDS (putative Periplasmic binding protein)	–/–
Reads supporting (aligned to +/- strand): ref base C (51/139); new base T (0/85); total (51/227)
Fisher's exact test for biased strand distribution p-value = 1.19e-09
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 288 total aligned reads displayed.

CTATCGCTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTC  >  AE015451_phaM1/5671519‑5671789

  aTCGCTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAANCTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTg                                                                                                                                      <  1:18024404/137‑1 (MQ=255)

  aTCGCTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTg                                                                                                                                      <  1:17071004/137‑1 (MQ=255)

     gCTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGGGCGGG‑GGTGCCCCC‑‑‑CAAGGGAACAACCCCATCtttttta                                                                                                                                   >  2:17154950/1‑133 (MQ=255)

      cTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGGGGG‑GGTGCCCGC‑‑‑TAAAAGAGCAAACCCATCttttgtaa                                                                                                                                  >  2:9081174/1‑132 (MQ=255)

      cTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCAc                                                                                                                                  <  2:15464828/137‑1 (MQ=255)

             aGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGGGCGGG‑GGTGTCCAG‑‑‑GAAGAGAGCAACCCCATCGTTGGGACCAGTTTc                                                                                                                           >  2:14811863/1‑137 (MQ=255)

              gcggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCg                                                                                                                          <  1:11768991/137‑1 (MQ=255)

              gcggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCg                                                                                                                          <  2:13271138/137‑1 (MQ=255)

              gcggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCg                                                                                                                          <  1:8470570/137‑1 (MQ=255)

              gcggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCg                                                                                                                          <  1:15466470/137‑1 (MQ=255)

              gcggcgCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCg                                                                                                                          <  1:17584653/137‑1 (MQ=255)

                 gcgcAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCg                                                                                                                       <  2:11115229/137‑1 (MQ=255)

                    cAAGGATATCGAGGTACATTTGCTGCGATCCACCCAAGCCTGATGTGCCCACAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCTGCTTCGCCGgtg                                                                                                                    <  2:7566336/137‑1 (MQ=39)

                        gATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGgtgccaa                                                                                                                <  2:8114326/137‑5 (MQ=255)

                               aGGCGCAGCTGCTTCGATCCACCCAATCTTGATTTGCCCGCAAAACTCATTGTCTTGTTTGGTCTTGTGGGAGTGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGTTTTTCGTGGGCAc                                                                                                         <  2:8070087/137‑1 (MQ=255)

                                 gcgcGTTTGTTGCGATCCCCCCAATCCTGGTGTGCCCACAAACCTCCTTTTTTTTCGAGGTCTTGGGGGAGTGGG‑CGGGCCCGC‑‑‑GAAGAGAGCAACCCGATTGTTGGAACTGGTTTCGCGGGTGTTTGTGGGCACGc                                                                                                       <  2:9371640/137‑1 (MQ=255)

                                     aTCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATTGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGc                                                                                                   <  2:651734/137‑1 (MQ=255)

                                        tgttgagATCCACCCAAGCCGGATGTGCCCACAAACCCCGTTTTCTTGCGCGGTCTTGTGGGAACGGG‑CGTGCCCGC‑‑‑GAAGGGAGCAACCCGATCGTTGGCACTGGCTTCGCTGGTTTTCGCGGGCACGCCCGCTcc                                                                                                <  2:3921679/131‑1 (MQ=25)

                                         gcggcGACCCACCCAACCTGGAGGTTCCCGCAAAACTCATTGTCTTGCGCGGTTTTGGGGGAGTGGG‑TGTGCCCGA‑‑‑GAGGAGAGAAACCCGATTGTTGGCACTGGTTTCGCCGGTGTTCGTGGGCACGCCCGCTccc                                                                                               <  2:18392668/134‑1 (MQ=255)

                                         acgacGATCCACCCACGACCGATGTTCCCACAAAACCCAGTGTCTTGCGCGGTCTGGTGGGGGCGGG‑CGTGCCCGC‑‑‑GAAGAGAACAACCCAATCGTGGGCACCGGCTTCGCCGGTGTTCGCGGACACACCCGCTccc                                                                                               <  2:15089590/133‑1 (MQ=25)

                                             cGATCCACCCAAGCCTGGTGTTCCCGCAAACCTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑TGTGCCCGC‑‑‑GAAGAGAGCAACCCGATTGTTGGCACTGGTTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACAt                                                                                           <  2:5495629/137‑1 (MQ=37)

                                             cGATCCACCCAAGCCTGATTTGCCCGCAAAACTCATTGTCTTTCACGGTCTGTTGGGAGCGGG‑CGTCCCCGC‑‑‑GAAGAGAGCAACCCGCTCTTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAt                                                                                           <  2:4307283/137‑1 (MQ=25)

                                                  cACCCAAGCCGGATGTGCCCGCAAAACTCATGGTCTTGCGCGGTTTTGTGGGAGCGGT‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCACCTGTGTTCGCGGGCACGCCCGCTCCCACATCGATa                                                                                      <  2:3759893/137‑1 (MQ=38)

                                                  aaCCCAAGCCTGATGTGCCAGAAAAACTCATTGTTTTGCGCGGTTTGGTGGGAGCGGG‑TGTGCCCGG‑‑‑GAAGAGAGCAACCCGGTCGTGGGCACCGGCTTCACCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATa                                                                                      <  1:13867102/136‑1 (MQ=25)

                                                     ccAAGCCTGATGTTCCCGCAAAACTCGTTGTCTTGCGTGGTCTTGTGGGAGTGGG‑TGTGCCTGCGGAGAAGAGAGCACCCCGATTGTTGGCACCGGTTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATa                                                                                      <  2:17109758/137‑1 (MQ=25)

                                                       aaGCTTGATGTGCCCACAAAACTCATTGTCTTGCGCGTTTTTGTGGGAGTGGG‑CGTGCCTGC‑‑‑GAAGAGAGAAACCCGATTGTTGGCACCGGCTTCGCCGGGGTTCGCGGGCACACCCGCTCCCACATCGATAatgca                                                                                 <  2:708380/137‑1 (MQ=255)

                                                       aaGCCTAAGGTGCCCGAAAAACCCATTTTCTTGCGCGGCCTTGTGGGAGTGGG‑CGTGCCCGC‑‑‑GAAGGAATCAACCCGCTCGTTGGCACCGGCTCCGCCGTTGTTCTCGGCCACGCCCGCTCCCACATCGATAatgca                                                                                 <  2:14017928/137‑1 (MQ=255)

                                                         tcttggtgtaCCCACAAAATTCGTTGTCTTGTGTGGTTTTGTGGGAGTGGG‑CGTGCCCGA‑‑‑GAAGAGAACAACCCGATTGTGGGCACCGGCTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  1:14405189/131‑1 (MQ=21)

                                                         tcttGATGTGCCCGAAAAACTCTTTGTCTTGCGCGTTTTTGTGGGAGCGGG‑TGTGCCCGA‑‑‑GAAGAGAGCAACCCGATCGTGGGCACCGGCTTCGCCGGTGTTCGCGGACACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:18650838/134‑1 (MQ=37)

                                                         tctggATGTTCCCGCAAAACCCCTTGTCTTGCACGGTCTTGGGGGAGCGGG‑CGTACCCGC‑‑‑GAAGAGAGCAACCAGATCGTTGGCACCGGCTTCTCTGGTGTTCGTGGGCACTCCCGCTCCCACATCGATAatgcatg                                                                               <  2:10432794/133‑1 (MQ=255)

                                                         tccTGATGTGCCCACAAAACTCATGTCCTTGCGTGGTTTTGTGGGAGCGGG‑CGTGCCTGC‑‑‑GAGGAGAGCACCCCGATTGTTGGCACCGGCTTCGCTGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:13839986/136‑1 (MQ=25)

                                                         gccgggtgtgCCCGAAAACCTCGTTGTCTTGCGTGGTTTGGGGGGGGGGGG‑TGTGCCGGA‑‑‑GAGGAGAGAAACCCGGTTGTGGGCACTGGTTTCGCTGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:6106627/131‑1 (MQ=255)

                                                         gCCTGGTGTGCCCGCAAAACTCATTTTCTTGCGCGGTCTTGGGGGAGTGGG‑CGTGCCCGC‑‑‑GAAGAGAACAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:15655438/137‑1 (MQ=38)

                                                         gCCTGGTGTGCCCGAAAAACGCGTTGTCTTGTGTGGTTTTGGGGGGGTGGG‑TGTTCCCGA‑‑‑GAAGAGAGCAACCCGATTGTGGGCACTGGCTTCGCTGGTGTTTGTGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  1:18481600/137‑1 (MQ=255)

                                                         gCCTGGTGTGCCCGAAAAACCCATTTTCTTGCGTGTTCTTTTGGGAGCGGT‑CGTGCCCGC‑‑‑GAGGAGAGCAACCCGATTGTTGGCACCGGCTTCGCTGTTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:533881/137‑1 (MQ=25)

                                                         gCCTGATTTTCCCACAAACCTCGTTTTCTTTCGCGTTCTTGTGGGAGTGGT‑CGTGCCAGC‑‑‑GAAGGGAGCAACCCGCTTGTTGGCACCGGTTTCTCCGGTGTTCTCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:1287260/137‑1 (MQ=21)

                                                         gCCTGATTTGCCCAAAAAACTCATTGTCTTTAGCGGTCTTGGGGGAACGGG‑CGTGCCCGC‑‑‑GAAGAGAACAACCCGATCTTTGGCACCGGCTTCACCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:18855558/137‑1 (MQ=25)

                                                         gCCTGATGTTCCCGCAAAACCCATTGCCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGTTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:16896820/137‑1 (MQ=255)

                                                         gCCTGATGTGCCCACAAAACTCATTGCCTTGCTCGTTCTTGTGGGAGCGGG‑CGGGCCCGC‑‑‑GAAGAGAGCACCCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:13703277/137‑1 (MQ=37)

                                                         gCCTGAGGGGCCCGAAAAACTCATTGTCTTGCGTGGTCTTGGGGGAGCGGG‑CGGGCCCGC‑‑‑GAAGAGAGCAACCCGATTGTGGACACCGGCTTCGCTGGTGTTCGCGGGCACCCCCGCTCCCACATCGATAatgcatg                                                                               <  1:13065610/137‑1 (MQ=255)

                                                         gCCGGATGTTCCCGCAAACCTTGTTTTTTTGTGTGTTTTTGTGGGAGCGGG‑CGTTCCCGC‑‑‑GAAGAGAGCAACCCGATTGTTGGCACCGGCTTCGCTGGTGTTCGCGGACACCCCCGCTCCCACATCGATAatgcatg                                                                               <  1:15429207/137‑1 (MQ=255)

                                                         ccttGATGTTCCCGCAAAACTCCTTGTCTTGCGCGGTTTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCACCCCGATCGTTGACACCGTCTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAatgcatg                                                                               <  2:11480829/134‑1 (MQ=37)

                                                           cggATGTACCCACAAAACTCATTTTCTTGCGCGGTTTTGTGGGAGGGGG‑CGTGCCCGA‑‑‑AAAGAGAGAAACCCGATTGTGGGCACCGTTTTCGCTGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGcc                                                                             <  2:1588828/135‑1 (MQ=25)

                                                                            cTCATTGTCTTGCGCGGTCTTGTGGGAGTGGG‑CGTGCCTGC‑‑‑GAAGAGAGCAACCCGATCTTTGGCACCGGCTTCTCTGGTGTTCGCGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGcc                                                            <  2:3374003/137‑1 (MQ=255)

                                                                                ttttCTTGCGCGTTCTTTTGGGAGCGGG‑TGTGCCCGC‑‑‑GAAGAGAGCAACCCGCTTGTTGGCACCGGCTTCGCCGGTTTTCGCGGGCATGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgc                                                        <  2:12617618/134‑1 (MQ=25)

                                                                                 tGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGGGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAgcg                                                       <  2:14318186/137‑1 (MQ=255)

                                                                                      ggcgtGGTCTTGTGGGGTCGGG‑TGTGCCCGG‑‑‑GAGGAGAGAAACCCGATTGTGGGAACTGGTTTCGCTGGTGTTTGTGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTaa                                                  <  1:4123890/136‑1 (MQ=255)

                                                                                          tggTTTTGGGGGGGTGGG‑GGGGCCCGG‑‑‑GAGGAGAGCACCCCGGTTGTTGGCATTGTTTTCGCTGGTTTTCGCGGACACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAAtctc                                              <  1:489775/136‑1 (MQ=255)

                                                                                          ggttCTTGTGGGAGTGGT‑TGTCCCTGC‑‑‑GAGGAGAGCACCCCGATCGTGGGCACTGTTTTCGCCGGTGTTCGCGGACACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAAtctc                                              <  1:9881061/134‑1 (MQ=255)

                                                                                            tttttGTGGGGGGGGG‑TGTGCCGGG‑‑‑GAGGGGAGAAACCTGGTTGTTGGCACTGTTTTCGTTTTTTTTTGTGGGCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  1:9395522/134‑1 (MQ=255)

                                                                                            tttttGGGGGAGCGGG‑TGTGCCTGC‑‑‑GAAGAGAGCAACCCGATCGTGGGCACTGGTTTCGCCGGTGTTCGCGGGCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  1:10266354/134‑1 (MQ=38)

                                                                                            ttCTTGTGGGAGAGGG‑CGTGCCCGG‑‑‑GAAGAGAGCAACCCGATCGTTGACACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  2:11380488/136‑1 (MQ=255)

                                                                                            gTCTTGTGGGAGCGGG‑TGTGCCCGC‑‑‑GAAGAGAGCAACCCGGTCGTTGGCCCTGGCTTCTCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  2:8674059/137‑1 (MQ=38)

                                                                                            gTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAAGGCAACCCGATCTTTGGCACCGACTTCGCCGTTGTTCGCGGGCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  2:18503303/137‑1 (MQ=37)

                                                                                             gtttGTGGGGGTGGGTCGTGCCAGG‑‑‑GAAGAGGAAACCCCGGTTGTTGTCACTGGTTTTGCTGGTTTTCGCGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCg                                            <  1:1415047/135‑1 (MQ=255)

                                                                                               tggtgGGGGTGGG‑AGTTCCCGA‑‑‑GAAGAGAACAACCCGGTCGTTGAAACCGGCTTCGCCGGTTTTTGCGGGCAATCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGcc                                         <  1:3881780/135‑1 (MQ=255)

                                                                                                tgggggAGAGGG‑CGTGCCCGA‑‑‑GAAGAGAGCAACCCGATTGTTGGCACTGGTTTTGCTGGTGTGCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCa                                        <  1:11400137/134‑1 (MQ=37)

                                                                                                         ggg‑TGTGCCGGC‑‑‑GAAGGGAGCAACCTGATTGTTGGCACTGGTTTTGCCGTTTTTCGCGGCCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACata                               <  1:6038347/137‑1 (MQ=255)

                                                                                                          gg‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACatat                              <  2:13778004/137‑1 (MQ=255)

                                                                                                             cGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTGGACACCGGCTTCTCTGGTGTTTGCGGGCACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCg                            <  1:12729383/137‑1 (MQ=255)

                                                                                                               ttcccGC‑‑‑GAAGAGAGCAACCTGATCTTTGGCACTGTCTTCGCCGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:15312312/135‑1 (MQ=39)

                                                                                                               ttcccGA‑‑‑GAAGAAAGCACCCCGCTTTTTGACACCGTCTTTCCTGGTGTTTGTGGGCACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:6215411/135‑1 (MQ=255)

                                                                                                               tGCCTGC‑‑‑GAGGGGAGCAACCCGATCGTTGTCACTGGCTTTGCTGGTTTTCGTGGGAACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:2422520/137‑1 (MQ=255)

                                                                                                               tGCCTGA‑‑‑GAGGAGAGAAACCTGGTCGTGGGAACCGGTTTTGTGGGTGTTCGAGGGCATTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:10073780/137‑1 (MQ=255)

                                                                                                               tGCCCGG‑‑‑AAAGAGAGAAACCCGATTGTTGGCACTGGTTTCGCTGGTGTTCGCGGACACGCCCGCTCCCACATCGATAATGCATGCCCTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:16381434/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAGGAGAGCAACCTGGTCGTTGGCACTGGTTTTGCTGGTGTTCGCGGACACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:16021068/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAAGAGAGCACCCCGATCGTGGGCACCGGTTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:13096663/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAAGAGAGCAACCTGATTGTTGGCACTGGCTTTGCCGGTGTTCGCGGACATACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:5386344/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAAGAGAGCAACCCGGTTGTTGGCACCGGTTTCGCTGGTGTTTGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:7390748/137‑1 (MQ=39)

                                                                                                               tGCCCGC‑‑‑GAAGAGAGCAACCCGATTGTTGGCACCGGCTTCGCTGGTGTTCGCGGGCACACCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:12180692/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  2:16299340/137‑1 (MQ=255)

                                                                                                               tGCCCGC‑‑‑GAAGAGAACAACCCGATCGTTGGCACCGTCTTCGCCGGTGTTCGCGGACACTCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:3401801/137‑1 (MQ=255)

                                                                                                               tGCCCGA‑‑‑GAAGAGAGCAACCCGATCGTTGGAACTGTCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATACATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:18192706/137‑1 (MQ=39)

                                                                                                               tGCCCGA‑‑‑GAAGAGAACAACCCGATCGTTGGCACTGTCTTCGCAGGTGTTCGCGGACACCCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:16868519/137‑1 (MQ=255)

                                                                                                               tGCCCAT‑‑‑GAAGAGAACAACCCGATCGTGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGcc                          <  1:12526032/137‑1 (MQ=38)

                                                                                                                 cccGC‑‑‑GAAGAAAGAAACCTGGTTGTGGGCACTGGTTTCGCTGGTGTTTGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCg                        <  1:7701816/137‑1 (MQ=255)

                                                                                                                   cgc‑‑‑gAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCAAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTc                      >  1:4469175/1‑137 (MQ=255)

                                                                                                                    gc‑‑‑gAAGAGAGCAACCCGGTCGTTGGCACTGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCt                     <  2:13323727/137‑1 (MQ=255)

                                                                                                                    gc‑‑‑gAAGAGAGCAACCCGATTGTGGGCACCGGTTTCGCCGGTGTTCGTGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCt                     <  2:10136889/137‑1 (MQ=255)

                                                                                                                           agagagCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCTATAATGCATGCCTTTACCGTACACCGGCCTAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTc                 >  1:8768744/1‑137 (MQ=255)

                                                                                                                            gagagCAACCCGATCGTTGGCANCGGTTTCGCTGGTGTTCGCGGGCACGCCCGCGCCCACATCGATAATGCATGCCTTGACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCt                <  2:8504021/137‑1 (MQ=255)

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

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

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

                                                                                                                                 cAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGTGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATc           >  2:2134042/1‑137 (MQ=255)

                                                                                                                                 cAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATc           >  2:4494735/1‑137 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                       gATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACCCCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCctgctg     >  2:3158236/1‑137 (MQ=255)

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

                                                                                                                                         tCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTgctgct   <  2:5434507/137‑1 (MQ=255)

                                                                                                                                         tCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGAACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTgctgct   <  2:10015642/137‑1 (MQ=255)

                                                                                                                                          cGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTc  >  1:5571573/1‑137 (MQ=255)

                                                                                                                                          cGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTc  <  2:4699537/137‑1 (MQ=255)

                                                                                                                                          cGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTc  >  2:17640275/1‑137 (MQ=255)

CTATCGCTGGGCGAGCGGCGCAAGGATATCGAGGCGCATCTGCTGCGATCCACCCAAGCCTGATGTGCCCGCAAAACTCATTGTCTTGCGCGGTCTTGTGGGAGCGGG‑CGTGCCCGC‑‑‑GAAGAGAGCAACCCGATCGTTGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACATCGATAATGCATGCCTTTACCGTACACCGGCCGAGCGCGTAATCTCCGGCCATGCCACATATCGCCCGTCTGTTCTGGATCCTGCTGCTC  >  AE015451_phaM1/5671519‑5671789

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