BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	835,750	C→A	5.2%	S111I (AGC→ATC)	PP_0718 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	835,750	0	C	A	5.2%	371.5 / 4.7	115	S111I (AGC→ATC)	PP_0718	sulfate transporter
Reads supporting (aligned to +/- strand): ref base C (57/52); new base A (3/3); total (60/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 115 total aligned reads displayed.

GCAACTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGC  >  NC_002947/835610‑835898

gcgaCTGGCTTGTGAATACGGTCGGCGCGATGCCTATCATGTTCCCATAGCGGATGGGCGTGCACACCGACCCGGTCATAAGTCCCAACCGTAGCAGCCCGGCCGTCACTCAGCCCGCCCCGGCTACCAGCGCCATCATGCtggcgatgg                                                                                                                                             <  2:93900/147‑1 (MQ=255)

gCAACTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATTATGCtggcgatgg                                                                                                                                             <  2:601638/150‑1 (MQ=255)

gCAACTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatgg                                                                                                                                             <  1:2150741/150‑1 (MQ=255)

  aaCTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatggcg                                                                                                                                           >  2:330775/1‑150 (MQ=255)

    cTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATCGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCAGGGCAGCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTggcgatggcgat                                                                                                                                         >  1:734322/1‑150 (MQ=255)

                   ggTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTc                                                                                                                          <  2:470913/150‑1 (MQ=255)

                       aGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTg                                                                                                                      <  2:75272/150‑1 (MQ=255)

                        gggCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTgg                                                                                                                     >  2:1063145/1‑150 (MQ=255)

                        gggCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTgg                                                                                                                     <  2:214144/150‑1 (MQ=255)

                            gATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTggcggc                                                                                                                 <  2:802991/150‑1 (MQ=255)

                             aTGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTTGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGgcggcg                                                                                                                >  1:1307431/1‑150 (MQ=255)

                              tGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTggcg                                                                                                                   <  2:2225490/146‑1 (MQ=255)

                                ccGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTAt                                                                                                             <  1:330775/150‑1 (MQ=255)

                                 cGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGACCGCTGCGGGTCGCTGGCGGCGTAtt                                                                                                            >  1:705468/1‑150 (MQ=255)

                                    tCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGca                                                                                                         <  2:2172440/150‑1 (MQ=255)

                                     cATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGcac                                                                                                        <  1:1945080/150‑1 (MQ=255)

                                          aCCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGcaccacca                                                                                                   <  2:1961999/150‑1 (MQ=255)

                                           cccATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatggcga                                                                                                                                          <  2:1057283/110‑1 (MQ=255)

                                           cccATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatggcga                                                                                                                                          >  1:1057283/1‑110 (MQ=255)

                                           cccATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGcaccaccac                                                                                                  >  1:1068442/1‑150 (MQ=255)

                                            ccATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatggc                                                                                                                                            <  2:231171/107‑1 (MQ=255)

                                            ccATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcgatggc                                                                                                                                            >  1:231171/1‑107 (MQ=255)

                                             cATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTAt                                                                                                             >  1:1293059/1‑137 (MQ=255)

                                             cATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAg                                                                                                >  2:1367214/1‑150 (MQ=255)

                                              aTAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCGGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGc                                                                                               <  1:2082438/150‑1 (MQ=255)

                                                aGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCa                                                                                             >  2:2194476/1‑150 (MQ=255)

                                                   ggATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGGGCAGCAGCCCGGCAATCACACATAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGa                                                                                          >  2:121253/1‑150 (MQ=255)

                                                    gATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCg                                                                                                                              >  1:1647060/1‑113 (MQ=255)

                                                    gATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCg                                                                                                                              <  2:1647060/113‑1 (MQ=255)

                                                    gATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGAt                                                                                         >  1:1899260/1‑150 (MQ=255)

                                                     aTGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGGAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATa                                                                                        <  1:2016046/150‑1 (MQ=255)

                                                        ggCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGgc                                                                                     >  1:926244/1‑150 (MQ=255)

                                                         gCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGgcg                                                                                    >  2:87965/1‑150 (MQ=255)

                                                         gCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGgcg                                                                                    >  1:229688/1‑150 (MQ=255)

                                                          cTTGGCCAGCAACCCGGTGATGAAGCCCACATGGAACACATCGGCATATACACAGAACGCTCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGgcgc                                                                                   <  2:734322/150‑1 (MQ=255)

                                                             ggACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggc                                                                                                                                                  <  1:915789/84‑1 (MQ=255)

                                                             ggACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggc                                                                                                                                                  >  2:915789/1‑84 (MQ=255)

                                                               aCAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCa                                                                              >  2:1964882/1‑150 (MQ=255)

                                                               aCAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCa                                                                              >  1:752123/1‑150 (MQ=255)

                                                                 aGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTggcg                                                                                                                   >  2:723694/1‑111 (MQ=255)

                                                                 aGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTggcg                                                                                                                   <  1:723694/111‑1 (MQ=255)

                                                                 aGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTTcaccac                                                                                                     <  1:450577/125‑1 (MQ=255)

                                                                 aGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGcaccac                                                                                                     >  2:450577/1‑125 (MQ=255)

                                                                    aaCTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcga                                                                                                                                                <  1:551793/79‑1 (MQ=255)

                                                                    aaCTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCtggcga                                                                                                                                                >  2:551793/1‑79 (MQ=255)

                                                                    aaCTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCACCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGcc                                                                         >  2:1152927/1‑150 (MQ=255)

                                                                      cTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGa                                                                       <  2:752123/150‑1 (MQ=255)

                                                                      cTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGa                                                                       <  2:713005/150‑1 (MQ=255)

                                                                      cTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGa                                                                       >  2:1242342/1‑150 (MQ=255)

                                                                      cTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGa                                                                       >  1:162829/1‑150 (MQ=255)

                                                                       tCGGTGATAAAGACCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGTCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAg                                                                      >  1:51953/1‑150 (MQ=255)

                                                                         ggTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAg                                                                                       <  1:54826/131‑1 (MQ=255)

                                                                         ggTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAg                                                                                       >  2:54826/1‑131 (MQ=255)

                                                                           tGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCg                                                                  <  2:926244/150‑1 (MQ=255)

                                                                            gagaacgCCCAGGCGCAGCAGCCCGGCCATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCgg                                                                 <  2:705468/144‑1 (MQ=255)

                                                                            gATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCgg                                                                 <  1:848962/150‑1 (MQ=255)

                                                                            gATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCgg                                                                 >  2:767287/1‑150 (MQ=255)

                                                                             aTAAAGCCCAGTCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCTCTGCGCGTCGCTGGCGGCGTATTCCACCACCACAGCCAGTATAGGCGCGGCCAGCGCCGAGTCCggg                                                                <  2:387841/150‑1 (MQ=255)

                                                                               aaaGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTTGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCAAGGCCGGGcc                                                              >  2:1811377/1‑150 (MQ=255)

                                                                                aaGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGccc                                                             >  1:2092350/1‑150 (MQ=255)

                                                                                 aGCCGATGCGCAGCAGCCCGCCAATCACACAGAACGCGCCGGCGACCAGCGCCATTAGGCGGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTCCACCCCCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCa                                                            <  2:859092/150‑1 (MQ=255)

                                                                                  gCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGTGACCAGCGCCATCATGCTGGCGCTGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAg                                                           >  1:132848/1‑150 (MQ=255)

                                                                                        gcgcAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAg                                                     >  2:1253123/1‑150 (MQ=255)

                                                                                          gcagcaGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGa                                                   <  1:1367214/150‑1 (MQ=255)

                                                                                                 cccGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCt                                            <  1:1152927/150‑1 (MQ=255)

                                                                                                   cGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGcgcg                                                                                  <  2:593259/110‑1 (MQ=255)

                                                                                                   cGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGcgcg                                                                                  >  1:593259/1‑110 (MQ=255)

                                                                                                   cGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTgg                                          <  2:1380894/150‑1 (MQ=255)

                                                                                                   cGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTgg                                          >  1:1256355/1‑150 (MQ=255)

                                                                                                   cGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTgg                                          <  1:1242342/150‑1 (MQ=255)

                                                                                                       aaTCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGAcc                                      <  2:1068442/150‑1 (MQ=255)

                                                                                                         tCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGa                                    >  1:158082/1‑150 (MQ=255)

                                                                                                          cacacaGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGaa                                   <  2:542698/150‑1 (MQ=255)

                                                                                                              caGAACGCCCCGGCGACCAGCGCCATCATGCTGGCTATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAgg                               <  1:1407729/150‑1 (MQ=255)

                                                                                                              caGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAgg                               >  2:1260831/1‑150 (MQ=255)

                                                                                                              caGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAgg                               <  1:87965/150‑1 (MQ=255)

                                                                                                                    gCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAg                         <  2:162829/150‑1 (MQ=255)

                                                                                                                    gCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAg                         >  1:941003/1‑150 (MQ=255)

                                                                                                                     ccccGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCgg                                                                                 <  1:1439088/93‑1 (MQ=255)

                                                                                                                     ccccGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCgg                                                                                 >  2:1439088/1‑93 (MQ=255)

                                                                                                                     ccccGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCa                                                      >  1:1047894/1‑120 (MQ=255)

                                                                                                                     ccccGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCa                                                      <  2:1047894/120‑1 (MQ=255)

                                                                                                                        cGGCGACCAGCGCCATCATGCTGGCGATGGCGTTCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGccc                     >  1:656206/1‑149 (MQ=255)

                                                                                                                             gccatcgcCAGCATGATGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGc                             >  1:1177414/6‑137 (MQ=255)

                                                                                                                             gccATCGCCAGCATGATGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGc                             <  2:1177414/136‑1 (MQ=255)

                                                                                                                              ccatcgcCAGCATGATGGCGCTGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCt                                            <  2:1773368/117‑1 (MQ=255)

                                                                                                                              ccATCGCCAGCATGATGGCGCTGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGa                                                                       >  1:1762825/1‑94 (MQ=38)

                                                                                                                              ccATCGCCAGCATGATGGCGCTGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCt                                            >  1:1773368/1‑121 (MQ=255)

                                                                                                                              ccAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCTGGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAggg               >  1:404010/1‑150 (MQ=255)

                                                                                                                              ccAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAggg               >  1:1940992/1‑150 (MQ=255)

                                                                                                                              ccAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAggg               >  2:2150369/1‑150 (MQ=255)

                                                                                                                               cAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTgg                                          >  1:1614345/1‑122 (MQ=255)

                                                                                                                               cAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTCGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTgg                                          <  2:1614345/122‑1 (MQ=255)

                                                                                                                                    ccATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGgatga         >  1:1127448/1‑150 (MQ=255)

                                                                                                                                     catcatGCTGGCGATGGCGATCGCCCGCTGCGGGTCGTCGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGgatgat        <  1:612530/150‑1 (MQ=255)

                                                                                                                                     catcatGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGgatgat        >  1:859542/1‑150 (MQ=255)

                                                                                                                                     catcatGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGgatgat        >  1:309057/1‑150 (MQ=255)

                                                                                                                                       tcatGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATgg      <  2:2092350/150‑1 (MQ=255)

                                                                                                                                         atGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTa    <  2:1127448/150‑1 (MQ=255)

                                                                                                                                           gCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGc  <  2:1467189/150‑1 (MQ=255)

GCAACTGGCTGATCAATACGGTCAGGGCGATGCCGTTCATGTACCCATAGCGGATGGGCTTGGACAGCAACTCGGTGATAAAGCCCAGGCGCAGCAGCCCGGCAATCACACAGAACGCCCCGGCGACCAGCGCCATCATGCTGGCGATGGCGATCGCCCGCTGCGGGTCGCTGGCGGCGTATTGCACCACCACAGCCAGGATAGGCGCGGCCAGCGCCGAGTCCGGGCCCAGCACCAGGATACGGCTGGGACCGAACAGGGCATAGGCCAGCAGGGGGATGATGGTAGC  >  NC_002947/835610‑835898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`