BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,499,508	+C	intergenic (+85/‑79)	trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	1,499,506	1	.	C	93.0%	967.1 / 78.8	318	intergenic (+83/‑81)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base . (21/1); new base C (188/102); total (215/103)
Fisher's exact test for biased strand distribution p-value = 1.82e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GTT‑T‑CGGG‑C‑AC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGGGTGCCCGCGAAAAAAGAGCCACGGTGCCCGGCACCGGCTTTGCCAGT‑GTT‑TGCGGG‑C‑Ac‑g‑cc                                                                                                                                   >  1:16450621/1‑137 (MQ=255)

       ccATTAAGAGTCGGAATGAAGCCTTGGCACCACAGCGGGACTGAGGCTGCCAAATGGTGGGA‑GCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCGTGGCACCGGCTTTGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTccc                                                                                                                             <  2:10116964/137‑1 (MQ=18)

         gttCAGAGCCTGAACGAAGCTTTGGCACCACAGAAGGACTGGGGCTGCAAAATGTGGGGA‑GCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCTGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCcac                                                                                                                           <  1:13371314/136‑1 (MQ=17)

         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGTGTGTCCGCGAAAAAGGCGCCGCGGGGCATGGCACCGGATTTGACAGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCcaa                                                                                                                           >  2:2557943/1‑136 (MQ=25)

         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGTCTGGCACCGGCTTTGCCGGT‑GTT‑TGCGGG‑C‑AC‑GCCCGCTCCcac                                                                                                                           >  2:11918265/1‑137 (MQ=21)

              gagCCTGAACGAGGCCTTGGCACCACAGCAGGACTGAGGTTGACAAACGGGGGGG‑GTGTGTGTGCCCGCGAAACGGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTCCCACAGACt                                                                                                                      <  2:12583346/137‑1 (MQ=17)

              gagCCGGAACGAAGCCTTGGCACCACATCAGGATTGAGGTTGCCAAATGGTGGGA‑GCGGGCGTGCCCGCGAAACAGTTGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACt                                                                                                                      <  2:7601525/137‑1 (MQ=17)

                   tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCATGCCCGCGAAAAAGGCGCCGCGGTGCATTGCAACGGCTTTGTCGGT‑GTTCG‑CGGG‑T‑AC‑GCCCGCTCCCACAGACTGGTTg                                                                                                                 >  2:7349167/1‑137 (MQ=21)

                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGGGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACCAGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGAATGATTGc                                                                                                                >  2:5023688/1‑137 (MQ=18)

                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCAGGAAACAGGCGCCGCGGTGCATGGCAACGGTCTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGATCCCAAAGACTGATTga                                                                                                                >  2:16402286/1‑136 (MQ=25)

                       cGAATCCTTGACACCACAGCAGAACTGAGGCCACCAAACGGTGGGG‑GCGGACGTGCCTGCAAAACAGGCGCCGCGGTACATGGCACCGGTTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCAt                                                                                                             <  2:17391876/137‑1 (MQ=16)

                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGGAACAGGAGCCGAGGTGCCTGGCACCGGCTTTGCCGGT‑GTTCC‑CGGG‑C‑AC‑GCCCGCTCCCACAGAATGATTGCCAt                                                                                                             >  2:8047876/1‑137 (MQ=21)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGGGTGCCCACGAAACAAGCGTCGCGGTGCCTTGCACCGGCCTTTTCGGT‑GTTCG‑CGGG‑C‑AC‑GACCGCTCCCAAAAACAGATTTCCAtt                                                                                                            >  1:4041925/1‑137 (MQ=255)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGTCCACAAAACAGGCGCCACAGTGCCCTGCAACAGCCTTTCCGGT‑GTTCT‑CGGGAC‑CC‑C‑CCGCCCCCCCAGAACGATTGcccct                                                                                                            >  1:3128091/1‑134 (MQ=255)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGGCCGCGAAACAAGCACCGCGGTGCATGGCAACAGATTTGCCGGT‑GTTCT‑CGGG‑C‑AC‑GACCACTCCCACAGAATGACTGCCagt                                                                                                            >  1:17589470/1‑135 (MQ=255)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGT‑GTT‑TGCGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                            >  2:15096648/1‑137 (MQ=21)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGCACCGGCTTTGCCGGTGGTC‑G‑CGGG‑C‑AC‑GACCGCTCACAAAGAATGAGTGACAtt                                                                                                            >  2:4689714/1‑137 (MQ=18)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAACAGCGCAGAGGTGTATGGCCCCGGCCTTGCCAGTGGTC‑A‑CGGG‑C‑ACAC‑CCACTCCCCAAGACTGATTggcctt                                                                                                            >  2:12292264/1‑132 (MQ=16)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAAAAGAGCCGCGGTGCATGGCCCCGGGTTTGCCTGT‑GTTCT‑CAGG‑C‑ACAC‑CCACCCCCACAGAACGATTGCCAtt                                                                                                            >  1:17212110/1‑137 (MQ=21)

                          atctttGACACCACAACAGAACTGAGGCTACCAAAAGGTGGGA‑GCAGGCGTGCCCGTGAAACGGGCGCCGTGGGGCATGGCACCGGTTTTGCCGGT‑GTT‑CGTGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTaa                                                                                                          <  2:312703/133‑1 (MQ=12)

                            ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAAAGGAGCCGCAGGGCAGGGCAACGGATTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCCCAGACTGATTGCCATTAAga                                                                                                        >  2:3359670/1‑136 (MQ=21)

                            ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCACGAAACAAGCGCCACGGTGCCTGGCACCAGCTTTTCCAGT‑GTTCG‑CCGG‑C‑AC‑GCCCCCCCCCACAAAACGATTGCCCttcaac                                                                                                        >  2:7098666/1‑133 (MQ=255)

                              ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCgg                                                                                                      >  2:3511199/1‑137 (MQ=18)

                                 gCACCACAGCAGGACTGAGGCTGCCAAATGGTGGGA‑GTGGGCGTGCCTGCGAAACGGGCGCCGTGGTGCATGGCACCGGCTTTGCTGGT‑GTTCG‑TGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  1:16925188/137‑1 (MQ=21)

                                 gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGTGTGCCCGCGAAACAGGCACCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:1703177/137‑1 (MQ=21)

                                 gCACCACAGCAGGACTGAGGCTGCCAAA‑TGTGGGAGGTGGTCGTGCCCGCAAAACAGGCACCTCGGGACAGGACACCGGCTTTGCCGGT‑GTT‑TGCGGA‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  1:3582824/137‑1 (MQ=16)

                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCt                                                                                                  >  2:14180485/1‑137 (MQ=21)

                                         gcagGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCACCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCtc                                                                                           >  2:3700121/1‑136 (MQ=21)

                                         gcagGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                                           >  2:8330667/1‑137 (MQ=21)

                                         gcagGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCCCAAAACAAGCACCACAGTGCATGGCACCGGCTTTTCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCACCCCCACAAACTGATTGTCATTAAGCGGGACTATCGCtc                                                                                           >  2:13180200/1‑136 (MQ=12)

                                           agGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGT‑GTTCG‑CGGG‑T‑AC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                         >  2:17918428/1‑137 (MQ=18)

                                                  aGGCTGCCAAACGGTGGGA‑GCGGGGGTGCCCGCGAAACAGGCGCCGCGGTGTCTGGCACCGGATTTGCCGGT‑GTTCT‑CGGG‑C‑AC‑GCCCGCTCCCACAGAATGATTGCCACTAAACGGGCCTATCGCTACAGTGACgg                                                                                  >  2:880318/1‑137 (MQ=21)

                                                  aGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAAAAAGCGCCGCTGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑CC‑GACCGATCCCAAAGACCGCTCTCCATTAAACGGGCACATCGCTACAATGACgg                                                                                  >  2:18264003/1‑137 (MQ=9)

                                                  aGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCCCAAAACAAGCGCCGCGGTGCCTGGCACCGGCTTTTTCTGT‑GTTCT‑CGGG‑C‑AC‑GCCCGCTCCCCCAAAATGATTGCCCTTCAGCGGGCCTCTCGCTCCAATGACgg                                                                                  >  2:16583295/1‑137 (MQ=255)

                                                  aGGCTGCCAAACGGTGGGA‑GCGGGAGTGCCCGCGAAACAGGAGCCACGGTGCATGGCCCCCGCTTTTCCGGT‑TTTCG‑CCGG‑G‑AC‑GCCCGCTCCCACAGAACGGGTGCCATTCAGCAGGACTAGCGCTACCGTGAcag                                                                                  >  1:7474550/1‑135 (MQ=255)

                                                     cttccAAATGGTGGGG‑GCGGGCGTGCCTGCGAAACGGGCGCTGCGGTGAATGGCACTGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                               <  2:9298235/134‑1 (MQ=18)

                                                     cTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACCAGCGCCCCGGTGCACGGCACCCGCTTTTCCGGT‑GTTCT‑CTGG‑T‑AC‑GCCCACTCCCCCAGACCGATTGCCATTCAGCGGGCCTATCGCTACAGTGACGGacc                                                                               >  2:16304470/1‑137 (MQ=16)

                                                     cTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCAAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑TGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCCTTAAGAGGGCATATCGCTAAAGTGACGGccc                                                                               >  2:18249873/1‑137 (MQ=18)

                                                     cTGCCAAACGGTGGGA‑GCGGGCGTGCCCACGAAACAGGCACCACGGTGCATGGCACCGGATTTGCCGGT‑GTTCA‑CGGG‑C‑AA‑GACCCATCCCACAAAATGATTGCCATTAAGCGGGCCTATCGCTCCAATGACGGccc                                                                               >  2:2116738/1‑137 (MQ=255)

                                                     cTGCCAAACGGTGGGA‑GCGGGAGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCAGCTTTGCCAGT‑GTTCG‑CAGG‑C‑ACAC‑CCACTCCCACAGAATGATTGCCATTCAGCGGGTCTACCTCTACAGTGACGGccc                                                                               >  2:731397/1‑137 (MQ=255)

                                                     cTGCCAAACGGTGGGA‑GCGGGAGTGCCCGCGAAACAGGAGCCACGGTGCATGGCCACAGATTTGCCGGT‑GTTCG‑CGGG‑C‑ACAC‑CCGCTCCCAAAGAATGATTGACACTAAACGGGCCCATCGCTACAGGGACGGccc                                                                               >  2:14343223/1‑137 (MQ=16)

                                                          aaaCGGTGGGA‑GCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGTCGGT‑GTTCG‑CGGG‑C‑AA‑GCCCGCTCCCACAAAACGATTGCCATTCCACGGGCCTATCGCTCCAATGACGGCCCGgttt                                                                          >  2:12708698/1‑137 (MQ=9)

                                                           aaCGGTGGGA‑GCGGGCGTGTCCGCAAAACAGGCGCCGCGGTGCATGGCACCGGATTTGTCGGT‑GTTCG‑CGGG‑C‑CC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGGCTATCGCTACAGTGACGGCACGttttt                                                                         >  2:14991187/1‑137 (MQ=255)

                                                           aaCGGTGGGA‑GCGGGCGTGCCCGCGAAACAAGCGCCACGGTGTCTGGCACCGGCTTTTCCGGT‑GTTCT‑CGGG‑C‑ACAC‑CCGCTCCCCCCAACCGATTGCCCTTCAGCGGGTCTCTCGCTCCAATGACGGCCCCttttt                                                                         >  2:16812437/1‑137 (MQ=9)

                                                                  ggA‑GCGGGCGTACCCGAGAAACGGGCGCCGTGGTGCATGGCACCGTTTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                  <  1:13921118/137‑1 (MQ=25)

                                                                      gCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGATCCCAAAAAATGATTGCCACTAAACGGGCCCATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCg                                                               >  2:17764594/1‑137 (MQ=18)

                                                                         ggCGTTCCCGCGAAACAGTCGCCGCGGTGCGTGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTa                                                            <  2:18901299/137‑1 (MQ=25)

                                                                             tccccGCGAAAAAGGCGCCGCGGTGCATGGCACCTGCTTTGCCGGT‑GTTCG‑TGGG‑C‑AC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                                                        <  1:10522407/135‑1 (MQ=255)

                                                                                 cgcgAAACGGGCGCCGTGGTACGTGGCACCGGCTTTGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:2520232/137‑1 (MQ=25)

                                                                                 cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCTGGG‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:13772008/137‑1 (MQ=25)

                                                                                 cgcgAAACAGGCGCCGCGGTACATGACACTGGCTTTACTGGT‑GTTCG‑CGGG‑C‑AC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:18942267/137‑1 (MQ=255)

                                                                                 cgcaAAACAGGCGCCGCGGTGCAGGACACCGGTTTTGCGGGT‑GTT‑TGCGGG‑C‑AT‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:12463548/137‑1 (MQ=37)

                                                                                 cgaaaaaCAGGCGCCGCGGTGCATGGNACCGGTTTTTCTGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:18752229/133‑1 (MQ=37)

                                                                                 agaaaaaCGGGCGCCGCGGGAAATGGCACTGGTTTTGCTGGT‑GTTCG‑TGGG‑C‑AC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:6544651/133‑1 (MQ=255)

                                                                                  tcgAAACAGGCGCCGCGTGGCATTGCGCCGTCTTTGCCGGT‑GTTCT‑AGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAATCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:3911524/136‑1 (MQ=255)

                                                                                  gcgAAACAGGCGCCGGGGTGCGTGGCACCGGTTTTGCTGGG‑GTTCG‑TGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  1:8970433/137‑1 (MQ=25)

                                                                                        cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                             <  2:4310713/137‑1 (MQ=37)

                                                                                            cgctgtggTGCGTGGCACCGGTTTTGCTGGT‑GTTCT‑CGGA‑C‑AC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:13180200/131‑1 (MQ=255)

                                                                                            cgctgtggTGCATGGCACTGTTTTTGCTGGT‑GTTCG‑TGGG‑C‑AC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:712035/131‑1 (MQ=255)

                                                                                            cgccggGGTGCGTGGCACCGGTTTTTCCGGT‑GTT‑T‑GAGG‑CAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:5814654/137‑1 (MQ=255)

                                                                                            cgccgcGGTGCGGGGCACCGGCTTTTCTGGT‑GTT‑CGTGGA‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:10943708/137‑1 (MQ=255)

                                                                                            cgccgcGGCGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACTGCCCGTTTTTTGTTGCCTCGCTCACGAATCACTACTcccctc                                         <  2:13475873/137‑4 (MQ=21)

                                                                                              ctgtggTGCATGGCACCGGCTTTGCCGGG‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       <  2:6305387/133‑1 (MQ=25)

                                                                                                   gTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                  >  1:14106601/1‑137 (MQ=37)

                                                                                                        tggcaccGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                             >  1:10094774/7‑137 (MQ=37)

                                                                                                         ggcaccGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  2:5607122/6‑137 (MQ=37)

                                                                                                            accGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                      <  1:18054088/58‑1 (MQ=11)

                                                                                                            accGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  1:8411261/3‑137 (MQ=38)

                                                                                                            accGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  1:5428500/3‑137 (MQ=38)

                                                                                                            accGGCTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  1:3726759/3‑137 (MQ=38)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                      <  1:11756771/55‑1 (MQ=11)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTACCGAATCCTGGCTACCTTAGAACGCTATCAAGCAGATCTg                    >  2:10509169/1‑137 (MQ=38)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGACCGCTATCAAGCAGATCTg                    >  1:10173760/1‑137 (MQ=38)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:15361591/1‑137 (MQ=38)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:9391067/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:12975891/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:7857970/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:14954298/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:15687564/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:2188346/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:16775714/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:8565938/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:12946056/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:10986664/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:14576963/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:15535845/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:2331679/1‑137 (MQ=39)

                                                                                                                 cTTTGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:15936869/1‑137 (MQ=39)

                                                                                                                   ttGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:11065676/1‑137 (MQ=39)

                                                                                                                   ttGCCGGT‑GTTCG‑CGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTGaa                  >  2:13989900/1‑137 (MQ=39)

                                                                                                                   ttGCCGGT‑GTT‑CGCGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:11600167/1‑137 (MQ=39)

                                                                                                                            tgttcgcGGG‑C‑AC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc           <  1:17374385/131‑1 (MQ=255)

                                                                                                                            tgttcgcGGG‑C‑AC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc           <  1:16812437/131‑1 (MQ=255)

                                                                                                                            tgttcgcGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACAGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc           <  2:2897970/131‑1 (MQ=255)

                                                                                                                             gttcgcGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg          >  1:8097101/6‑137 (MQ=255)

                                                                                                                             gttcgcGGG‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg          >  1:11073181/6‑137 (MQ=255)

                                                                                                                              ttcgcGGG‑C‑AC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa         <  1:11134095/133‑1 (MQ=255)

                                                                                                                                   ggg‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTCGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:9808486/1‑137 (MQ=255)

                                                                                                                                   ggg‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:13802608/1‑137 (MQ=255)

                                                                                                                                     g‑C‑AC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc  >  1:14766411/1‑137 (MQ=255)

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGA‑GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GTT‑T‑CGGG‑C‑AC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

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