BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,299	G→A	5.3%	G861G (GGC→GGT)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,299	0	G	A	5.3%	1256.0 / 36.6	392	G861G (GGC→GGT)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base G (172/199); new base A (11/10); total (183/209)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 399 total aligned reads displayed.

TCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAC  >  NC_002947/3803150‑3803448

 ctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGc                                                                                                                                                      <  2:1780539/150‑1 (MQ=255)

 ctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGc                                                                                                                                                      <  2:8363376/150‑1 (MQ=255)

  tatcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGcc                                                                                                                                                     >  2:8269926/3‑150 (MQ=255)

         tGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCa                                                                                                                                              >  1:7660989/1‑150 (MQ=255)

          ggTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCaa                                                                                                                                             <  2:1013796/150‑1 (MQ=255)

           gTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAAt                                                                                                                                            >  2:7365952/1‑150 (MQ=255)

             cgcgCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATgg                                                                                                                                          <  2:1570151/150‑1 (MQ=255)

             cgcgCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATgg                                                                                                                                          <  1:3940269/150‑1 (MQ=255)

               cgcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGt                                                                                                                                        >  1:1236154/1‑150 (MQ=255)

                gcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGtt                                                                                                                                       >  2:1150504/1‑150 (MQ=255)

                gcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGtt                                                                                                                                       >  2:43634/1‑150 (MQ=255)

                gcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGtt                                                                                                                                       <  1:876566/150‑1 (MQ=255)

                gcTCGACTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGtt                                                                                                                                       <  1:6045462/150‑1 (MQ=255)

                  tCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTg                                                                                                                                     >  1:2982166/1‑150 (MQ=255)

                   cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                                                                                                    <  1:1580233/150‑1 (MQ=255)

                   cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATCCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                                                                                                    <  1:2864857/150‑1 (MQ=255)

                    gatgcaaTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                                                   >  2:2813634/7‑150 (MQ=255)

                    gATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCGCTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                                                   >  2:3160871/1‑150 (MQ=255)

                    gATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                                                   >  1:8071615/1‑150 (MQ=255)

                    gATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                                                   >  1:6872824/1‑150 (MQ=255)

                       tGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCAt                                                                                                                                <  1:3057219/150‑1 (MQ=255)

                       tGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCAt                                                                                                                                <  2:2697749/150‑1 (MQ=255)

                        acATTGAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCAta                                                                                                                               <  1:8247615/149‑2 (MQ=255)

                         cATTGAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCAtag                                                                                                                              <  1:1117086/150‑3 (MQ=255)

                         cATTGAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCAtag                                                                                                                              <  2:5681134/150‑3 (MQ=255)

                         cATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGc                                                                                                                              >  2:7066121/1‑150 (MQ=255)

                         cATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGc                                                                                                                              <  1:897669/150‑1 (MQ=255)

                           tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTg                                                                                                                            >  2:8476877/1‑150 (MQ=255)

                             gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                                                          <  2:5336902/150‑1 (MQ=255)

                             gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                                                          <  1:3774651/150‑1 (MQ=255)

                             gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                                                          <  2:6831933/150‑1 (MQ=255)

                              aGTACCGCCAGGCCTACCGCGCATGACACCGGGTGACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACg                                                                                                                         <  1:8693130/150‑1 (MQ=255)

                                   cGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATc                                                                                                                    >  1:535150/1‑150 (MQ=255)

                                   cGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATc                                                                                                                    >  2:8215116/1‑150 (MQ=255)

                                     ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTGTGCCCATGCTGACGGCATCCg                                                                                                                  >  2:460517/1‑150 (MQ=255)

                                     ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGCATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                                                  >  2:8655649/1‑150 (MQ=255)

                                      cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGt                                                                                                                                                   <  1:2737038/116‑1 (MQ=255)

                                      cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                 <  1:6230804/150‑1 (MQ=255)

                                      cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                 <  1:515149/150‑1 (MQ=255)

                                         gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                                              <  2:8747845/150‑1 (MQ=255)

                                         gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                                              <  1:8546073/150‑1 (MQ=255)

                                         gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                                              >  2:7808352/1‑150 (MQ=255)

                                                 cgcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCt                                                                                                      <  2:5528162/150‑1 (MQ=255)

                                                   cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                    <  1:4933545/150‑1 (MQ=255)

                                                   cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                    <  1:4589202/150‑1 (MQ=255)

                                                       acacCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGtct                                                                                                <  1:2561074/150‑1 (MQ=255)

                                                           cGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGaa                                                                                            <  2:3951690/150‑1 (MQ=255)

                                                            gggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAAt                                                                                           <  1:1021436/150‑1 (MQ=255)

                                                            gggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAAt                                                                                           <  2:6416876/150‑1 (MQ=255)

                                                              gTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGc                                                                                         <  1:4169939/150‑1 (MQ=255)

                                                                tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc                                                                                       <  1:4218005/150‑1 (MQ=255)

                                                                 gccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCcca                                                                                      <  1:7381731/145‑1 (MQ=255)

                                                                   caccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACCTTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacc                                                                                    >  2:4805122/4‑150 (MQ=255)

                                                                   caccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacc                                                                                    >  1:1401802/4‑150 (MQ=255)

                                                                    ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                   <  2:7968278/150‑1 (MQ=255)

                                                                          ggTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                                                             <  2:8134667/150‑1 (MQ=255)

                                                                           gTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTg                                                                            >  1:4319164/1‑150 (MQ=255)

                                                                              ttGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCg                                                                         >  1:5491139/1‑150 (MQ=255)

                                                                               tGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGa                                                                        <  1:2411287/150‑1 (MQ=255)

                                                                                    tACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCgg                                                                   >  2:4837408/1‑150 (MQ=255)

                                                                                     aCTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACGTCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGc                                                                  >  2:3248908/1‑150 (MQ=255)

                                                                                     aCTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGc                                                                  <  1:6558763/150‑1 (MQ=255)

                                                                                       tCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCg                                                                >  1:360787/1‑150 (MQ=255)

                                                                                       tCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCg                                                                >  2:897358/1‑150 (MQ=255)

                                                                                         cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                              >  2:4296849/1‑150 (MQ=255)

                                                                                         cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                              >  1:1638063/1‑150 (MQ=255)

                                                                                          ctGCCGTGGTTGAAAGC‑ATCGGCAACCTCGCGCGGGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                             <  1:2150662/149‑1 (MQ=255)

                                                                                               gttgttGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGAc                                                        <  2:5937433/150‑1 (MQ=255)

                                                                                                    tgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGaa                                                   <  2:387300/150‑1 (MQ=255)

                                                                                                     gAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAc                                                  <  1:4708761/150‑1 (MQ=255)

                                                                                                     gAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAc                                                  <  1:64112/150‑1 (MQ=255)

                                                                                                     gAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGG‑CGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAcc                                                 >  1:836137/1‑150 (MQ=255)

                                                                                                      aaaGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAcc                                                 >  1:5137489/1‑150 (MQ=255)

                                                                                                         gC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCa                                              <  2:4565078/150‑1 (MQ=255)

                                                                                                          c‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCAt                                             <  2:5195120/150‑1 (MQ=255)

                                                                                                             tCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCg                                           >  2:952161/1‑150 (MQ=255)

                                                                                                                 cAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCa                                       <  1:7041551/150‑1 (MQ=255)

                                                                                                                   aCCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAgg                                     <  1:7759698/150‑1 (MQ=255)

                                                                                                                     cTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGca                                   <  2:3640734/150‑1 (MQ=255)

                                                                                                                      tCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcac                                  >  2:1306465/1‑150 (MQ=255)

                                                                                                                       cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca                                 >  2:6248180/1‑150 (MQ=255)

                                                                                                                       cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca                                 <  2:655884/150‑1 (MQ=255)

                                                                                                                        gcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAg                                >  2:7324214/1‑150 (MQ=255)

                                                                                                                          gcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGcc                              >  1:1885499/1‑150 (MQ=255)

                                                                                                                                tCACGACCGCCGACATCGGGTGGCCGTTACCAAGGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCtt                        <  2:7196561/150‑1 (MQ=255)

                                                                                                                                 cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                       <  2:2279914/150‑1 (MQ=255)

                                                                                                                                  aCTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCg                      <  2:4096839/150‑1 (MQ=255)

                                                                                                                                  aCTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCg                      <  1:4079064/150‑1 (MQ=255)

                                                                                                                                   cTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCgg                     <  2:1466595/150‑1 (MQ=255)

                                                                                                                                   cTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCgg                     <  1:7186292/150‑1 (MQ=255)

                                                                                                                                      ccgccgACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTc                  >  2:2580877/1‑150 (MQ=255)

                                                                                                                                            aCATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCa            >  1:5785476/1‑150 (MQ=255)

                                                                                                                                             cATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCAt           <  1:3242728/150‑1 (MQ=255)

                                                                                                                                              aTCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGTCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATg          <  2:3353271/150‑1 (MQ=255)

                                                                                                                                              aTCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATg          <  2:3842796/150‑1 (MQ=255)

                                                                                                                                                cGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGc        <  1:446536/150‑1 (MQ=255)

                                                                                                                                                    tGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTac    >  2:511/1‑150 (MQ=255)

                                                                                                                                                      gCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTacac  <  2:7741892/150‑1 (MQ=255)

                                                                                                                                                      gCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTacac  <  2:2150025/150‑1 (MQ=255)

                                                                                                                                                      gCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTacac  <  2:3003824/150‑1 (MQ=255)

TCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAC  >  NC_002947/3803150‑3803448

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