BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,106,909	A→C	5.4%	A55A (GCA→GCC) ‡	PP_3610 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,106,909	0	A	C	5.4%	875.9 / 27.3	260	A55A (GCA→GCC) ‡	PP_3610	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (125/121); new base C (7/7); total (132/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

Only 100 of 262 total aligned reads displayed.

AGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGCCCAGGC  >  NC_002947/4106761‑4107058

 gTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACt                                                                                                                                                     <  2:346616/150‑1 (MQ=255)

  tCCCTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTg                                                                                                                                                    <  2:2808931/150‑1 (MQ=255)

    ccTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTggg                                                                                                                                                  >  2:5248457/1‑150 (MQ=255)

     cTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGc                                                                                                                                                 <  1:4512013/150‑1 (MQ=255)

           gACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGttc                                                                                                                                           <  2:2107335/150‑3 (MQ=255)

               gTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCaaaaa                                                                                                                                       >  1:5352770/1‑150 (MQ=255)

                 ggCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCaaaaaaa                                                                                                                                     <  1:3968203/150‑1 (MQ=255)

                   ccACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                                                                                                   >  1:3581785/1‑150 (MQ=255)

                    cACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCCTTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTg                                                                                                                                  >  2:4790966/1‑150 (MQ=255)

                    cACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTg                                                                                                                                  >  2:877731/1‑150 (MQ=255)

                    cACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCGTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGCGCAGGCCCAGGCCCTGGGCCGGTCAAAAAAAGTg                                                                                                                                  >  2:3088294/1‑150 (MQ=255)

                       cTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTggcg                                                                                                                               >  1:1005860/1‑150 (MQ=255)

                       cTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTggcg                                                                                                                               >  1:5420800/1‑150 (MQ=255)

                       cTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCGGTCACAAAAAGTggcg                                                                                                                               >  2:4826696/1‑150 (MQ=255)

                          ccAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGcc                                                                                                                            >  1:1441776/1‑150 (MQ=255)

                           cAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCa                                                                                                                           >  1:575927/1‑150 (MQ=255)

                            aaCCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCAt                                                                                                                          <  2:5530936/150‑1 (MQ=255)

                             aCCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCAtt                                                                                                                         <  1:2232370/150‑1 (MQ=255)

                             aCCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCAtt                                                                                                                         >  1:901769/1‑150 (MQ=255)

                             aCCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCAtt                                                                                                                         <  1:403466/150‑1 (MQ=255)

                              cccGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTg                                                                                                                        <  1:4458994/150‑1 (MQ=255)

                              cccGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTg                                                                                                                        <  2:1249548/150‑1 (MQ=255)

                               ccGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTgg                                                                                                                       >  1:1392284/1‑150 (MQ=255)

                                 gCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGcc                                                                                                                     <  2:3279082/150‑1 (MQ=255)

                                 gCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGcc                                                                                                                     >  2:4876760/1‑150 (MQ=255)

                                 ccAACCTCGCCGAATGCCAGGCGCTCGACCAAGCCGCCTTGCTCAACGCCTAGACAGGGCAACTCCCAGCAGCATTCGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGCCCCTGTAAAAAAAAGTGGCGGCCATTGGcc                                                                                                                     <  2:4058958/149‑1 (MQ=255)

                                  cAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCt                                                                                                                    >  1:5150683/1‑150 (MQ=255)

                                    aCCTGGCCGAATGCCCGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTgg                                                                                                                  <  2:591561/150‑1 (MQ=255)

                                      cTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCa                                                                                                                >  1:2805289/1‑150 (MQ=255)

                                        ggCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACt                                                                                                              <  1:2073923/150‑1 (MQ=255)

                                         gCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAATGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTg                                                                                                             >  2:4476337/1‑150 (MQ=255)

                                          ccGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTgg                                                                                                            <  1:28666/150‑1 (MQ=255)

                                          ccGAATGCCTGGCGGTGGACCAAGCCGCCTTCGTCAAAGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCACGCTGCCAGTGCGGAGCAGGCCCACGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTgg                                                                                                            <  2:3625392/150‑1 (MQ=255)

                                           cGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTggg                                                                                                           >  1:3828992/1‑150 (MQ=255)

                                            gAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGc                                                                                                          <  1:3183271/150‑1 (MQ=255)

                                              aTGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCaa                                                                                                        <  2:3837896/150‑1 (MQ=255)

                                               tGCCTGGCGGTGGACCAAGCCGCCTTGGTTAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAg                                                                                                       >  1:4442437/1‑150 (MQ=255)

                                               tGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAg                                                                                                       >  2:80249/1‑150 (MQ=255)

                                                   tGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGc                                                                                                   <  1:5563540/150‑1 (MQ=255)

                                                    ggCGGTGGACCCAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCt                                                                                                  <  2:442936/150‑1 (MQ=255)

                                                     gCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTg                                                                                                 <  2:5446024/150‑1 (MQ=255)

                                                       ggTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGACACTGGGCAAGTGGCTGGa                                                                                               <  2:862512/150‑1 (MQ=255)

                                                         ttgACCAAGGCGGCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGc                                                                                             >  2:2401989/3‑150 (MQ=255)

                                                         tGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGc                                                                                             >  1:1040418/1‑150 (MQ=255)

                                                         tGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGc                                                                                             <  1:4298978/150‑1 (MQ=255)

                                                              cAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAgcc                                                                                        >  1:2703014/1‑150 (MQ=255)

                                                                 gccgccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAgccgcc                                                                                     <  1:708853/150‑1 (MQ=255)

                                                                       ttGGTCAACGCCTTGACAGGGCAACTACCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCg                                                                               >  2:145080/1‑150 (MQ=255)

                                                                        tGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGt                                                                              <  1:1478097/150‑1 (MQ=255)

                                                                           tCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTgc                                                                           >  1:2934182/1‑150 (MQ=255)

                                                                             aaCGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTgcgc                                                                         >  2:1157778/1‑150 (MQ=255)

                                                                             aaCGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTgcgc                                                                         >  2:1595221/1‑150 (MQ=255)

                                                                               cGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCa                                                                       >  2:2986241/1‑150 (MQ=255)

                                                                                  cTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCa                                                                    >  2:3546398/1‑150 (MQ=255)

                                                                                     gACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAgt                                                                 <  2:4609452/150‑1 (MQ=255)

                                                                                        aGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGtt                                                              >  1:631052/1‑150 (MQ=255)

                                                                                          ggCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCg                                                            >  2:2546215/1‑150 (MQ=255)

                                                                                           gCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGa                                                           <  1:3808787/150‑1 (MQ=255)

                                                                                           gCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGa                                                           <  1:4221366/150‑1 (MQ=255)

                                                                                            cAACTGCCAGCAGCATTGGGAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAg                                                          <  2:1588342/150‑1 (MQ=255)

                                                                                             aaCTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGc                                                         >  1:5406427/1‑150 (MQ=255)

                                                                                                tGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAAt                                                      >  1:2122331/1‑150 (MQ=255)

                                                                                                 gCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAAtg                                                     >  2:1972912/1‑150 (MQ=255)

                                                                                                  ccAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAAtgc                                                    >  1:2775194/1‑150 (MQ=255)

                                                                                                       agcaTTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTgcac                                               >  1:3715198/1‑150 (MQ=255)

                                                                                                          aTTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTgcacgca                                            >  1:2378451/1‑150 (MQ=255)

                                                                                                          aTTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTgcacgca                                            >  1:936254/1‑150 (MQ=255)

                                                                                                                agagCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGTCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATc                                      <  2:664130/150‑1 (MQ=255)

                                                                                                                 gagCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATcc                                     <  2:5406427/150‑1 (MQ=255)

                                                                                                                    cACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGAt                                  <  2:1403507/150‑1 (MQ=255)

                                                                                                                     aCCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATa                                 <  1:1282687/150‑1 (MQ=255)

                                                                                                                       cTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATAcc                               <  1:1795284/150‑1 (MQ=255)

                                                                                                                        tCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCg                              >  1:1976341/1‑150 (MQ=255)

                                                                                                                         cGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGt                             <  2:437990/150‑1 (MQ=255)

                                                                                                                           cTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTgc                           <  2:127675/150‑1 (MQ=255)

                                                                                                                           cTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTgc                           <  1:5482661/150‑1 (MQ=255)

                                                                                                                            tGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTgcg                          <  1:4599197/150‑1 (MQ=255)

                                                                                                                              ccAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCTTTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCa                        >  2:4550567/1‑150 (MQ=255)

                                                                                                                               cAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAg                       >  1:2450503/1‑150 (MQ=255)

                                                                                                                               cAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAg                       >  2:3968615/1‑150 (MQ=255)

                                                                                                                               cAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAg                       >  1:5274813/1‑150 (MQ=255)

                                                                                                                                 gTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCt                     >  2:2831276/1‑150 (MQ=255)

                                                                                                                                  tGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTg                    <  2:2068190/150‑1 (MQ=255)

                                                                                                                                   gCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTgg                   >  2:1383923/1‑150 (MQ=255)

                                                                                                                                    cGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTggg                  >  2:3195659/1‑150 (MQ=255)

                                                                                                                                     ggAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCCGCTGggc                 <  1:483474/150‑1 (MQ=255)

                                                                                                                                     ggAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGggc                 >  2:2876629/1‑150 (MQ=255)

                                                                                                                                        acAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGggcggc              >  1:5591499/2‑150 (MQ=255)

                                                                                                                                         cAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCa             >  2:4164560/1‑150 (MQ=255)

                                                                                                                                          aGGCCCAGTGCCTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCAt            >  1:166241/1‑150 (MQ=255)

                                                                                                                                          aGGCCCAGTGCCTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGAGGCAt            >  2:3824965/1‑150 (MQ=255)

                                                                                                                                             cccAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCg         <  1:3589689/150‑1 (MQ=255)

                                                                                                                                              ccAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGc        >  1:4924815/1‑150 (MQ=255)

                                                                                                                                              ccAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGc        >  2:3124681/1‑150 (MQ=255)

                                                                                                                                               cAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGcc       >  2:2755708/1‑150 (MQ=255)

                                                                                                                                                aGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGccc      >  1:120943/1‑150 (MQ=255)

                                                                                                                                                aGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGccc      >  2:636776/1‑150 (MQ=255)

                                                                                                                                                 ggCACTGGGCCTGTCAAACAAAGGGGCGGCCATTGGCCTGGCACTGGGCACGGGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGCCCa     <  2:3581790/150‑1 (MQ=255)

                                                                                                                                                 ggCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGCCCa     <  2:5150683/150‑1 (MQ=255)

                                                                                                                                                 ggCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGCCCa     <  2:2268796/150‑1 (MQ=255)

AGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAACCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGAGTCGTTGCGCCAGCAAGTGTTCGAGCAATGCTGCACGCACCCATCCGATACCGTGCGCAGCTGGGCGGCATTCGCCCAGGC  >  NC_002947/4106761‑4107058

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