BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,351	G→A	7.1%	T844I (ACC→ATC) ‡	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,351	0	G	A	7.1%	554.6 / 21.8	183	T844I (ACC→ATC) ‡	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base G (92/78); new base A (2/11); total (94/89)
Fisher's exact test for biased strand distribution p-value = 8.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 188 total aligned reads displayed.

ATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTGGGGATC  >  NC_002947/3803202‑3803500

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

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

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

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

           ttgccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                                            <  1:2170854/143‑1 (MQ=255)

           ttgccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                                            <  1:3021213/143‑1 (MQ=255)

           tagccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                                            <  1:3255659/143‑1 (MQ=255)

            tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc                                                                                                                                           >  1:1457912/1‑150 (MQ=255)

              ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                                         <  2:3898604/150‑1 (MQ=255)

              ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                                         <  1:2823616/150‑1 (MQ=255)

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

                ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                       >  1:2575615/1‑150 (MQ=255)

                accGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacca                                                                                                                                       <  1:2047916/148‑1 (MQ=255)

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

                   gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                                                                                                    >  2:4592070/1‑150 (MQ=255)

                      ggTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCt                                                                                                                                 <  2:1732259/150‑1 (MQ=255)

                         aTTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc                                                                                                                              >  1:717244/1‑150 (MQ=255)

                         aTTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc                                                                                                                              >  1:4455298/1‑150 (MQ=255)

                             aaaTATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTAcg                                                                                                                          <  2:61277/150‑2 (MQ=255)

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

                                tACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCgg                                                                                                                       <  2:4840898/150‑1 (MQ=255)

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

                                     cATGCCGTTGTTGAAATC‑ATCGGCAACCTCGCGGGTGTTCACTACCGCCTACATCTGGTGGGCTTTACCAATGTTTTTGCCCATGCTGTCGTCTTCCGGGTCCACTCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGat                                                                                                                  <  1:1601852/150‑2 (MQ=255)

                                     cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                                                                                  <  1:3477964/150‑1 (MQ=255)

                                     cATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                                  <  1:378038/150‑1 (MQ=255)

                                     cATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                                  <  1:1505374/150‑1 (MQ=255)

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

                                      aTGCCGTTGTTGAAAGC‑ATCGGCAACATCGCTCGTGGTCACTACCGCCGAAATCGGGTGGCCGTTACCAATGGGTTTGCCAATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGACGATGCCGACCCGGCCGAAc                                                                                                                 >  2:1863307/1‑150 (MQ=255)

                                        gCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAccc                                                                                                               <  1:1692125/150‑1 (MQ=255)

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

                                          cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                                             >  2:3417204/1‑150 (MQ=255)

                                            ttgttgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGAcc                                                                                                           >  2:1632493/1‑150 (MQ=255)

                                               ttgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGa                                                                                                        >  1:1618627/1‑150 (MQ=255)

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

                                                      c‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCAt                                                                                                 >  1:3809434/1‑150 (MQ=255)

                                                         tCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCg                                                                                               <  1:2902676/150‑1 (MQ=255)

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

                                                              aaCCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAg                                                                                          <  1:3926053/150‑1 (MQ=255)

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

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

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

                                                                       cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                                 <  1:1190677/150‑1 (MQ=255)

                                                                       cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                                 <  1:4774573/150‑1 (MQ=255)

                                                                            tCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCtt                                                                            >  2:842192/1‑150 (MQ=255)

                                                                             cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                                                                           >  2:3830490/1‑150 (MQ=255)

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

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

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

                                                                                           tCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATgg                                                             >  1:1753812/1‑150 (MQ=255)

                                                                                             gggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCg                                                           >  2:2683971/1‑150 (MQ=255)

                                                                                             gggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCg                                                           >  1:3742253/1‑150 (MQ=255)

                                                                                                tGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTac                                                        <  2:2576147/150‑1 (MQ=255)

                                                                                                 ggCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTaca                                                       <  1:1248112/150‑1 (MQ=255)

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

                                                                                                    cGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAcct                                                    >  1:604768/1‑150 (MQ=255)

                                                                                                        aCCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCtt                                                <  1:800229/150‑1 (MQ=255)

                                                                                                        aCCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCtt                                                <  2:2207175/150‑1 (MQ=255)

                                                                                                         ccAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTg                                               >  1:828802/1‑150 (MQ=255)

                                                                                                          cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                              >  2:1052057/1‑150 (MQ=255)

                                                                                                          cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                              >  2:3667576/1‑150 (MQ=255)

                                                                                                          cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                              <  1:3188332/150‑1 (MQ=255)

                                                                                                          cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                              >  1:3071266/1‑150 (MQ=255)

                                                                                                           aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTTGGCTCGCACCATGGCGTACACCTCCTTGAg                                             >  2:2878929/1‑150 (MQ=255)

                                                                                                           aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                             <  2:2784337/150‑1 (MQ=255)

                                                                                                           aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                             <  2:790338/150‑1 (MQ=255)

                                                                                                           aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                             <  1:1748551/150‑1 (MQ=255)

                                                                                                           aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                             <  2:5078093/150‑1 (MQ=255)

                                                                                                            aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                       >  1:4260704/1‑59 (MQ=255)

                                                                                                            aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                       <  2:4260704/59‑1 (MQ=255)

                                                                                                            aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAgg                                            >  1:2139912/1‑150 (MQ=255)

                                                                                                            aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAgg                                            <  2:4580501/150‑1 (MQ=255)

                                                                                                             tGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGt                                           <  1:4579391/150‑1 (MQ=255)

                                                                                                              gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                                          >  1:2641545/1‑150 (MQ=255)

                                                                                                              gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                                          >  2:50998/1‑150 (MQ=255)

                                                                                                                     cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                                   >  2:4505312/1‑150 (MQ=255)

                                                                                                                     cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                                   >  1:1082863/1‑150 (MQ=255)

                                                                                                                     cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                                   <  1:960927/150‑1 (MQ=255)

                                                                                                                      ccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCgg                                  <  2:3234732/150‑1 (MQ=255)

                                                                                                                      ccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCgg                                  >  1:2277805/1‑150 (MQ=255)

                                                                                                                      ccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCgg                                  <  1:5095134/150‑1 (MQ=255)

                                                                                                                          gCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCaa                              >  2:3437363/1‑150 (MQ=255)

                                                                                                                            tGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCaaaa                            <  2:1176475/150‑1 (MQ=255)

                                                                                                                              aCGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAac                          >  1:4191374/1‑150 (MQ=255)

                                                                                                                              aCGGCAACCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAac                          >  1:3095269/1‑150 (MQ=255)

                                                                                                                                ggCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAacac                        >  1:3438380/1‑150 (MQ=255)

                                                                                                                                 gCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACAcc                       >  1:4530063/1‑150 (MQ=255)

                                                                                                                                   aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTg                     >  2:905156/1‑150 (MQ=255)

                                                                                                                                   aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTg                     >  1:4223546/1‑150 (MQ=255)

                                                                                                                                   aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGACAAAAACACCTg                     >  2:4025418/1‑150 (MQ=255)

                                                                                                                                    tCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGc                    >  2:473516/1‑150 (MQ=255)

                                                                                                                                      cGGGACAACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGAACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGccc                  >  2:1776081/1‑150 (MQ=255)

                                                                                                                                       gggACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGcccc                 >  2:2261349/1‑150 (MQ=255)

                                                                                                                                             aCGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACg           >  2:3252087/1‑150 (MQ=255)

                                                                                                                                             aCGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTGGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACg           >  2:2101884/1‑150 (MQ=255)

                                                                                                                                                cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACTTTCGGCTCGCACCATGGCGTACACCTCCTCGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg        <  1:4592070/150‑1 (MQ=255)

                                                                                                                                                cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg        <  2:2955152/150‑1 (MQ=255)

                                                                                                                                                cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg        <  2:5077800/150‑1 (MQ=255)

                                                                                                                                                 ccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTgg       <  2:1457912/150‑1 (MQ=255)

                                                                                                                                                  cTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTggg      >  2:2735012/1‑150 (MQ=255)

                                                                                                                                                  cTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCGTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTggg      >  2:2735119/1‑150 (MQ=255)

ATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTGGGGATC  >  NC_002947/3803202‑3803500

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