BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,714,201	C→T	5.6%	H98Y (CAC→TAC)	PP_5012 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,714,201	0	C	T	5.6%	611.1 / 16.6	196	H98Y (CAC→TAC)	PP_5012	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (102/83); new base T (6/5); total (108/88)
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.71e-01

Only 100 of 199 total aligned reads displayed.

TCATCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGACTCGCCAA  >  NC_002947/5714052‑5714350

  aTCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTTGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCAc                                                                                                                                                     >  2:487356/1‑150 (MQ=255)

  aTCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCAc                                                                                                                                                     <  1:2938657/150‑1 (MQ=255)

   tCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACg                                                                                                                                                    >  2:1542730/1‑150 (MQ=255)

   tCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACg                                                                                                                                                    >  2:533489/1‑150 (MQ=255)

     cTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGc                                                                                                                                                  <  2:644757/150‑1 (MQ=255)

       gCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGa                                                                                                                                                >  1:4323118/1‑150 (MQ=255)

              gAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGcc                                                                                                                                         <  1:4666304/150‑1 (MQ=255)

                   gggCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCTCATCTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGtgca                                                                                                                                    <  2:2869446/150‑4 (MQ=255)

                    ggCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAgcg                                                                                                                                           <  2:1200632/142‑1 (MQ=255)

                    ggCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTg                                                                                                                                   >  1:3701463/1‑150 (MQ=255)

                    ggCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTg                                                                                                                                   >  1:3393630/1‑150 (MQ=255)

                      ccTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCt                                                                                                                                 <  2:1854563/150‑1 (MQ=255)

                       cTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTc                                                                                                                                <  1:541003/150‑1 (MQ=255)

                          aTGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGAc                                                                                                                             <  1:1733714/150‑1 (MQ=255)

                          aTGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGAc                                                                                                                             <  1:1733859/150‑1 (MQ=255)

                          aTGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGAc                                                                                                                             <  2:3086815/150‑1 (MQ=255)

                             cTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTg                                                                                                                          >  2:4297324/1‑150 (MQ=255)

                               cgccgcCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCAGTACGGCGCTGTCGCCGTGCAGCTCGACCTgcg                                                                                                                        <  2:1631443/150‑3 (MQ=38)

                               cgccgcCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCAGTACGGCGCTGTCGCCGTGCAGCTCGACCTgcg                                                                                                                        <  2:1531100/150‑3 (MQ=38)

                                 ccgccCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCg                                                                                                                      >  1:1629007/1‑150 (MQ=255)

                                 ccgccCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCg                                                                                                                      <  2:3663897/150‑1 (MQ=255)

                                 ccgccCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCg                                                                                                                      <  2:1282160/150‑1 (MQ=255)

                                 ccgccCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCACAGGTCGAGCAGTACGGCGCTGTCGCCGTGCAGCTCGACCTgcggg                                                                                                                      <  1:1328374/150‑5 (MQ=37)

                                   gccCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGgt                                                                                                                    <  2:2512855/150‑1 (MQ=255)

                                    cccATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGgtg                                                                                                                   <  1:1084631/150‑1 (MQ=255)

                                     ccATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGgtgt                                                                                                                  <  2:4666076/150‑1 (MQ=255)

                                      cATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGtgtg                                                                                                                 >  2:2192078/1‑150 (MQ=255)

                                       aTTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGtgc                                                                                                                >  2:2021288/1‑150 (MQ=255)

                                          ggCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGtgctgc                                                                                                             <  1:2192078/150‑1 (MQ=255)

                                            cAGGGCCAGCTCGACGGCCGCCTGCGCGCCCCCGCCGCTACCCTGGCGCAACTGGCCATGGTCAAGGACAAGACCGCTGTGCCGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACCGCTCGACCTGTTCGGTGTCCTGCAg                                                                                                           <  2:1568571/150‑1 (MQ=255)

                                             aGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGgtg                                                                                                                   >  2:2163585/1‑141 (MQ=255)

                                               ggCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGAc                                                                                                        <  1:5209776/150‑1 (MQ=255)

                                                    ggTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAACCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGa                                                                                                   >  2:3003291/1‑150 (MQ=255)

                                                     gTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAg                                                                                                  <  2:1817987/150‑1 (MQ=255)

                                                        gACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTg                                                                                               >  1:348934/1‑150 (MQ=255)

                                                         aCTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTgg                                                                                              >  1:5354769/1‑150 (MQ=255)

                                                             ccgccTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTg                                                                                          <  2:2480537/150‑1 (MQ=255)

                                                               gccTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTggg                                                                                        >  1:4534964/1‑150 (MQ=255)

                                                                ccTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGa                                                                                       <  2:5138926/150‑1 (MQ=255)

                                                                ccTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGa                                                                                       <  2:2049462/150‑1 (MQ=255)

                                                                ccTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGa                                                                                       >  1:3622739/1‑150 (MQ=255)

                                                                   gcgcgcCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCa                                                                                    >  2:1515038/1‑150 (MQ=255)

                                                                         ccccGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGGGTGCTGCAGGACCGGGAGCTGGACTGGGAGCACGAGCt                                                                              >  1:2760250/1‑150 (MQ=255)

                                                                                    gCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGGGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGc                                                                   <  2:171112/150‑1 (MQ=255)

                                                                                     ccTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCt                                                                  >  1:1416637/1‑150 (MQ=255)

                                                                                      cTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCtt                                                                 >  1:981299/1‑150 (MQ=255)

                                                                                       tGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTg                                                                >  2:60849/1‑150 (MQ=255)

                                                                                        ggggcAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACTGCGCCGTACTGCTCGTCGTGTTCGGTGTGCTGCAGGACGTGGAGCTGGAGTGGGAGCTCGAGCTGCAACGCTGGCTTgg                                                               >  2:2579959/4‑150 (MQ=255)

                                                                                        ggCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTgg                                                               >  2:427814/1‑150 (MQ=255)

                                                                                         gcgcAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGGTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGc                                                              <  2:4188843/150‑1 (MQ=255)

                                                                                              aCTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGt                                                         >  1:781583/1‑150 (MQ=255)

                                                                                               cTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTc                                                        <  1:5300296/150‑1 (MQ=255)

                                                                                                    cTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCAc                                                   >  2:3217799/1‑150 (MQ=255)

                                                                                                    cTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCAc                                                   >  1:1639934/1‑150 (MQ=255)

                                                                                                      tGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACgg                                                 <  1:609784/150‑1 (MQ=255)

                                                                                                         tCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCtt                                                                 >  1:4348705/1‑131 (MQ=255)

                                                                                                          cAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAAt                                             >  2:5007619/1‑150 (MQ=255)

                                                                                                            aGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGc                                           >  2:4374888/1‑150 (MQ=255)

                                                                                                             ggACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCt                                          <  1:1681471/150‑1 (MQ=255)

                                                                                                             ggACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCt                                          >  1:4408316/1‑150 (MQ=255)

                                                                                                             ggACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCt                                          >  2:3755216/1‑150 (MQ=255)

                                                                                                              gACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTg                                         >  2:2101269/1‑150 (MQ=255)

                                                                                                              gACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTg                                         >  2:2427908/1‑150 (MQ=255)

                                                                                                                 aaGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggcc                                      >  2:450447/1‑150 (MQ=255)

                                                                                                                 aaGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggcc                                      >  1:3184410/1‑150 (MQ=255)

                                                                                                                 aaGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggcc                                      >  2:757165/1‑150 (MQ=255)

                                                                                                                  aGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccg                                     <  2:1426622/150‑1 (MQ=255)

                                                                                                                   gACCGCTGTGCTGCACAGCCCGCAGGTCGTGCTGCACGGCGACAGCCCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCGGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccgg                                    >  1:2496892/1‑150 (MQ=255)

                                                                                                                   gACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccgg                                    >  2:3968466/1‑150 (MQ=255)

                                                                                                                   gACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccgg                                    >  2:920327/1‑150 (MQ=255)

                                                                                                                    aCCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccggc                                   <  1:3122980/150‑1 (MQ=255)

                                                                                                                    aCCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTAGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCAGAAACGCTGGCTAGGCCCGGTCGCCACGGCAATGCTggccggc                                   >  2:2318339/1‑150 (MQ=255)

                                                                                                                     ccGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTggccggcc                                  <  2:5279809/150‑1 (MQ=255)

                                                                                                                      cGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCAGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCAGGCCGGCCa                                 >  1:1844498/1‑150 (MQ=255)

                                                                                                                         tgtgCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCatat                              >  1:1569393/1‑150 (MQ=255)

                                                                                                                         tgtgCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCatat                              <  1:5418453/150‑1 (MQ=255)

                                                                                                                           tgctgcACAGCCCGCAGGTCGTGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTTCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATcc                            >  1:3032407/1‑150 (MQ=255)

                                                                                                                               gcACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCt                        <  1:1542672/150‑1 (MQ=255)

                                                                                                                               gcACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCt                        <  1:1542791/150‑1 (MQ=255)

                                                                                                                                  caGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTgcg                     >  2:3523818/1‑150 (MQ=255)

                                                                                                                                   aGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTgcgc                    >  1:1135834/1‑150 (MQ=255)

                                                                                                                                       cGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctc                >  2:4771161/1‑150 (MQ=255)

                                                                                                                                       aacAGGTCGAGCAGTACGGCGCTGTCGCCGTGCTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctc                >  1:1467068/3‑150 (MQ=38)

                                                                                                                                        gCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctcg               >  1:5020169/1‑150 (MQ=255)

                                                                                                                                        gCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctcg               >  1:5148835/1‑150 (MQ=255)

                                                                                                                                          aGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctcgct             <  2:5148835/150‑1 (MQ=255)

                                                                                                                                          aGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGcgctcgct             <  1:4297067/150‑1 (MQ=255)

                                                                                                                                           ggTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTg            <  1:2831514/150‑1 (MQ=255)

                                                                                                                                           ggTCGAGCAGTACGGCGCTGTCGCCGTGCAGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTg            >  1:229334/1‑150 (MQ=38)

                                                                                                                                            gTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCT‑GACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGa          >  2:1614453/1‑150 (MQ=255)

                                                                                                                                            gTCGAGCAGTACGGCGCTGTCGCCGTGCTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTgg           >  1:2467762/1‑150 (MQ=39)

                                                                                                                                            gTCGAGCAGTACGGCGCTGTCGCCGTGCAGCTCGACCTGCGGGCTGTGCAGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTgg           >  1:4153817/1‑150 (MQ=25)

                                                                                                                                               gAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGcc                                                             >  1:388226/1‑97 (MQ=255)

                                                                                                                                               gAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGcc                                                             <  2:388226/97‑1 (MQ=255)

                                                                                                                                               gAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGACt        >  1:2089688/1‑150 (MQ=255)

                                                                                                                                               gAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCACGCTGGACt        >  1:2153122/1‑150 (MQ=255)

                                                                                                                                                 gCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGACTCg      >  2:3379932/1‑150 (MQ=255)

                                                                                                                                                  cTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGACTCGc     <  2:2545592/150‑1 (MQ=255)

                                                                                                                                                   tGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTTGc                                                              >  2:1088034/1‑92 (MQ=255)

                                                                                                                                                   tGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTTGc                                                              <  1:1088034/92‑1 (MQ=255)

TCATCCTGCCCGATGAAGAGGGCCTGATGCTCGCCGCCCATTGGCAGGGCGAGGTCGACTGCCGCCTGCGCGCCCCCGCCGGTAGCCTGGCGCAACTGGCCTTGGTCAAGGACAAGACCGCTGTGCTGCACAGCCCGCAGGTCGAGCTGCACGGCGACAGCGCCGTACTGCTCGACCTGTTCGGTGTGCTGCAGGACCTGGAGCTGGACTGGGAGCACGAGCTGCAACGCTGGCTTGGCCCGGTCGCCACGGCAATGCTGGCCGGCCATATCCGCCTGCGCGCTCGCTGGACTCGCCAA  >  NC_002947/5714052‑5714350

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: `‑`