BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,932,640	+C	intergenic (‑107/+43)	rluA CDS (23S rRNA (746) and tRNA (32) pseudouridine synthase) ← / ← minE CDS (Cell division topological specificity factor)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	1,932,638	1	.	C	89.0%	1328.8 / 165.3	393	intergenic (‑105/+45)	rluA CDS (23S rRNA (746) and tRNA (32) pseudouridine synthase)/minE CDS (Cell division topological specificity factor)	–/–
Reads supporting (aligned to +/- strand): ref base . (39/4); new base C (200/147); total (241/152)
Fisher's exact test for biased strand distribution p-value = 9.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 414 total aligned reads displayed.

TCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑G‑C‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAT  >  AE015451_phaM1/1932502‑1932769

tCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGCAAGACCGCTCCCACCCAAGGCCTGCGCCCCCCGTGGG‑AG‑A‑GG‑G‑A‑TT‑g                                                                                                                                        >  1:7956639/1‑137 (MQ=255)

   aaGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑G‑c‑cc                                                                                                                                   <  2:16234937/137‑2 (MQ=255)

   aaGAGGATCTGGACATTCGACAGAGGCGTTGAAGTTCTCACGAAAAGAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGCAAGCCCGCTCCTACAAAGGATGTGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑G‑c‑cc                                                                                                                                   <  2:1689027/137‑2 (MQ=38)

      agGATCTGGACATTCGACAGCGGCATTGCAGTTCTCACGAAAATAGCAGTAAACGCTGGGACCGCTTTGCGGCCCTTCGCGGGCAAGCCCACTCCTACAAAGGATGTGCTTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑G‑c‑ccgcg                                                                                                                                <  2:1894579/137‑5 (MQ=38)

          tttGGACGTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCTCCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑CCC‑CGCGAAgg                                                                                                                             <  1:8128228/135‑1 (MQ=37)

           cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGACCGCTCCTACAAAGGATATGCGTCCCCTGTGGG‑AG‑C‑GG‑G‑CTTT‑T‑C‑CGCGAAggg                                                                                                                            >  2:10218721/1‑137 (MQ=38)

           cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCAGGCAAACCCGCTCCCCCAAAAGGTATGTGTCCCCCGTAGG‑AG‑C‑GG‑GTT‑TT‑C‑C‑CGCGAAggg                                                                                                                            >  1:13849651/1‑137 (MQ=255)

            tGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTTTCGGGCAAGCCCTCTCCCCCAAAAGACATGTGTCCCCTGTGGG‑AG‑CGGG‑G‑T‑TT‑C‑C‑CACGAAgggg                                                                                                                           >  1:1880822/1‑136 (MQ=255)

                cATTCGACAGCGGCATGTCAGGTCTCACGAAAATCGCAGTAAACGCTTGTGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TTAC‑C‑CGCGAATGGCTA‑Ca                                                                                                                      <  2:15037756/137‑1 (MQ=37)

                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTCCAAAGGATATGCGTCCCCTGTAGG‑GG‑A‑GG‑G‑CTTT‑C‑C‑CGCAAAGGGCTG‑CAAAGCAg                                                                                                                >  2:13680226/1‑137 (MQ=37)

                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAATCCGCTCCTCCAAAGGATATGCGTCCCCCGTGGG‑AG‑C‑GG‑G‑CTTT‑C‑C‑CACGAAGGGCTG‑CAAGGCAg                                                                                                                >  2:16445059/1‑137 (MQ=38)

                      aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCGCTCCTACAAAAGATATGCGTCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CACAAAGGGATG‑CAAGGCAg                                                                                                                >  2:4047316/1‑137 (MQ=38)

                      aCAGCGGCAGTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGGAAGCCCGCTCCTAAAAAGTATATGTTTCTCCTGCAGG‑AG‑C‑GG‑G‑CTTT‑C‑C‑CGAGAAGGGCTG‑CAAGGCAg                                                                                                                >  2:3140741/1‑137 (MQ=255)

                       cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGGTATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGc                                                                                                               <  2:6377016/137‑1 (MQ=255)

                       cAGCGGAATTGAAGGTCTCACGAAAAGAGCAGGAAACGCTGGGTCCGTTTTGCGGCCCTTCGAGGGCAAGCCCGTTCCTACAAGGGGCGTGCGTCCCTTGGGGG‑AG‑C‑GG‑T‑C‑TTAC‑C‑CGCGAAGGGCTG‑CAAGGCAGc                                                                                                               <  1:8740157/137‑1 (MQ=255)

                           ggCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCACGGGCAAGCCCGCTCCCCCAAAAGATATGCGTCCCCCGTAGG‑AA‑C‑GG‑G‑CTTT‑T‑C‑CCCAAAAGGCCT‑CAAAGCAGCCCCa                                                                                                           >  1:9734612/1‑137 (MQ=255)

                                tGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGTAAAACCACCCCCACAAAAGAGCTGTGTCCCCCATAGG‑AG‑A‑GG‑GAT‑TT‑C‑C‑CCCGAAAGGCTGCCAAGCCACcccccactg                                                                                                       >  1:8323188/1‑132 (MQ=255)

                                 gCAGGTCTCACGAAAATAGCAGTAAACTCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATCCGTCCCCTTTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGt                                                                                                     <  2:15340399/137‑1 (MQ=38)

                                  cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCACTCCCCCAAAAGAGATGCGTCCCCCGTAGG‑AGAG‑GG‑G‑C‑TT‑C‑C‑CACGAAAGGACG‑CAAAGCAACCCCcaactttc                                                                                                    >  2:13946206/1‑132 (MQ=19)

                                       ctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGACCGCTCCTCCAAAGGATATGCGTCCCCTGTGGG‑GG‑C‑GG‑GAT‑TT‑C‑C‑CGCGAAGGGCTG‑CAAGGAAGCCCCAAATTGTCGACTa                                                                                               >  2:8829179/1‑137 (MQ=37)

                                        tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGTAAGACCCCTCCCACAAAAGAGACGCGTCCCCCGTAGG‑AG‑A‑GG‑G‑CTTT‑C‑C‑CGCGAAGGGCTG‑CCAAGCAGCCCCAAACTGTCGAcccc                                                                                              >  1:2514833/1‑134 (MQ=255)

                                          aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCCCCAAAGGATGTGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑T‑TT‑GCC‑CGCGAAGGGATG‑CAAAGCAGCCCCCAATTGTCGACTAcgc                                                                                            >  2:16480224/1‑137 (MQ=25)

                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATGTGCGTCCCCTGTAGG‑GG‑C‑GG‑G‑C‑TT‑TCC‑CGCGAAGGGCTG‑TCAGGCAGCCCCAAATTGTCGACCACGCGGa                                                                                         >  2:861123/1‑137 (MQ=37)

                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAGCCCCCTCCTCCAAAAGGCTTGCGTCCCCCTTCGGAAG‑G‑GG‑A‑T‑TT‑C‑C‑CGCGAAAGGCTG‑CAAGGAAACACCAAACTGTCTAATACGCGGa                                                                                         >  2:11483467/1‑137 (MQ=255)

                                             aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCACAGGAAAGACCACCCCCCCCAAAGGCCTTTTTCCCCTGTTAG‑AG‑A‑GGTT‑T‑TT‑C‑C‑CGCAAAGGGCTT‑CCAAGAAGCCCCAAAATGTCGACTACCCGGa                                                                                         >  1:10252075/1‑137 (MQ=255)

                                                  aGCAGGAACCACTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCTCTCCTACAAAGGAGATGCGTCCCCTTTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                                                                    <  2:12057012/137‑1 (MQ=37)

                                                        aaaCGCTGGGGCCGCTTTGCGGCCTTTCGCGGGCAAGCCCGCTCCCACAAAGAATATGCTTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GTC‑CGCGAAGTGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                              <  2:4963391/137‑1 (MQ=38)

                                                        aaaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCTCTCCTACAAAGAATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                              <  2:601119/137‑1 (MQ=255)

                                                        aaaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTACTACAAAGGATATGCGTCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑G‑CTCGCGAAGGGCGG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                              <  2:894170/137‑1 (MQ=38)

                                                        aaaCGCTGGGGCCGCTTGGCGGCCCTTCGTGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTGGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                              <  2:10534525/137‑1 (MQ=38)

                                                           cgttGGGGCCGCTTTGCGGCCCTTTGTGGGAAATCCCACTCCTACAAAGGATATGCGCCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                           <  1:1957849/134‑1 (MQ=37)

                                                           cGCTGGGGCCGCTTTGCGGCCCTTCGTGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                           <  2:13387568/137‑1 (MQ=255)

                                                             ctttggCCGCTTTACGGCCCTCCGCGGGCAAGCCCTCTCCTACAAAGGATATGCGTCCCCTTTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                         <  2:11282689/133‑1 (MQ=38)

                                                             cTGGGGCCCCTTTGCGGCCCTTCGCGGACAAACCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                         <  1:9083262/137‑1 (MQ=39)

                                                                   ccGCTTGGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATGTACTTCCCCTTTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                   <  1:15848951/137‑1 (MQ=37)

                                                                     gttttGCGGCCCTTCTTGGGCAAGCCCGCTCCTACAAAGGATATTCGTCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                 <  1:2290447/135‑1 (MQ=38)

                                                                       tttGCGGCCCTTCGCGGGCAAGCCCACCCCCCCAAAGGAGATGCGTCCCCCGTAGG‑AG‑C‑GG‑GAT‑TT‑C‑C‑CACAGAAGGCTG‑CAAAGCAACACCAAACCGTCGAACACGCGGATCAACGAGCCGGGAAGGTGATGGTCa                                                               >  2:12863957/1‑137 (MQ=255)

                                                                        ttGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                              <  1:15398221/137‑1 (MQ=255)

                                                                        ttGCGGCCCTTCGCGGGCAAGCCCACTCCTACAAAGGATATGCGTCCCCTTTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                              <  2:3666380/137‑1 (MQ=255)

                                                                         tGCGGCCCTTCGCGGGCAAGACCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑CTTT‑C‑C‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACACGGATCAGCGATCCGGCAGGGGGGTGTTCAGc                                                             >  2:699365/1‑137 (MQ=37)

                                                                               ccTTCGCGGGCAAGCCCGCTCCTACAAAGGAGATGAGTCCCATGTAGG‑AA‑C‑GG‑G‑C‑CT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTAAGTGGATCAGCGATCCGGCCGGGGGATGTTCCGCTCCAGa                                                       >  2:7299966/1‑137 (MQ=37)

                                                                                 ttCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAt                                                     <  2:1442769/137‑1 (MQ=255)

                                                                                     cgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                 <  1:7153122/137‑1 (MQ=255)

                                                                                     cgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                 <  2:10850483/137‑1 (MQ=255)

                                                                                      gggCAACCCCGCTCCTACAAAGGATATGCTTCCCCTGTGGG‑AG‑T‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                <  1:100500/137‑1 (MQ=38)

                                                                                         aaaGCCCGCTCCTACAAAGGATATGCTTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagcagc                                             <  1:10889378/136‑1 (MQ=255)

                                                                                          aaGCCCGCTCCTACAAAGTATATGGGTCCCCTGTAGG‑AG‑C‑GG‑A‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAGGGCAGCCCCAAATTGTCTACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAACTCCAGAATCTAGCAGCt                                            <  2:11153555/137‑1 (MQ=255)

                                                                                           aGCCCGCTCCTACAAAGGATATGCGTCCCTTGTAGG‑AG‑C‑GG‑G‑T‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTg                                           <  1:10462737/137‑1 (MQ=255)

                                                                                                 cTCCGACAAAGGATATCCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTgg                                     <  1:4323154/137‑1 (MQ=255)

                                                                                                     tACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                   <  2:11855406/135‑1 (MQ=255)

                                                                                                       cAAAGGATATGCGTCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                               <  1:11892871/137‑1 (MQ=255)

                                                                                                        aaaGGGTATGTGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                              <  2:539560/137‑1 (MQ=255)

                                                                                                        aaaGGAGATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                              <  1:16877651/137‑1 (MQ=255)

                                                                                                         aaGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  2:6873705/1‑137 (MQ=255)

                                                                                                         aaGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  1:11180033/1‑137 (MQ=255)

                                                                                                         aaGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                             >  2:4343288/1‑137 (MQ=255)

                                                                                                          aGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                            <  2:13137684/137‑1 (MQ=255)

                                                                                                          aGGATATGCGTCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                            <  1:16411249/137‑1 (MQ=255)

                                                                                                            gATATGCGTCCCTTTGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                          <  1:5522640/137‑1 (MQ=39)

                                                                                                            gATATGCGTCCCCTGTAGG‑AG‑C‑GT‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                          <  2:14515358/137‑1 (MQ=255)

                                                                                                               gtGCGTCCCCTGTAGG‑AG‑C‑GG‑A‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGTATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTGTCCAGCTCGat                       <  1:11780192/136‑1 (MQ=39)

                                                                                                               atGCNTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCANAATCNAGCAGCTGCCCTGGTTTTCCAGCTCGat                       >  2:7087938/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGt                                                                                                     <  2:4576137/57‑1 (MQ=18)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:12792841/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:16827844/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:417122/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:1155776/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:15144995/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:5894161/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:14883927/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  2:7698194/1‑137 (MQ=255)

                                                                                                                 gCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                     >  1:12368677/1‑137 (MQ=255)

                                                                                                                  cGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                    >  2:15955159/1‑137 (MQ=255)

                                                                                                                  cGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                    >  2:6126331/1‑137 (MQ=255)

                                                                                                                  cGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                    <  1:3646524/137‑1 (MQ=255)

                                                                                                                  cGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                    >  2:5483171/1‑137 (MQ=255)

                                                                                                                   ttCCCCTGTAGG‑AG‑T‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                   <  1:52406/136‑1 (MQ=255)

                                                                                                                   gTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAGTTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                   >  1:14193324/1‑137 (MQ=255)

                                                                                                                    tCCCCTGGAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                  <  1:14938171/137‑1 (MQ=255)

                                                                                                                       ccTGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat               >  1:4651780/1‑137 (MQ=255)

                                                                                                                       ccTGTAGG‑AG‑C‑GG‑G‑C‑TT‑ACC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat               <  1:1287846/137‑1 (MQ=255)

                                                                                                                         tGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:2707944/1‑137 (MQ=255)

                                                                                                                         tGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:2092131/1‑137 (MQ=255)

                                                                                                                         tGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:8529222/1‑137 (MQ=255)

                                                                                                                         tGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  1:5409369/1‑137 (MQ=255)

                                                                                                                         tGTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCAAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca             >  2:8103975/1‑137 (MQ=255)

                                                                                                                          gTAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat            >  2:5191927/1‑137 (MQ=255)

                                                                                                                           tAGG‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc           >  1:16803805/1‑137 (MQ=255)

                                                                                                                             gg‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt         >  1:6065595/1‑137 (MQ=255)

                                                                                                                              g‑AG‑C‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt        >  2:1239/1‑137 (MQ=255)

                                                                                                                                aG‑C‑GG‑G‑C‑TT‑ACC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg       <  1:13582798/137‑1 (MQ=255)

                                                                                                                                   c‑GG‑G‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACGACGCGGATCAGCGATCCGGCAGGGTGATGTGCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGcc     <  1:4006454/137‑1 (MQ=255)

                                                                                                                                     gg‑g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCg    >  1:5071392/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:2405960/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:15714016/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:1557652/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:12515538/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  <  1:3401449/137‑1 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  1:2927247/1‑137 (MQ=255)

                                                                                                                                        g‑C‑TT‑GCC‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt  >  2:9333246/1‑137 (MQ=255)

TCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGG‑AG‑C‑GG‑G‑C‑TT‑G‑C‑CGCGAAGGGCTG‑CAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAT  >  AE015451_phaM1/1932502‑1932769

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