BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,460	C→T	8.7%	intergenic (‑237/+10)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,460	0	C	T	8.7%	312.1 / 16.2	115	intergenic (‑237/+10)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (58/46); new base T (5/5); total (64/51)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 118 total aligned reads displayed.

GCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGC  >  NC_002947/5743311‑5743605

gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGTCCGCGAAc                                                                                                                                                    >  1:398230/1‑150 (MQ=255)

        aCGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGcctcac                                                                                                                                            >  1:1757508/1‑150 (MQ=255)

             acaATGGGTGGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                                                                                       >  2:2294155/1‑150 (MQ=255)

              caATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTACAGGTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                                      <  1:30175/150‑1 (MQ=255)

                              cAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAATCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAt                                                                                                                      <  2:2142093/150‑1 (MQ=255)

                               ccgagatctggaaggcattggattgaggtgtgtccattcgcgggcgcactCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGcctcacct                                                                                                                                          <  1:2545282‑M1/79‑1 (MQ=255)

                               ccgagatctggaaggcattggattgaggtgtgtccattcgcgggcgcactCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGcctcacct                                                                                                                                          >  2:2545282‑M1/51‑129 (MQ=255)

                               aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                                                                                 <  1:2526666/138‑1 (MQ=255)

                               aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                                                                                 >  2:2526666/1‑138 (MQ=255)

                                gCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                                                                    >  2:1527275/1‑150 (MQ=255)

                                  cTTCGCCAGGCAGGGGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                                                                  <  2:475962/150‑1 (MQ=255)

                                      gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAAc                                                                                                                                                    >  2:790591/1‑112 (MQ=255)

                                      gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAAc                                                                                                                                                    <  1:790591/112‑1 (MQ=255)

                                       ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                                                             >  1:1031908/1‑150 (MQ=255)

                                                tggtCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCCTGCTCCGAAATACCG‑G‑TCGCTAGCTACCTTACCTGAATCCAATTCCTTCCAGATCTCGGTGGCATATTcgcg                                                                                                    >  1:1594402/1‑150 (MQ=255)

                                                  gtCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                                                                    >  2:1742672/1‑132 (MQ=255)

                                                  gtCGGGTTGTGTGTCCATTCGCGGCCGCGCTCCCTCCCACAGGTACAGCGCAAGCCTTGAACCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGaa                                                                                                  <  2:1387240/150‑1 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             <  1:1488175/150‑1 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  2:2226147/1‑150 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  2:1432876/1‑150 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  1:2479059/1‑150 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTAACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  1:2688795/1‑150 (MQ=255)

                                                       gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑AGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  2:1901665/1‑150 (MQ=255)

                                                        cagtctgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAAAGGCCTAACCTCAAGCCAATGCCTTCCAGGTCTCGTTGGCATTGTCGCGAATGGTGc                                                                                            >  1:2410746/6‑150 (MQ=37)

                                                         tgtgtgTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                                                                           <  1:2294155/150‑1 (MQ=255)

                                                          aggtgaggccgttcgcgggcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGTGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCgg                                                                                                                 >  2:2545282‑M2/80‑129 (MQ=255)

                                                          aggtgaggccgttcgcgggcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGTGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCgg                                                                                                                 <  1:2545282‑M2/50‑1 (MQ=255)

                                                                  aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                                                               <  1:1524530/121‑1 (MQ=255)

                                                                  aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                                                               >  2:1524530/1‑121 (MQ=255)

                                                                          gcgcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCacc                                                                          >  2:775946/1‑148 (MQ=255)

                                                                          gcgcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                                                                          >  2:1827988/1‑150 (MQ=255)

                                                                          gcgcgcTCGCCCCCACAGGTACAGCGCAAGCCTGAAAGCAAGCTAGACA‑CCTG‑TGGGCTCGGGCGTGCTCGCGAACGGCCTCACCTCAATCCAATGCTTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                                                                          <  1:1105970/150‑1 (MQ=255)

                                                                            gcggggGCTCCCACAGGCACAGCGCAAGTGTTGTAGCTTGCGGGACA‑CTTG‑TGGGCGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGccc                                                                        >  2:686294‑M1/7‑150 (MQ=255)

                                                                              gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGcctcac                                                                                                                                            <  2:1732422/80‑1 (MQ=255)

                                                                              gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGcctcac                                                                                                                                            >  1:1732422/1‑80 (MQ=255)

                                                                              gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑At                                                                     >  1:986262/1‑150 (MQ=255)

                                                                                tcgctcCCACAGGTACAGCGCACGCGTTGAAGCTTTCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATcc                                                                   <  2:2688795/150‑1 (MQ=255)

                                                                                tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATGCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATcc                                                                   >  1:1697161/1‑150 (MQ=255)

                                                                                cCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCaa                                                                                                                                  <  1:1883249‑M1/87‑1 (MQ=255)

                                                                                cCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCaa                                                                                                                                  >  2:1883249‑M1/2‑88 (MQ=255)

                                                                                 cccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                                      >  1:2854660‑M2/57‑85 (MQ=255)

                                                                                 cccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                                      <  2:2854660‑M2/29‑1 (MQ=255)

                                                                                      tcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                                                             >  1:1663821‑M2/52‑105 (MQ=255)

                                                                                      tcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                                                             <  2:1663821‑M2/54‑1 (MQ=255)

                                                                                      CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATccgtgccg                                                             <  1:2106932‑M1/150‑1 (MQ=255)

                                                                                      CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGAGCGCGCCCGCGAATGGcctcac                                                                                                                                            >  1:2780454‑M1/1‑72 (MQ=255)

                                                                                      CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGAGCGCGCCCGCGAATGGcctcac                                                                                                                                            <  2:2780454‑M1/72‑1 (MQ=255)

                                                                                       CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGATCGGCCTCACCTCAATCCAATACCTTCCAGATCTCGGTGGCATATTCGCGCATGGTGCGGTCCGGCGAGAACCAGCCC‑ATccgtgccgt                                                            <  1:549189‑M1/150‑1 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACgg                                                                                                                                                  <  2:858124‑M2/17‑1 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACgg                                                                                                                                                  >  1:858124‑M2/50‑66 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             <  2:1651507‑M2/70‑1 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                             >  1:1651507‑M2/50‑119 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGTAAGGAAAGCGGGCCGACGAGAACCAGCCC‑ATccgtgccg                                                             >  2:2547727‑M2/50‑150 (MQ=255)

                                                                                        ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATccgtgccg                                                             >  2:1101647‑M2/50‑150 (MQ=255)

                                                                                          acaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATccgtgccgtg                                                           >  1:1168619‑M2/48‑150 (MQ=255)

                                                                                          AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGtt                                                         >  2:886634‑M1/1‑150 (MQ=255)

                                                                                           caggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGaa                                                                                                  >  1:1398614‑M2/47‑111 (MQ=255)

                                                                                           caggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGaa                                                                                                  <  2:1398614‑M2/65‑1 (MQ=255)

                                                                                             TACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCaa                                                      >  1:1432876‑M1/1‑150 (MQ=255)

                                                                                             TACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCaa                                                      <  1:2533337‑M1/150‑1 (MQ=255)

                                                                                               cAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGACTCACCTCAATCCAAGACCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATa                                                    >  1:1456428/1‑150 (MQ=255)

                                                                                                 tctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATa                                                    <  1:1101647‑M2/110‑1 (MQ=255)

                                                                                                 gcgcAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATAcc                                                  >  1:2502632/1‑150 (MQ=255)

                                                                                                   tcgcaagcttcaaggcttgcgagacacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTcc                                                                                                                          <  1:2142629‑M2/41‑1 (MQ=255)

                                                                                                     aaGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCTGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCt                                              >  2:2819860/1‑150 (MQ=255)

                                                                                                     aaGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGAACAACGATATCTTGCCC‑ATCCGTGCCGTGTTCAATACCGCCt                                              >  2:245467/1‑150 (MQ=255)

                                                                                                      caagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAAt                                                                                                                                 >  1:946005‑M2/36‑69 (MQ=255)

                                                                                                         cTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCCACACCGCCTTgcg                                          >  1:2773573/1‑150 (MQ=255)

                                                                                                         cTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTgcg                                          >  1:1060056/1‑150 (MQ=255)

                                                                                                         cTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTgcg                                          >  2:46943/1‑150 (MQ=255)

                                                                                                         cTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTgcc                                          >  2:1408003/1‑149 (MQ=255)

                                                                                                          ttGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                                                    >  1:2438019/1‑108 (MQ=255)

                                                                                                            gAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTc                                                        >  1:804693/1‑133 (MQ=255)

                                                                                                            gAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTc                                                        <  2:804693/133‑1 (MQ=255)

                                                                                                            gAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGcca                                       >  1:537205/1‑150 (MQ=255)

                                                                                                             aaGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                      <  1:1527275/150‑1 (MQ=255)

                                                                                                               gCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                                                    <  1:2279430/103‑1 (MQ=255)

                                                                                                               gCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                                                    >  2:2279430/1‑103 (MQ=255)

                                                                                                                 ttGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACt                                  >  2:1232237/1‑150 (MQ=255)

                                                                                                                     gagaCA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                                                                >  1:1897985/1‑69 (MQ=255)

                                                                                                                     gagaCA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                                                                <  2:1897985/69‑1 (MQ=255)

                                                                                                                     gagaCA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTg                              <  2:32529/150‑1 (MQ=255)

                                                                                                                     gagaCA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATACGCGAATGGTGCGGTCAGACGCGAATTAGCCG‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTg                              >  1:690529/1‑150 (MQ=255)

                                                                                                                       gaCA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGGCGGA‑GTGGGACAGCCCAAACCGCTCCGTGTTCAATACCGCCTTGCGCCCCCACTCCTGCg                            <  2:2410746/150‑1 (MQ=255)

                                                                                                                         tgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATccgtgccgt                                                            <  1:2623062‑M2/102‑1 (MQ=255)

                                                                                                                         ca‑cCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgt                          <  2:2502632/150‑1 (MQ=255)

                                                                                                                            acctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtg                         >  2:1903028‑M2/14‑150 (MQ=255)

                                                                                                                              ctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtg                       <  1:775946‑M2/139‑1 (MQ=255)

                                                                                                                              ctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtg                       >  2:2977647‑M2/12‑150 (MQ=255)

                                                                                                                               tgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtgt                      >  1:1104033‑M2/11‑150 (MQ=255)

                                                                                                                                     agcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCcaca                <  2:537205‑M2/146‑1 (MQ=255)

                                                                                                                                        gGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                                                     <  1:563925‑M2/62‑1 (MQ=255)

                                                                                                                                        gGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                                                     >  2:563925‑M2/2‑63 (MQ=255)

                                                                                                                                         ggCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATc            <  2:1104033/150‑1 (MQ=255)

                                                                                                                                            gTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTc         <  1:46943/150‑1 (MQ=255)

                                                                                                                                            gTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTc         <  2:1697161/150‑1 (MQ=255)

                                                                                                                                            gTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCACACC‑ACCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTc         >  1:1898551/1‑150 (MQ=255)

                                                                                                                                               cccGCGAACGGCCCCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGAc      >  2:1847683/1‑150 (MQ=255)

                                                                                                                                                 cgcgAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCTCCACTCCTGCGGTGTGTGCCACAGATCCTCGAccc    <  1:245467/150‑1 (MQ=255)

                                                                                                                                                   cgAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACccgc  <  2:1031908/150‑1 (MQ=255)

GCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACA‑CCTG‑TGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCC‑ATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGACCCGC  >  NC_002947/5743311‑5743605

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