BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,293	C→T	5.2%	P863P (CCG→CCA)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,293	0	C	T	5.2%	1287.6 / 35.5	402	P863P (CCG→CCA)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base C (177/203); new base T (11/10); total (189/213)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 409 total aligned reads displayed.

GCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGC  >  NC_002947/3803144‑3803442

gCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATc                                                                                                                                                       <  2:3751324/150‑1 (MQ=255)

gCGTTTTCTTCCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATc                                                                                                                                                       <  2:5040418/150‑1 (MQ=255)

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

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

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

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

     ttctttcAACTGGTCGCGCTCGATTGCGTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTg                                                                                                                                                  <  2:2595867/150‑1 (MQ=255)

      tctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTgg                                                                                                                                                 <  1:4986829/150‑1 (MQ=255)

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

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

            aaCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATTCTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTa                                                                                                                                           >  2:4579299/1‑150 (MQ=255)

            aaCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTa                                                                                                                                           >  2:7349851/1‑150 (MQ=255)

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

                ggTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCaa                                                                                                                                       >  2:32511/1‑150 (MQ=255)

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

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

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

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

                          gATTGCATCGAGTACCGCCAGGCCTACCGCGCGTGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                                             >  1:6441347/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                aTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCt                                                                                                                       >  2:2453651/1‑150 (MQ=255)

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

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

                                  cGAGTACCGCCAGGCCTACCACGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGa                                                                                                                     <  2:1785023/150‑1 (MQ=255)

                                   gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                                                    >  2:5183801/1‑150 (MQ=255)

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

                                    aGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACg                                                                                                                   <  1:1498238/150‑1 (MQ=255)

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

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

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

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

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

                                                  tACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCAc                                                                                                     >  1:1644290/1‑150 (MQ=255)

                                                  tACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCAc                                                                                                     >  1:263249/1‑150 (MQ=255)

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

                                                        gcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTg                                                                                               <  2:3914303/150‑1 (MQ=255)

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

                                                            gACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGtc                                                                                           >  2:8703506/1‑150 (MQ=255)

                                                               acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                        >  2:5622266/1‑150 (MQ=255)

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

                                                                  tggTTACCGGCGAATTTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAAt                                                                                     <  2:1617028/149‑1 (MQ=255)

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

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

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

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

                                                                   ggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATg                                                                                    >  2:6679506/1‑150 (MQ=255)

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

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

                                                                        ccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccac                                                                               >  2:5872473/5‑150 (MQ=255)

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

                                                                            cGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGt                                                                           <  2:481781/150‑1 (MQ=255)

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

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

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

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

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

                                                                                        aaTACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAccc                                                               <  2:6062576/150‑1 (MQ=255)

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

                                                                                          tACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCAgg                                                             <  1:1598880/150‑1 (MQ=255)

                                                                                            cTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGTCGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGcc                                                           >  2:1508085/1‑150 (MQ=255)

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

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

                                                                                                aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                       <  2:5241399/150‑1 (MQ=255)

                                                                                                aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                       <  1:2150716/150‑1 (MQ=255)

                                                                                                aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                       <  2:7599934/150‑1 (MQ=255)

                                                                                                aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                       <  1:8780707/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                               gC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCa                                        >  2:1023131/1‑150 (MQ=255)

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

                                                                                                                    cGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCgg                                    >  1:4576491/1‑150 (MQ=255)

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

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

                                                                                                                              gcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAg                          <  2:6915560/150‑1 (MQ=255)

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

                                                                                                                                      tCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCtt                  <  2:3008727/150‑1 (MQ=255)

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

                                                                                                                                        aCTACCGCCCACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCCGCCGAACCCGACCAGAACCTCAACGGCCAGGCACAGCCCCCCATCg                >  1:2446534/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

GCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGC  >  NC_002947/3803144‑3803442

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