BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,348,705	0	A	G	62.3%	53.8 / 110.6	114	intergenic (‑70/+404)	PP_3820/galU	group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand): ref base A (19/24); new base G (34/37); total (53/61)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 115 total aligned reads displayed.

TTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCA  >  NC_002947/4348588‑4348851

tttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAg                                                                                                                        <  2:130415/149‑1 (MQ=11)

  tCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                       <  3:68295/148‑1 (MQ=14)

  tCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                       <  3:262106/148‑1 (MQ=31)

   cTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                      <  7:163298/148‑1 (MQ=32)

      cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCG‑GGTAGCACAGTGGCaa                                                                                                                  <  1:137787/148‑1 (MQ=17)

       aTTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAg                                                                                                                 >  8:4900/1‑149 (MQ=17)

          aCGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                               <  6:348121/148‑1 (MQ=21)

          aCGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                               <  1:130279/148‑1 (MQ=34)

            ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                            <  3:179716/149‑1 (MQ=21)

            ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                            <  7:312887/149‑1 (MQ=34)

            ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                            <  3:174466/149‑1 (MQ=34)

             gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCTGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                           <  8:121093/149‑1 (MQ=21)

              ttACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCGACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAg                                                                                                          <  7:252530/149‑1 (MQ=21)

                aCCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                             >  2:390012/1‑144 (MQ=21)

                  cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtct                                                                                                      <  4:269808/149‑1 (MQ=21)

                  cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtct                                                                                                      <  2:248820/149‑1 (MQ=34)

                   aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                     >  1:61303/1‑149 (MQ=21)

                    aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                                                                    <  3:361732/149‑1 (MQ=21)

                    aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                                                                    <  4:51239/149‑1 (MQ=21)

                       agaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGGGGCAAGCGGACAGATCTCCGAg                                                                                                 >  4:299705/1‑149 (MQ=255)

                        gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                <  3:361468/149‑1 (MQ=21)

                        gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  4:200729/1‑149 (MQ=21)

                        gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  5:165218/1‑149 (MQ=21)

                        gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                <  7:260645/149‑1 (MQ=21)

                         aCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGTCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAggg                                                                                               >  8:295766/1‑149 (MQ=21)

                         aCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAggg                                                                                               <  1:5659/149‑1 (MQ=21)

                           ggCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGt                                                                                              >  6:166916/1‑148 (MQ=21)

                           ggCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  >  4:299048/1‑144 (MQ=34)

                           ggCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  <  3:299048/144‑1 (MQ=34)

                             cGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAt                                                                                           <  3:299705/149‑1 (MQ=34)

                               aCGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTc                                                                                         <  6:165218/149‑1 (MQ=21)

                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     <  4:141197/149‑1 (MQ=21)

                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     >  5:12726/1‑149 (MQ=21)

                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     <  6:109516/149‑1 (MQ=34)

                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     >  5:329179/1‑149 (MQ=34)

                                   cTACTCGCGGCATAACCTGTACGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCCCAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     >  8:236668/1‑149 (MQ=255)

                                         gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCtt                                                                               >  6:243295/1‑149 (MQ=21)

                                         gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCtt                                                                               >  7:350082/1‑149 (MQ=34)

                                            gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTGCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                             >  1:299549/1‑149 (MQ=255)

                                            gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                                      >  7:363413/1‑139 (MQ=21)

                                            gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                                      <  8:363413/139‑1 (MQ=21)

                                              aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                          <  3:187632/149‑1 (MQ=21)

                                               tAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                          >  1:140256/1‑148 (MQ=21)

                                                  ccTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                      <  7:346826/149‑1 (MQ=21)

                                                     gTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑Gc‑c                                                                >  8:216798/1‑149 (MQ=21)

                                                      ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑Gc                                                                  <  8:101490/147‑1 (MQ=21)

                                                      ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑Gc                                                                  <  3:108809/147‑1 (MQ=34)

                                                      ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑Gc                                                                  <  4:267264/147‑1 (MQ=34)

                                                        cGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTg                                                              <  5:278971/148‑1 (MQ=34)

                                                         ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTc                                                            >  7:55916/1‑149 (MQ=21)

                                                         ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTc                                                            <  6:227208/149‑1 (MQ=34)

                                                          gCCCTTGGTGGCGCGGTTAACCTCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGTCGCGCGGTAGCACAGTGGCTAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCg                                                           <  7:181546/149‑1 (MQ=255)

                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGa                                                         <  7:4900/149‑1 (MQ=21)

                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     <  2:224109/124‑1 (MQ=35)

                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGa                                                         >  7:101829/1‑149 (MQ=34)

                                                             cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGAt                                                        >  1:208937/1‑149 (MQ=34)

                                                                ggtggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTa                                                     >  7:97788/1‑149 (MQ=21)

                                                                ggtggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATc                                                       <  2:208937/147‑1 (MQ=34)

                                                                 gtggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGAGGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTAc                                                    <  8:350082/149‑1 (MQ=21)

                                                                  tggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACg                                                   >  7:267609/1‑149 (MQ=21)

                                                                   ggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                      <  8:80280/132‑1 (MQ=21)

                                                                   ggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                      >  7:80280/1‑132 (MQ=21)

                                                                   ggCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGt                                                  <  8:55916/149‑1 (MQ=21)

                                                                     cgcgGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTca                                                <  8:226622/149‑1 (MQ=21)

                                                                      gcgGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑Ct                                                               <  8:304358/133‑1 (MQ=21)

                                                                      gcgGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑Ct                                                               >  7:304358/1‑133 (MQ=21)

                                                                       cgGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTcaca                                              <  8:101829/149‑1 (MQ=34)

                                                                        ggTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAg                                             <  1:390012/149‑1 (MQ=21)

                                                                              ccGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTcc                                       >  3:4921/1‑149 (MQ=34)

                                                                                gtcACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGt                                     <  7:295766/147‑1 (MQ=21)

                                                                                 ccACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTg                                    >  2:144024/1‑149 (MQ=21)

                                                                                 ccACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTg                                    <  4:4921/149‑1 (MQ=34)

                                                                                  cACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGc                                   <  3:200729/149‑1 (MQ=21)

                                                                                  cACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGc                                   <  1:144024/149‑1 (MQ=21)

                                                                                    cTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                 <  8:97788/149‑1 (MQ=21)

                                                                                      tactacGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTAACAGCGTTCCGTGCAAg                                >  5:42764/1‑148 (MQ=21)

                                                                                       actacGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  5:128993/1‑86 (MQ=255)

                                                                                       actacGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                <  6:128993/86‑1 (MQ=255)

                                                                                       actacGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTa                              >  1:170679/1‑149 (MQ=34)

                                                                                           cGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGAt                                                        <  3:242657/119‑1 (MQ=25)

                                                                                           cGGCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAGCTCCCAGGGTAATTCGCGCGACCTGCGGGCTTAT‑‑GC‑CGGTCGGAGCGACGGCACAGCGGTCCGGGCAAGGATggg                          >  8:247887/1‑149 (MQ=11)

                                                                                           cGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTAtt                            <  6:329179/147‑1 (MQ=34)

                                                                                            ggCTT‑CGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTggg                         >  7:242581/1‑149 (MQ=21)

                                                                                                  cGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTTGCAAGCGTACAGATCTCCCAGGGTATTTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCtt                      <  7:185540/147‑1 (MQ=21)

                                                                                                   ggCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTgg                   >  5:266355/1‑149 (MQ=21)

                                                                                                      tGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTAATGGGCGGTGTCGGATATACGACACAGGGTTCCGTCAAAGTTTTGGGCTTTgg                   >  1:132081/1‑149 (MQ=11)

                                                                                                        aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAt               <  7:99543/148‑1 (MQ=21)

                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                                                                  <  8:332945/82‑1 (MQ=255)

                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGtt                                         >  5:153843/1‑120 (MQ=25)

                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAt                                                                                                        >  8:114098/1‑60 (MQ=255)

                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAt                                                                                                        <  7:114098/60‑1 (MQ=255)

                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg            >  5:246105/1‑149 (MQ=34)

                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg            >  1:336210/1‑149 (MQ=34)

                                                                                                            cTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCt                                                      >  7:348293/1‑105 (MQ=25)

                                                                                                               cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGGGATTTGTGCCa       >  6:239091/1‑149 (MQ=21)

                                                                                                               cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGcc        <  6:266355/148‑1 (MQ=21)

                                                                                                               cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCa       >  4:49675/1‑149 (MQ=21)

                                                                                                                  gCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCt    >  4:358806/1‑149 (MQ=21)

                                                                                                                    tCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCa  >  1:60020/1‑149 (MQ=34)

                                                                                                                      aTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCa  <  2:336210/147‑1 (MQ=34)

TTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTT‑CGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT‑‑GC‑CTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCA  >  NC_002947/4348588‑4348851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`