BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,955,992	T→G	5.1%	K657T (AAG→ACG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,955,992	0	T	G	5.1%	1312.5 / 39.2	383	K657T (AAG→ACG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand): ref base T (193/170); new base G (12/8); total (205/178)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 390 total aligned reads displayed.

TCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCGGG  >  NC_002947/5955843‑5956141

 cTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCtt                                                                                                                                                      >  1:8331080/1‑150 (MQ=255)

 cTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCtt                                                                                                                                                      <  1:1448638/150‑1 (MQ=255)

   ggTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGt                                                                                                                                                    <  1:5675905/150‑1 (MQ=255)

    gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTc                                                                                                                                                   <  2:6734617/150‑1 (MQ=255)

     ttGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCg                                                                                                                                                  >  2:491832/1‑150 (MQ=255)

     ttGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCg                                                                                                                                                  >  2:2177657/1‑150 (MQ=255)

       gCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                                <  1:2109164/150‑1 (MQ=255)

        ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                               <  1:3205872/150‑1 (MQ=255)

        ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                               >  2:5500465/1‑150 (MQ=255)

         cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCGTGTCGAGGt                                                                                                                                              >  2:2742356/1‑150 (MQ=255)

          ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgacc                                                                                                                                             <  1:523260/150‑4 (MQ=255)

             cGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcg                                                                                                                                          <  1:7489417/150‑1 (MQ=255)

              ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCTAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                         <  1:1006762/150‑1 (MQ=255)

              ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                         <  2:818422/150‑1 (MQ=255)

              ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                         <  2:7470795/150‑1 (MQ=255)

              ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                         >  1:1522090/1‑150 (MQ=255)

              ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                         >  2:253379/1‑150 (MQ=255)

                  tgctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcca                                                                                                                                     <  2:8275924/150‑1 (MQ=255)

                     tgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatcc                                                                                                                                  >  2:7771195/1‑150 (MQ=255)

                     tgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatcc                                                                                                                                  >  2:8488765/1‑150 (MQ=255)

                     tgcGCACGGAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCGTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatcc                                                                                                                                  <  2:6062350/150‑1 (MQ=255)

                          aCGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCg                                                                                                                             >  2:3239256/1‑150 (MQ=255)

                              aGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgcgc                                                                                                                         <  2:6577388/150‑1 (MQ=255)

                               gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTTCGCTCCATCCAGTCGGCGCt                                                                                                                        >  2:4497818/1‑150 (MQ=255)

                               gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                                                        >  2:2332844/1‑150 (MQ=255)

                                cgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTg                                                                                                                       <  1:1610541/150‑1 (MQ=255)

                                  cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTggn                                                                                                                     <  1:1203431/150‑2 (MQ=255)

                                  cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGa                                                                                                                     <  1:2937765/150‑1 (MQ=255)

                                   ccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                                                    <  1:5878774/150‑1 (MQ=255)

                                   ccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                                                    <  2:93370/150‑1 (MQ=255)

                                    cGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACa                                                                                                                   <  2:5428900/150‑1 (MQ=255)

                                     gTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAg                                                                                                                  >  2:5264662/1‑150 (MQ=255)

                                      tCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAgg                                                                                                                 <  2:3503035/150‑1 (MQ=255)

                                         ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAg                                                                                                              <  1:2468977/150‑1 (MQ=255)

                                                  gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCTCCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGctct                                                                                                     >  2:1935478/1‑150 (MQ=255)

                                                   tCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCtt                                                                                                    >  2:2582135/1‑150 (MQ=255)

                                                   tCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCtt                                                                                                    >  2:2050043/1‑150 (MQ=255)

                                                   tCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCtt                                                                                                    >  2:6603550/1‑150 (MQ=255)

                                                      aGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcgc                                                                                                 <  1:537455/150‑1 (MQ=255)

                                                         cGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcgccgc                                                                                              <  1:5022625/150‑1 (MQ=255)

                                                             ggTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCAtt                                                                                          <  2:7373654/150‑1 (MQ=255)

                                                               tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAAGAGTTTCTTGCCTTCCACCTGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGATCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                                                        <  2:1526093/150‑1 (MQ=255)

                                                               tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAggcct                                                                                                                                            >  1:6881844/1‑95 (MQ=255)

                                                               tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                                                        <  1:4778850/150‑1 (MQ=255)

                                                                 ttGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAgg                                                                                      <  1:449252/150‑1 (MQ=255)

                                                                  tGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGa                                                                                     <  2:1039605/150‑1 (MQ=255)

                                                                   gTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGaa                                                                                    >  1:8578601/1‑150 (MQ=255)

                                                                    tCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAg                                                                                   >  2:4682678/1‑150 (MQ=255)

                                                                      ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTa                                                                                 >  1:6621800/1‑150 (MQ=255)

                                                                      ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTa                                                                                 <  1:6316719/150‑1 (MQ=255)

                                                                            gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAggcgct                                                                           >  1:3989520/1‑145 (MQ=255)

                                                                            gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAacac                                                                           <  1:643910/150‑1 (MQ=255)

                                                                            gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCCCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAacac                                                                           >  2:3921059/1‑150 (MQ=255)

                                                                                ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCcgc                                                                       >  2:5614727/1‑150 (MQ=255)

                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGTTCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCgt                                                                     >  2:7153636/1‑150 (MQ=255)

                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCgt                                                                     <  2:5983702/150‑1 (MQ=255)

                                                                                          ccTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                                                             <  1:5587554/150‑1 (MQ=255)

                                                                                           cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGACCACCCGCGTGTGCTCGAg                                                            >  1:1707120/1‑150 (MQ=255)

                                                                                           cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAg                                                            >  1:1900873/1‑150 (MQ=255)

                                                                                            ttcttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                           >  2:4532410/1‑150 (MQ=255)

                                                                                            ttcttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                           >  1:5186133/1‑150 (MQ=255)

                                                                                            ttcttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                           >  2:1241348/1‑150 (MQ=255)

                                                                                            ttcttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCAGTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                           >  2:6798955/1‑150 (MQ=255)

                                                                                               ttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGACGTAGAACACCCGCGTGTGCTCGAGGAAg                                                        >  2:7889633/1‑150 (MQ=255)

                                                                                                tcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGc                                                       >  1:3795738/1‑150 (MQ=255)

                                                                                                   cccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAagcagc                                                    >  2:7617744/1‑150 (MQ=255)

                                                                                                   cccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGc                                                    >  2:8541641/1‑150 (MQ=255)

                                                                                                      gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCaa                                                 >  1:2564241/1‑150 (MQ=255)

                                                                                                           aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc                                            >  2:1787298/1‑150 (MQ=255)

                                                                                                           aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc                                            <  1:1274357/150‑1 (MQ=255)

                                                                                                              aGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAg                                         <  2:6174444/150‑1 (MQ=255)

                                                                                                              aGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAg                                         <  1:1094923/150‑1 (MQ=255)

                                                                                                               gTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAgc                                        <  2:6935706/150‑1 (MQ=255)

                                                                                                                 ttCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAgcgc                                      <  1:4429035/150‑1 (MQ=255)

                                                                                                                 ttCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAgcgc                                      <  1:3944362/150‑1 (MQ=255)

                                                                                                                  tCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCa                                     >  1:6584193/1‑150 (MQ=255)

                                                                                                                    ttGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCAcc                                   >  2:4160803/1‑150 (MQ=255)

                                                                                                                      gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTg                                 <  1:6025832/150‑1 (MQ=255)

                                                                                                                       ccTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTgg                                >  1:4240547/1‑150 (MQ=255)

                                                                                                                        cTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGa                               <  2:7537264/150‑1 (MQ=255)

                                                                                                                         ttCCACCGGGAAGCAGGTCTCGGCACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAt                              >  1:2807772/1‑150 (MQ=255)

                                                                                                                         ttCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAt                              >  1:4555599/1‑150 (MQ=255)

                                                                                                                         ttCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAt                              >  2:140323/1‑150 (MQ=255)

                                                                                                                            cACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgtt                           <  2:433419/150‑1 (MQ=255)

                                                                                                                             aCCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttg                          <  1:3453400/150‑1 (MQ=255)

                                                                                                                             aCCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttg                          >  1:6237851/1‑150 (MQ=255)

                                                                                                                               cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTtgtg                        >  1:6525046/1‑150 (MQ=255)

                                                                                                                                gggAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc                       <  2:225495/150‑1 (MQ=255)

                                                                                                                                gggAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc                       <  1:2285483/150‑1 (MQ=255)

                                                                                                                                gggAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc                       >  1:7914528/1‑150 (MQ=255)

                                                                                                                                   aaGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGaa                    <  1:4652011/150‑1 (MQ=255)

                                                                                                                                      gcaacctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTAGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCGGGATGTTGCGCGAAAcc                 >  1:238619/6‑150 (MQ=255)

                                                                                                                                          tctcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCg             >  1:7646574/1‑150 (MQ=255)

                                                                                                                                             cGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGaa          <  2:6311906/150‑1 (MQ=255)

                                                                                                                                              gACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAAt         <  2:33442/150‑1 (MQ=255)

                                                                                                                                                 acGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCa                                                                                            <  1:3121082/64‑1 (MQ=255)

                                                                                                                                                 acGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCa                                                                                            >  2:3121082/1‑64 (MQ=255)

                                                                                                                                                 acGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGcc      <  1:7060066/150‑1 (MQ=255)

                                                                                                                                                  cGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGccc     <  2:2885920/150‑1 (MQ=255)

                                                                                                                                                   gCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCg    >  1:893320/1‑150 (MQ=255)

TCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCGGG  >  NC_002947/5955843‑5956141

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