BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,819,861	G→T	5.0%	intergenic (+56/‑419)	rpoS → / → PP_1624	RNA polymerase sigma S factor/group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,819,861	0	G	T	5.0%	360.4 / 2.4	115	intergenic (+56/‑419)	rpoS/PP_1624	RNA polymerase sigma S factor/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand): ref base G (54/55); new base T (2/4); total (56/59)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 116 total aligned reads displayed.

TGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTG  >  NC_002947/1819712‑1819987

tGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAg                                                                                                                                <  2:1318659/150‑1 (MQ=255)

     gggTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATg                                                                                                                           <  1:2004171/150‑1 (MQ=255)

       gTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGgc                                                                                                                         >  1:682750/1‑150 (MQ=255)

             cAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTg                                                                                                                   >  2:1356365/1‑150 (MQ=255)

              aaaTCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                                  <  2:1555447/150‑1 (MQ=255)

                aTCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTgcgc                                                                                                                <  1:1215793/150‑1 (MQ=255)

                 tCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGcgcg                                                                                                               >  2:2247008/1‑150 (MQ=255)

                    aGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                            >  1:723417/1‑150 (MQ=255)

                    aGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                            <  2:1924510/150‑1 (MQ=255)

                      gTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgc                                                                                                          <  1:2247008/150‑1 (MQ=255)

                        ggAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgcgc                                                                                                        >  1:1458832/1‑150 (MQ=255)

                         gAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCa                                                                                                       >  2:1366942/1‑150 (MQ=255)

                           aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAAc                                                                                                     >  1:1485920/1‑150 (MQ=255)

                           aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAAc                                                                                                     <  1:336438/150‑1 (MQ=255)

                           aGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCCACCCGCTGTTCCATTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAAc                                                                                                     <  1:1687801/150‑1 (MQ=38)

                             ggCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAcccc                                                                                                   <  2:682750/150‑1 (MQ=255)

                              gCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCg                                                                                                  >  1:2078154/1‑150 (MQ=255)

                              gCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGACCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCg                                                                                                  >  1:1741238/1‑150 (MQ=255)

                               cTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGc                                                                                                 <  2:1997817/150‑1 (MQ=255)

                                tCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCt                                                                                                <  1:2096455/150‑1 (MQ=255)

                                  aaGCGCTTGCGCGAGATCCTTGGCCAGCCTGGATGCTCGGGTCAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                              <  2:1776364/150‑1 (MQ=25)

                                  aaGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                              >  2:851413/1‑150 (MQ=255)

                                    gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGt                                                                                            >  2:847733/1‑150 (MQ=255)

                                    gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGt                                                                                            >  1:1872822/1‑150 (MQ=255)

                                       cTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                                         <  1:1367961/150‑1 (MQ=255)

                                         tGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                                       <  1:705861/150‑1 (MQ=255)

                                         tGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                                       <  1:1366942/150‑1 (MQ=255)

                                               agaTCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                              <  1:577197/137‑1 (MQ=255)

                                               agaTCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                              >  2:577197/1‑137 (MQ=255)

                                                 aTCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGACCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                               <  2:1741238/150‑1 (MQ=255)

                                                  tCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                              <  2:1872822/150‑1 (MQ=255)

                                                   ccTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTg                                                                             >  2:1587354/1‑150 (MQ=255)

                                                    cTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGTTTGTGGTGGCCAAGTgg                                                                            <  2:2078154/150‑1 (MQ=255)

                                                     ttGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                        >  2:1810272/1‑137 (MQ=255)

                                                     ttGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                        <  1:1810272/137‑1 (MQ=255)

                                                     ttGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGGGAGG‑T‑CGGGGAG‑ATATCGCGAGCTCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGGGTGGGGTAAGTggg                                                                           >  1:165875/1‑150 (MQ=25)

                                                       gAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTg                                                                         >  1:2286763/1‑150 (MQ=255)

                                                        aaaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGTAACCCGCTCCGCCGCGGTGGCTAAGTGGGTgc                                                                        >  2:382539/1‑149 (MQ=37)

                                                         aaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                          >  1:285559/1‑147 (MQ=255)

                                                         aaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                          <  2:285559/147‑1 (MQ=255)

                                                         aaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                                       >  2:1498085/1‑150 (MQ=255)

                                                             aaTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCACGGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                                   >  2:317922/1‑150 (MQ=255)

                                                              aTGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                                  <  1:814365/150‑1 (MQ=255)

                                                                  tCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                                <  2:2293293/132‑1 (MQ=255)

                                                                  tCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                                >  1:2293293/1‑132 (MQ=255)

                                                                     gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           <  1:2191569/150‑1 (MQ=255)

                                                                     gTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           <  1:1498085/150‑1 (MQ=255)

                                                                      tCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCATGCCTGgcgcac                                                                                                                     <  2:403213/91‑3 (MQ=255)

                                                                      tCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCATGCCTGgcgcac                                                                                                                     >  1:403213/1‑89 (MQ=255)

                                                                        cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                        >  1:800650/1‑150 (MQ=255)

                                                                        cAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                        <  2:1706285/150‑1 (MQ=255)

                                                                         aGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           <  1:1440093/146‑1 (MQ=255)

                                                                         aGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           >  2:1440093/1‑146 (MQ=255)

                                                                           cGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                     <  1:847733/150‑1 (MQ=255)

                                                                            gAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                                    <  1:851413/150‑1 (MQ=255)

                                                                             aGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                                   <  2:40684/150‑1 (MQ=255)

                                                                                 aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCATGcc                                                                                                                             <  1:936096/72‑2 (MQ=255)

                                                                                 aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                                     <  1:2131857/112‑1 (MQ=255)

                                                                                  ttGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                              >  1:608499/1‑150 (MQ=255)

                                                                                  ttGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                              >  1:1409513/1‑150 (MQ=255)

                                                                                       ccAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTAGGTGGCTTGGAGATGAGAGTCCTCTACACACCCGGc                                         >  2:21202/1‑150 (MQ=255)

                                                                                             gctagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           <  2:921432/126‑1 (MQ=255)

                                                                                             gctagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                           >  1:921432/1‑126 (MQ=255)

                                                                                             gctagcGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                                   >  2:48310/1‑150 (MQ=255)

                                                                                                 gcgcTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                               <  1:1120536/150‑1 (MQ=255)

                                                                                                   ggtGTAACAAAAGACGCCCGGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                             <  2:723417/148‑1 (MQ=255)

                                                                                                     tGTAACAAAAAACGCCCCGGTGATC‑T‑CGTGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                                   <  2:1736727/78‑3 (MQ=21)

                                                                                                     tGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                 >  2:1990026/1‑128 (MQ=255)

                                                                                                     tGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                 <  1:1990026/128‑1 (MQ=255)

                                                                                                     tGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                                   >  1:1736727/1‑76 (MQ=25)

                                                                                                       tAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                                     >  2:891588/1‑138 (MQ=255)

                                                                                                       tAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAACTCCTCTACACACCCGGCAAgg                                     <  1:891588/138‑1 (MQ=255)

                                                                                                        acacaaaaaCGCCCC‑GAGATC‑ACCGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                               >  2:244322/5‑112 (MQ=11)

                                                                                                        acacaaaaaCGCCCC‑GAGATC‑ACCGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                               <  1:244322/108‑1 (MQ=11)

                                                                                                        aaCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATTCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGcc                        >  2:2312106/1‑150 (MQ=255)

                                                                                                         aCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                        >  2:463224/1‑117 (MQ=255)

                                                                                                         aCAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                        <  1:463224/117‑1 (MQ=255)

                                                                                                           aaaaacCGCCCCGGTGTTT‑T‑CCGTGCG‑TTTTTGTGTCCCCCATGCATGGCGCCTTGCGCGTATGCTCAACCCGCTTGGTTGTGGTGGCCAATTGGGTCACTTGGGGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa                     <  2:1035747/150‑1 (MQ=11)

                                                                                                           aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                     >  2:2206470/1‑118 (MQ=255)

                                                                                                           aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                     <  1:2206470/118‑1 (MQ=255)

                                                                                                           aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa                     <  1:382539/150‑1 (MQ=255)

                                                                                                           aaaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa                     >  2:2336117/1‑150 (MQ=255)

                                                                                                            aaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                    <  1:1356365/150‑1 (MQ=255)

                                                                                                            aaaaaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                    <  2:608499/150‑1 (MQ=255)

                                                                                                            aaaaaCGCCCC‑GAGATCAC‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                    >  1:1379764/1‑150 (MQ=11)

                                                                                                            aaaaaCGCCCC‑GAGATC‑ACCGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                               <  1:1518672/140‑1 (MQ=11)

                                                                                                            aaaaaCGCCCC‑GAGATC‑ACCGGGGCGTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                               >  2:1518672/1‑140 (MQ=11)

                                                                                                               aaCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCa                 >  1:244490/1‑150 (MQ=255)

                                                                                                                ccGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCGAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTTTTCAAGTTAATGACTCGGAGGTGACCGAGATCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa                <  2:285985/149‑1 (MQ=14)

                                                                                                                aCGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa                >  1:948238/1‑150 (MQ=255)

                                                                                                                    cccGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt            >  2:200501/1‑150 (MQ=255)

                                                                                                                      cGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                              >  2:2136634/1‑98 (MQ=255)

                                                                                                                      cGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                              <  1:2136634/98‑1 (MQ=255)

                                                                                                                        gTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTcg        >  1:2219560/1‑150 (MQ=39)

                                                                                                                            tc‑t‑cGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCggg    <  1:2312106/150‑1 (MQ=35)

                                                                                                                               t‑cGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTg  >  2:56951/1‑150 (MQ=35)

                                                                                                                                   gggCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                              <  2:550823/87‑1 (MQ=33)

                                                                                                                                   gggCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                              >  1:550823/1‑87 (MQ=33)

                                                                                                                                     gCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa                                       >  1:1492103/1‑108 (MQ=31)

                                                                                                                                     gCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa                                       <  2:1492103/108‑1 (MQ=31)

TGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑T‑CGGGGCG‑TTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTG  >  NC_002947/1819712‑1819987

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