BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,241,854	G→A	7.7%	A144V (GCG→GTG)	ldh ←	leucine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,241,854	0	G	A	7.7%	579.6 / 26.0	194	A144V (GCG→GTG)	ldh	leucine dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (95/84); new base A (5/10); total (100/94)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 198 total aligned reads displayed.

CAGCTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5241706‑5242002

cagcTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATcgccgc                                                                                                                                                     <  2:1905179/150‑1 (MQ=255)

     ccGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCAtgtg                                                                                                                                                >  2:4746077/1‑150 (MQ=255)

     ccGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCAtgtg                                                                                                                                                >  1:49138/1‑150 (MQ=255)

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

          cccACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGAt                                                                                                                                           <  1:5285595/150‑1 (MQ=255)

           ccACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATg                                                                                                                                          <  1:3165316/150‑1 (MQ=255)

            cACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgg                                                                                                                                         <  2:4195438/150‑1 (MQ=255)

            cACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgg                                                                                                                                         <  1:3905226/150‑1 (MQ=255)

             aCCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATggg                                                                                                                                        <  2:4064847/150‑1 (MQ=255)

             aCCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATggg                                                                                                                                        <  2:885392/150‑1 (MQ=255)

             aCCGCATGCAATTGCTCCGCCAGGACATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATggg                                                                                                                                        >  1:213330/1‑150 (MQ=255)

                  aTGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   <  2:2124524/150‑1 (MQ=255)

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

                           cTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgt                                                                                                                          >  1:2363810/1‑150 (MQ=255)

                            tCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                         >  1:1874834/1‑150 (MQ=255)

                            tCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGATgggggt                                                                                                                                     <  1:2728404/138‑4 (MQ=255)

                            tCCGCCAGGGCATAACCAACAGTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGACCGAGGTTGTGGCACGTATGCCCGCGAACACACCCATCGCCACATGCGGCGATgggggt                                                                                                                                     >  2:2728404/1‑135 (MQ=255)

                            cccGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGACCGTGGTTGTGGCACGTCTGCCCGCGAACACCCCCCTCGCCGCAGGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                         <  2:1949280/149‑1 (MQ=255)

                              cGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtggt                                                                                                                       >  1:1953997/1‑150 (MQ=255)

                               gCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtggtg                                                                                                                      >  1:2248384/1‑150 (MQ=255)

                                       aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCAc                                                                                                              >  1:3353705/1‑150 (MQ=255)

                                          aCCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTa                                                                                                           >  2:2329479/1‑150 (MQ=255)

                                          aCCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTa                                                                                                           >  2:2329644/1‑150 (MQ=255)

                                          aCCAACATTGCACAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTa                                                                                                           >  2:7794/1‑150 (MQ=255)

                                           ccAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAg                                                                                                          >  1:3250718/1‑150 (MQ=255)

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

                                            cAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCAGCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAgg                                                                                                         <  2:1161460/150‑1 (MQ=255)

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

                                                  tGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGc                                                                                                   >  1:2625832/1‑150 (MQ=255)

                                                   gCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCt                                                                                                  <  1:4509604/150‑1 (MQ=255)

                                                      cAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               <  2:2020317/150‑1 (MQ=255)

                                                      cAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               >  2:4890374/1‑150 (MQ=255)

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

                                                        aCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCGGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACCCCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTggg                                                                                             <  2:4454974/150‑1 (MQ=255)

                                                          cccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCa                                                                                           >  1:5245747/1‑150 (MQ=255)

                                                           tcagcgccGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                          <  1:1070953/144‑1 (MQ=255)

                                                            cTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAAt                                                                                         >  2:4238329/1‑150 (MQ=255)

                                                             tGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGcccgg                                                                                        <  2:398071/150‑5 (MQ=255)

                                                             tGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATg                                                                                        <  1:2080183/150‑1 (MQ=255)

                                                               gACCGCTACCCGCACGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGTTGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCa                                                                                      <  1:4678922/150‑1 (MQ=255)

                                                                aCCGCTACCCGCAAGCCTCTCAGGTCATCGCGGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAg                                                                                     <  1:808893/150‑1 (MQ=255)

                                                                 ccGCTACACGCAAGCCTCTCAGTTCATCGCTGCCCAGCCGGAAAGAGGTTGTGGCACGTATGCCCGCTAACACCCCCATCGCCGCATGTGGCTATGGGTCGCCCTATGCAGTGGTGCTTGTCACGTAGGGCGTGCTGTGTTCAATGCAGt                                                                                    <  2:4195252/150‑1 (MQ=255)

                                                                  cGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   >  2:3465100/1‑150 (MQ=255)

                                                                  cGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   >  2:4579042/1‑150 (MQ=255)

                                                                    cTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCa                                                                                 <  2:3250718/150‑1 (MQ=255)

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

                                                                     tACCCGCACGCCTCTCCGGTAATCGCTGCCCCGCCGGAAGGAGGTTGTGGCACTTATGCCCGCGAACACCCCCATCGCGGCTTGTGGCTATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                                <  1:955570/150‑1 (MQ=255)

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

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

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

                                                                        ccGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGTGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTc                                                                             >  2:1220027/1‑150 (MQ=255)

                                                                         cGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTcc                                                                            >  1:1718798/1‑150 (MQ=255)

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

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

                                                                                 tctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGa                                                                    >  1:4631862/1‑150 (MQ=255)

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

                                                                                  ctcAGGTCATCGCGGCCCAGCCGGAACGAGGTTGTGGCACGTAGGCCCGCGAACACCCCCATCGCCGCTTGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                   <  1:2963968/150‑1 (MQ=255)

                                                                                   tcATGTCATCGCTGCCCATCGGGATGGAGGTTGTGGCACGTATGCCCGCGAACACCCCCCTCGCCGCATTTGGCTATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAgg                                                                  <  1:4438256/150‑1 (MQ=255)

                                                                                    cAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGt                                                                 <  2:2466329/150‑1 (MQ=255)

                                                                                      ggTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGc                                                               <  1:4890423/150‑1 (MQ=255)

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

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

                                                                                            cGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGAGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGt                                                         >  2:1429863/1‑150 (MQ=255)

                                                                                             gCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgggg                                                                                                                                       >  2:1127936/1‑70 (MQ=255)

                                                                                               tGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCa                                                      >  2:3327107/1‑150 (MQ=255)

                                                                                                  ccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACgg                                                   >  1:1498227/1‑150 (MQ=255)

                                                                                                       cGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatg                                              >  1:4214644/1‑150 (MQ=255)

                                                                                                         gAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatgat                                            >  2:366892/1‑150 (MQ=255)

                                                                                                          aaCGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGatgatg                                           >  1:5383376/1‑150 (MQ=255)

                                                                                                            cGAGGTTGTGGCACGTATGCCCGCGAACACGCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  2:1040431/150‑1 (MQ=255)

                                                                                                            cGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  2:4580689/150‑1 (MQ=255)

                                                                                                            cGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  1:744169/150‑1 (MQ=255)

                                                                                                            cGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGGGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGAGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  1:2012005/150‑1 (MQ=255)

                                                                                                             gAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCTTGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCTGGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGTGGTGCCACTGTCCACGGCGTTGATGAAg                                        >  2:3165267/1‑150 (MQ=255)

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

                                                                                                                     ggCACGTTTGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATg                                <  1:102585/150‑1 (MQ=255)

                                                                                                                      gCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTc                                                                                                                >  2:2829196/1‑69 (MQ=255)

                                                                                                                      gCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTc                                                                                                                <  1:2829196/69‑1 (MQ=255)

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

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

                                                                                                                             tGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatga                                             >  2:552864/1‑129 (MQ=255)

                                                                                                                                ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGa                     >  1:1942406/1‑150 (MQ=255)

                                                                                                                                    gAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAg                            >  2:2208982/1‑139 (MQ=255)

                                                                                                                                    gAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAg                            <  1:2208982/139‑1 (MQ=255)

                                                                                                                                      acacCCCCATCGCCGCTTGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGc               <  1:932826/150‑1 (MQ=255)

                                                                                                                                        acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAg                            >  1:3683677/1‑135 (MQ=255)

                                                                                                                                        acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCgg             >  2:4181474/1‑150 (MQ=255)

                                                                                                                                        acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCgg             >  1:121039/1‑150 (MQ=255)

                                                                                                                                         gacccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc            >  1:4132430/3‑150 (MQ=255)

                                                                                                                                          ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAgg                                                                                                         >  2:2671984/1‑56 (MQ=255)

                                                                                                                                          ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTg                                                          <  1:2688457/103‑1 (MQ=255)

                                                                                                                                          ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGcc                                    >  1:3776252/1‑125 (MQ=255)

                                                                                                                                          ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGcc                                    <  2:3776252/125‑1 (MQ=255)

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

                                                                                                                                            ccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCATGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCaa         >  1:4194371/1‑150 (MQ=255)

                                                                                                                                            ccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCaa         >  1:315570/1‑150 (MQ=255)

                                                                                                                                            ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatga                                             <  1:2872902/114‑1 (MQ=255)

                                                                                                                                            ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatga                                             >  2:2872902/1‑114 (MQ=255)

                                                                                                                                              aTCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAg       <  2:5349721/150‑1 (MQ=255)

                                                                                                                                                   cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  <  1:602844/150‑1 (MQ=255)

CAGCTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5241706‑5242002

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