BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,919,611	C→T	5.4%	P701L (CCG→CTG)	yebT →	MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,919,611	0	C	T	5.4%	388.1 / 7.4	111	P701L (CCG→CTG)	yebT	MCE domain protein
Reads supporting (aligned to +/- strand): ref base C (55/50); new base T (3/3); total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 111 total aligned reads displayed.

TCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAA  >  NC_000913/1919468‑1919754

tCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAt                                                                                                                                           >  2:2612657/1‑150 (MQ=255)

  gCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                                         >  2:1216174/1‑150 (MQ=255)

  gCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                                         >  2:2036502/1‑150 (MQ=255)

  gCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                                         <  2:2574042/150‑1 (MQ=255)

  gCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                                         >  2:656681/1‑150 (MQ=255)

   cTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTc                                                                                                                                        <  2:277169/150‑1 (MQ=255)

     aTCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCgg                                                                                                                                      <  1:1144858/150‑1 (MQ=255)

     aTCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCgg                                                                                                                                      <  2:2293672/150‑1 (MQ=255)

       ctcGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCa                                                                                                                                    >  1:2420051/1‑150 (MQ=255)

             aTTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGc                                                                                                                              >  2:734345/1‑150 (MQ=255)

                       ggCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATc                                                                                                                    <  2:2666386/150‑1 (MQ=255)

                       ggCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATc                                                                                                                    <  2:847991/150‑1 (MQ=255)

                        gCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCTCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCt                                                                                                                   <  2:1765923/150‑1 (MQ=255)

                        gCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCt                                                                                                                   <  2:109431/150‑1 (MQ=255)

                        gCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCt                                                                                                                   <  2:1112311/150‑1 (MQ=255)

                        gCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCt                                                                                                                   <  2:2239322/150‑1 (MQ=255)

                           aaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGa                                                                                                                <  1:660113/150‑1 (MQ=255)

                               cAGACGCTTGAAATCATTACCGCAAGAAATGAATTAAATGCAAAGGCTGAACTTTACCCGGAATATGTCCAGACCTTTGATCGCGGTGGTACGCTATTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGAGTTGATCATCTTGATaat                                                                                                            >  2:1869158/1‑148 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  1:601423/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  2:779388/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  2:1498756/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  1:137068/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  2:453631/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  2:796048/1‑150 (MQ=255)

                               cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTAGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  1:413450/1‑150 (MQ=255)

                               cAGACGCAGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                                            >  2:2376222/1‑150 (MQ=255)

                                      tGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcctcc                                                                                                     >  1:1359041/1‑150 (MQ=255)

                                         aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGc                                                                                                  >  2:1361817/1‑150 (MQ=255)

                                         aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGc                                                                                                  >  2:755374/1‑150 (MQ=255)

                                         aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGc                                                                                                  >  2:823528/1‑150 (MQ=255)

                                         aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGc                                                                                                  >  2:895940/1‑150 (MQ=255)

                                            tGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGt                                                                                               >  1:270642/1‑150 (MQ=255)

                                            tGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGt                                                                                               <  1:948026/150‑1 (MQ=255)

                                                 aCCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACTGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATc                                                                                          >  1:243811/1‑150 (MQ=255)

                                                 aCCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACTGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATc                                                                                          >  1:243774/1‑150 (MQ=255)

                                                 aCCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACTGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATc                                                                                          >  1:755250/1‑150 (MQ=255)

                                                  ccGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCa                                                                                         >  1:1776028/1‑150 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  2:1429976/150‑1 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  2:1379349/150‑1 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  2:896852/150‑1 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  2:802402/150‑1 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  1:1926131/150‑1 (MQ=255)

                                                    gcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAAc                                                                                       <  1:1740214/150‑1 (MQ=255)

                                                          aaTGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGaa                                                                                 <  2:1972263/150‑1 (MQ=255)

                                                          aaTGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGaa                                                                                 <  2:193407/150‑1 (MQ=255)

                                                          aaTGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGaa                                                                                 <  2:2161776/150‑1 (MQ=255)

                                                             gAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCa                                                                              <  1:598059/150‑1 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:1572470/1‑150 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:1572473/1‑150 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:1572479/1‑150 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:690060/1‑150 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:690096/1‑150 (MQ=255)

                                                                 tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGcc                                                                          >  1:122429/1‑150 (MQ=255)

                                                                  aCAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCg                                                                         <  2:2122022/150‑1 (MQ=255)

                                                                  aCAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCg                                                                         <  2:2420051/150‑1 (MQ=255)

                                                                   cAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCgg                                                                        >  2:2604707/1‑150 (MQ=255)

                                                                       cAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCa                                                                    >  2:1689618/1‑150 (MQ=255)

                                                                        aaaGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCaa                                                                   >  1:1842411/1‑150 (MQ=255)

                                                                         aaGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAAt                                                                  >  2:2549479/1‑150 (MQ=255)

                                                                         aaGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAAt                                                                  >  2:51632/1‑150 (MQ=255)

                                                                            gCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCt                                                               <  1:1272673/150‑1 (MQ=255)

                                                                            gCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCt                                                               <  1:393944/150‑1 (MQ=255)

                                                                             cGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACTTCGAACCAGGCCGGGGCAATCCTc                                                              <  1:1177879/150‑1 (MQ=255)

                                                                             cGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTc                                                              <  1:588675/150‑1 (MQ=255)

                                                                                  cTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTcgccgc                                                         <  1:2464258/150‑1 (MQ=255)

                                                                                       cccGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACtt                                                    <  1:2536353/150‑1 (MQ=255)

                                                                                       cccGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACtt                                                    <  1:1914841/150‑1 (MQ=255)

                                                                                       cccGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACtt                                                    <  1:1914822/150‑1 (MQ=255)

                                                                                        ccGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACttt                                                   >  2:410276/1‑150 (MQ=255)

                                                                                                   ccAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACTGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAtga                                        <  2:755250/150‑1 (MQ=255)

                                                                                                   ccAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACTGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAaga                                        <  2:243774/150‑1 (MQ=255)

                                                                                                    cAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAagag                                       <  2:361880/150‑1 (MQ=255)

                                                                                                          tttGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGccacc                                 <  1:823528/150‑1 (MQ=255)

                                                                                                          tttGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGccacc                                 <  1:895940/150‑1 (MQ=255)

                                                                                                          tttGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGccacc                                 <  1:1361817/150‑1 (MQ=255)

                                                                                                              cTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTa                             <  1:2535398/150‑1 (MQ=255)

                                                                                                              cTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTa                             <  1:1216784/150‑1 (MQ=255)

                                                                                                              cTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTa                             <  1:2535407/150‑1 (MQ=255)

                                                                                                               tCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTAc                            >  2:1088837/1‑150 (MQ=255)

                                                                                                                    gtggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGAtt                       >  2:2532107/1‑150 (MQ=255)

                                                                                                                    gtggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGAtt                       >  2:2629793/1‑150 (MQ=255)

                                                                                                                    gtggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAACCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGAtt                       >  1:1886372/1‑150 (MQ=255)

                                                                                                                       gttcgcgCTTCTGAGTGGTCACAACGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCATTCCTCGCCGAGACTTTGAATTACAAGAGGCCACCATTACTGATTcgc                    <  2:251726/147‑1 (MQ=255)

                                                                                                                       gtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTcgc                    <  2:2159144/150‑1 (MQ=255)

                                                                                                                         aCGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGt                  >  1:1866339/1‑150 (MQ=255)

                                                                                                                         aCGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTGTGAATTACAAGAGGCCACCATTACTGATTCGCGt                  >  1:2001161/1‑150 (MQ=255)

                                                                                                                          cgcgCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGtt                 >  1:470597/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  1:1316093/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  1:719450/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  2:578665/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  2:541166/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  1:169883/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  1:2523649/1‑150 (MQ=255)

                                                                                                                           gcgcTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTa                >  1:1888318/1‑150 (MQ=255)

                                                                                                                              cTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCt             <  1:2541549/150‑1 (MQ=255)

                                                                                                                              cTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCt             <  2:2167956/150‑1 (MQ=255)

                                                                                                                              cTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCt             <  1:1178989/150‑1 (MQ=255)

                                                                                                                                      tGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCt     <  1:2705728/150‑1 (MQ=255)

                                                                                                                                         tCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTaa  <  2:2610996/150‑1 (MQ=255)

                                                                                                                                         tCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTaa  <  2:1806747/150‑1 (MQ=255)

TCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATTCGCGTTACCTGGATGGCTTAA  >  NC_000913/1919468‑1919754

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