BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,244,886	T→C	5.8%	P142P (CCA→CCG)	ldh ←	leucine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,244,886	0	T	C	5.8%	1124.2 / 35.2	348	P142P (CCA→CCG)	ldh	leucine dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (131/195); new base C (8/12); total (139/209)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 352 total aligned reads displayed.

CGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTCGAACAG  >  NC_002947/5244739‑5245035

cgcgCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTgg                                                                                                                                                     <  1:477064/150‑1 (MQ=255)

 gcgcCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGc                                                                                                                                                    >  1:2924457/1‑150 (MQ=255)

      cACCGCATGCAATTGCTCCGCCAGGGCATAACCAAGATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGATgg                                                                                                                                               <  2:1205744/150‑1 (MQ=255)

        ccGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACCCCCCCATCGCCACATGCGGCGATgggg                                                                                                                                             <  2:7154978/150‑2 (MQ=255)

        ccGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgggg                                                                                                                                             <  2:6216689/150‑2 (MQ=255)

          gCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGATgggggt                                                                                                                                           <  1:6633943/150‑4 (MQ=255)

          gCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGATgggggt                                                                                                                                           <  2:1010942/150‑4 (MQ=255)

                aaTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGAt                                                                                                                                     >  2:450339/1‑150 (MQ=255)

                aaTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGAt                                                                                                                                     <  2:4713559/150‑1 (MQ=255)

                  ttGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCTGGTCATCGCTGCCCAGCCGTAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGTTGGGTCGCCCGATGc                                                                                                                                   <  1:6491240/150‑1 (MQ=255)

                   tGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCa                                                                                                                                  <  2:2115493/150‑1 (MQ=255)

                   tGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCa                                                                                                                                  <  2:1587987/150‑1 (MQ=255)

                      tCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgggg                                                                                                                                             <  2:324520/136‑2 (MQ=255)

                      tCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgggg                                                                                                                                             >  1:324520/1‑135 (MQ=255)

                         gCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCCGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtggtg                                                                                                                            <  1:5911693/150‑1 (MQ=255)

                          ccAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGc                                                                                                                           >  2:6420337/1‑150 (MQ=255)

                           cAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgg                                                                                                                                               <  1:3068059/129‑1 (MQ=255)

                            aGGGCATACCACACATGGCCCAACCCGTGGACCGGTCCCGGCTAGCCTCTCAGGTAATCGGTGCCCCGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATTTGGCTATGGGTCGCCCGATGCAGTGGTGCTg                                                                                                                         <  1:2012479/150‑1 (MQ=255)

                            aGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTg                                                                                                                         <  2:6464015/150‑1 (MQ=255)

                            aGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCGGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTg                                                                                                                         <  2:5989368/150‑1 (MQ=255)

                             gggCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTgg                                                                                                                        <  1:3688003/150‑1 (MQ=255)

                             gggCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTgg                                                                                                                        <  1:4188415/150‑1 (MQ=255)

                                    aCCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGCCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTa                                                                                                                 >  1:5881717/1‑150 (MQ=255)

                                      cAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAgg                                                                                                               >  1:2338809/1‑150 (MQ=255)

                                      cAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgggg                                                                                                                                             <  1:1697016/120‑2 (MQ=255)

                                         cATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCg                                                                                                            >  2:239795/1‑150 (MQ=255)

                                               ccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                      <  2:1750666/150‑1 (MQ=255)

                                                 aaCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTCTGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                                    <  1:3224031/150‑1 (MQ=255)

                                                 aaCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                                    <  2:6757273/150‑1 (MQ=255)

                                                 aaCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                                    <  2:4119154/150‑1 (MQ=255)

                                                 aaCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                                    >  1:3364569/1‑150 (MQ=255)

                                                  aCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTggg                                                                                                   <  1:5759036/150‑1 (MQ=255)

                                                     ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                                <  1:5492613/150‑1 (MQ=255)

                                                     ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                                <  1:4635164/150‑1 (MQ=255)

                                                     ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTATGGGCaa                                                                                                <  2:2770485/150‑1 (MQ=255)

                                                      cTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCCAt                                                                                               >  2:2925088/1‑150 (MQ=255)

                                                         gACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                            <  2:254424/150‑1 (MQ=255)

                                                         gACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                            <  1:6081586/150‑1 (MQ=255)

                                                              cTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCa                                                                                       >  2:7050677/1‑150 (MQ=255)

                                                              cTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCa                                                                                       >  1:6199231/1‑150 (MQ=255)

                                                               tACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                                      >  2:6121929/1‑150 (MQ=255)

                                                               tACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                                      <  1:2755053/150‑1 (MQ=255)

                                                               tACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                                      <  2:3606236/150‑1 (MQ=255)

                                                                    gCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCa                                                                                 >  1:1014204/1‑150 (MQ=255)

                                                                    gCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCa                                                                                 >  2:4532462/1‑150 (MQ=255)

                                                                     cAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCaa                                                                                <  2:3554455/150‑1 (MQ=255)

                                                                     cAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCaa                                                                                >  2:7103674/1‑150 (MQ=255)

                                                                     cAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCaa                                                                                <  2:3809564/150‑1 (MQ=255)

                                                                        gCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGt                                                                             >  2:6269809/1‑150 (MQ=255)

                                                                        gCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGt                                                                             >  1:4429362/1‑150 (MQ=255)

                                                                          ctctCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCg                                                                           >  1:4525728/1‑150 (MQ=255)

                                                                           tctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGa                                                                          <  1:788295/150‑1 (MQ=255)

                                                                            ctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                         <  2:3364569/150‑1 (MQ=255)

                                                                            ctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                         <  1:133378/150‑1 (MQ=255)

                                                                            ctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                         <  2:4587453/150‑1 (MQ=255)

                                                                                gttCATCGCTGCCCAGCCGGAACGAGGTTGTGGAACGTATGCCCGCGAACACCCCAATCGCCGCATGTTGCGTTGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGc                                                                     <  1:7221833/148‑1 (MQ=255)

                                                                                ggTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                            >  2:1381186/1‑127 (MQ=255)

                                                                                 gTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGcc                                                                    <  1:2194975/150‑1 (MQ=255)

                                                                                  tCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCa                                                                   <  2:3722839/150‑1 (MQ=255)

                                                                                   cATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCAc                                                                  <  2:5476156/150‑1 (MQ=255)

                                                                                    aTCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACt                                                                 >  1:5575903/1‑150 (MQ=255)

                                                                                      cGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGt                                                               <  2:3890775/150‑1 (MQ=255)

                                                                                      cGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGt                                                               <  1:6280265/150‑1 (MQ=255)

                                                                                          gCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCAc                                                           <  1:5159684/150‑1 (MQ=255)

                                                                                           cccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACg                                                          <  2:1712233/150‑1 (MQ=255)

                                                                                            ccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACgg                                                         <  2:2521155/150‑1 (MQ=255)

                                                                                               gCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCga                                                      >  2:786180/1‑150 (MQ=255)

                                                                                                ccGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCTTGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgat                                                     <  1:2834317/150‑1 (MQ=255)

                                                                                                  ggAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatga                                                   <  1:490759/150‑1 (MQ=255)

                                                                                                   gAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatgat                                                  <  2:308818/150‑1 (MQ=255)

                                                                                                    aaCGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGatgatg                                                 <  2:6994165/150‑1 (MQ=255)

                                                                                                      cgcggTTGGGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                               <  1:7324873/147‑1 (MQ=255)

                                                                                                       gAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAg                                              >  2:1078860/1‑150 (MQ=255)

                                                                                                          gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                           <  2:1910369/150‑1 (MQ=255)

                                                                                                          gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                           <  2:4256848/150‑1 (MQ=255)

                                                                                                          gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                           >  2:6154983/1‑150 (MQ=255)

                                                                                                          gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                           <  2:539223/150‑1 (MQ=255)

                                                                                                          gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCACGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                           >  1:5543347/1‑150 (MQ=255)

                                                                                                           ttGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGcc                                          <  1:5394933/150‑1 (MQ=255)

                                                                                                              tgGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCAt                                       >  1:380231/1‑150 (MQ=255)

                                                                                                               ggCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATg                                      <  2:1935961/150‑1 (MQ=255)

                                                                                                                gCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGCTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGc                                     <  1:352772/150‑1 (MQ=255)

                                                                                                                  aCGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCaa                                   <  1:6795743/150‑1 (MQ=255)

                                                                                                                  aCGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCaa                                   <  2:4404995/150‑1 (MQ=255)

                                                                                                                        gCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATc                             >  2:5040526/1‑150 (MQ=255)

                                                                                                                         cccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCg                            >  1:1764014/1‑150 (MQ=255)

                                                                                                                         cccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCg                            <  1:2044375/150‑1 (MQ=255)

                                                                                                                          ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACt                                                                 <  2:6653059/112‑1 (MQ=255)

                                                                                                                          ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGa                           <  1:4102082/150‑1 (MQ=255)

                                                                                                                           cgcgAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGAt                          <  1:5923181/150‑1 (MQ=255)

                                                                                                                              gAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGaa                       <  1:925305/150‑1 (MQ=255)

                                                                                                                                 cacCCCCATCGCCGCATGTGGCGTTGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCTCGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCg                    <  2:5145249/150‑1 (MQ=255)

                                                                                                                                 cacCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCg                    <  2:1679412/150‑1 (MQ=255)

                                                                                                                                   cccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc                  <  1:5260536/150‑1 (MQ=255)

                                                                                                                                    ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAgg                                 <  1:1450912/134‑1 (MQ=255)

                                                                                                                                    ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGcc                 >  2:5642139/1‑150 (MQ=255)

                                                                                                                                       cATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCaaa              <  1:1224781/150‑1 (MQ=255)

                                                                                                                                         tCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAgg            >  1:637651/1‑150 (MQ=255)

                                                                                                                                            ccgcTTGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCt         <  2:5247736/150‑1 (MQ=255)

                                                                                                                                            ccgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCAc                                                           >  2:4926765/1‑100 (MQ=255)

CGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTCGAACAG  >  NC_002947/5244739‑5245035

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