BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,287	C→A	5.3%	S865S (TCG→TCT)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,287	0	C	A	5.3%	1333.8 / 43.1	393	S865S (TCG→TCT)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base C (175/197); new base A (11/10); total (186/207)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 416 total aligned reads displayed.

CTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGC  >  NC_002947/3803138‑3803434

 tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                                    >  2:3389462/1‑150 (MQ=255)

 tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                                    >  2:4630371/1‑150 (MQ=255)

 tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                                    >  1:1445449/1‑150 (MQ=255)

 tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                                    >  2:64112/1‑150 (MQ=255)

 tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                                    >  2:7695041/1‑150 (MQ=255)

   aaGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGAc                                                                                                                                                  >  2:8725330/1‑150 (MQ=255)

    aGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACa                                                                                                                                                 >  1:2279914/1‑150 (MQ=255)

     ggCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACAt                                                                                                                                                >  1:6284047/1‑150 (MQ=255)

        gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                                                                             <  1:678251/150‑1 (MQ=255)

        gTTATCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACGACCGCCGACATCgg                                                                                                                                             <  1:5399966/150‑1 (MQ=255)

         ttttCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCggg                                                                                                                                            >  1:5865009/1‑150 (MQ=255)

              tttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGcc                                                                                                                                       >  2:2681284/1‑150 (MQ=255)

                tcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGt                                                                                                                                     <  2:3730009/150‑1 (MQ=255)

                   aCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAc                                                                                                                                  >  1:4414077/1‑150 (MQ=255)

                   aCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAc                                                                                                                                  <  2:4832818/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                    gCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATg                                                                                                                 <  2:2926901/150‑1 (MQ=255)

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

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

                                        tgAGCACTGCGAGGCCTGCGGCACAGGAAACCGGATTGCCACCGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCAtagtga                                                                                                             >  1:8572556/2‑145 (MQ=255)

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

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

                                              ccGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCAt                                                                                                       <  2:6724210/150‑1 (MQ=255)

                                                 ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                                    >  1:5144622/1‑150 (MQ=255)

                                                 ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                                    >  2:8381523/1‑150 (MQ=255)

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

                                                    ggCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGa                                                                                                 <  1:5338418/150‑1 (MQ=255)

                                                    ggCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGa                                                                                                 >  1:3220810/1‑150 (MQ=255)

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

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

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

                                                           cgcgCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGcc                                                                                          >  1:4096839/1‑150 (MQ=255)

                                                            gcgcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGccc                                                                                         >  1:2904880/1‑150 (MQ=255)

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

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

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

                                                                      ccGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGa                                                                               <  1:1172893/150‑1 (MQ=255)

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

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

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

                                                                         ggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATg                                                                            <  2:704610/150‑1 (MQ=255)

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

                                                                              ccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccac                                                                       >  1:3524834/5‑150 (MQ=255)

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

                                                                                   gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                                  <  2:45760/150‑1 (MQ=255)

                                                                                   gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                                  <  2:2231026/150‑1 (MQ=255)

                                                                                      ggTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCACCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                                               <  2:2020185/150‑1 (MQ=255)

                                                                                        tATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGc                                                             <  2:8549005/150‑1 (MQ=255)

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

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

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

                                                                                             aaaTACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAcc                                                        >  1:5065788/1‑150 (MQ=255)

                                                                                              aaTATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTAcgc                                                       <  1:3967686/150‑3 (MQ=255)

                                                                                               aTACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCg                                                      >  2:505673/1‑150 (MQ=255)

                                                                                                tACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGGATGCCCACCAGTGGCTGCCGACCCgg                                                     >  2:1531051/1‑150 (MQ=255)

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

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

                                                                                                     cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGTGCGTGGTCACGACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                <  2:4021833/150‑1 (MQ=255)

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

                                                                                                     cATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                >  1:6972275/1‑150 (MQ=255)

                                                                                                        gCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAccc                                             >  1:4267922/1‑150 (MQ=255)

                                                                                                          cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                           <  2:8258583/150‑1 (MQ=255)

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

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

                                                                                                            ttgttgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGAcc                                         <  2:8704878/150‑1 (MQ=255)

                                                                                                                tgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGaa                                     >  2:3847569/1‑150 (MQ=255)

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

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

                                                                                                                 gAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGG‑CGTTAACAATGGGTTTGCCCATGCTTACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTTCCGACCCGGCCGAACCCGACCTGAAcc                                   >  1:836240/1‑150 (MQ=255)

                                                                                                                   aaGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCt                                  >  1:6583771/1‑150 (MQ=255)

                                                                                                                    aGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTc                                 >  2:8329027/1‑150 (MQ=255)

                                                                                                                    aGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTc                                 <  1:4195060/150‑1 (MQ=255)

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

                                                                                                                        aTCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATc                              >  1:942357/1‑150 (MQ=255)

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

                                                                                                                            gCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGcc                          <  1:2725412/150‑1 (MQ=255)

                                                                                                                            gCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGcc                          >  2:3925691/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                          ggTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCAcc            <  1:8716696/150‑1 (MQ=255)

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

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

                                                                                                                                                tACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc      >  2:4363911/1‑150 (MQ=255)

                                                                                                                                                tACCGCCGACATCGGGTGGCCGTTACAAATGGGTTTGCCCATGCAGACGGAATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc      >  1:1919567/1‑150 (MQ=255)

                                                                                                                                                  ccgccgACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTc    <  1:2813634/150‑1 (MQ=255)

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

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

                                                                                                                                                   cgccgACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCg   >  2:4784470/1‑150 (MQ=255)

CTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGC  >  NC_002947/3803138‑3803434

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