BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,816,587	C→G	8.0%	A114G (GCC→GGC)	PP_5092 →	NLP/P60 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,816,587	0	C	G	8.0%	348.1 / 18.6	112	A114G (GCC→GGC)	PP_5092	NLP/P60 protein
Reads supporting (aligned to +/- strand): ref base C (52/51); new base G (5/4); total (57/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 = 4.70e-01

Only 100 of 112 total aligned reads displayed.

TGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816441‑5816735

tGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTACGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGa                                                                                                                                                   <  2:1656023/150‑1 (MQ=255)

tGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCCTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGa                                                                                                                                                   <  2:1649870/150‑1 (MQ=255)

  tCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACg                                                                                                                                                 >  1:322812/1‑150 (MQ=255)

        ccATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGAcc                                                                                                                                           >  1:604816/1‑150 (MQ=255)

                 cTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCa                                                                                                                                  >  2:117480/1‑150 (MQ=255)

                    agcaaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCt                                                                                                                               >  2:464817/1‑150 (MQ=255)

                             gcGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCa                                                                                                                      <  2:828235/150‑1 (MQ=255)

                               gATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATg                                                                                                                    <  2:322812/150‑1 (MQ=255)

                                aTGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATgg                                                                                                                   >  2:1803679/1‑150 (MQ=255)

                                      tGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTTACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             >  2:1376555/1‑150 (MQ=255)

                                      tGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             <  1:117480/150‑1 (MQ=255)

                                      tGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             >  2:101395/1‑150 (MQ=255)

                                      tGAGCCGAGCGGTAAATGTGCTCGGTACTCCGTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGGCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             >  2:1970871/1‑150 (MQ=255)

                                        aGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCgg                                                                                                              <  1:1028455/147‑1 (MQ=255)

                                        aGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCgg                                                                                                              >  2:1028455/1‑147 (MQ=255)

                                           cGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCAc                                                                                                        <  2:2186449/150‑1 (MQ=255)

                                           cGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCCGGGCCat                                                                                                        <  2:1030502/150‑2 (MQ=255)

                                              gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCTTGTCGGTGGGCCACggg                                                                                                     <  2:1920035/150‑1 (MQ=255)

                                              gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  2:1440959/150‑1 (MQ=255)

                                              gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  1:1376555/150‑1 (MQ=255)

                                              gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  2:604816/150‑1 (MQ=255)

                                              gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCCCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  1:2162847/150‑1 (MQ=255)

                                                  tAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCaa                                                                                                                           >  2:2179277/1‑124 (MQ=255)

                                                  tAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCaa                                                                                                                           <  1:2179277/124‑1 (MQ=255)

                                                    aaTGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAg                                                                                               <  2:204466/150‑1 (MQ=255)

                                                         gCTCGGGAGTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGc                                                                                          <  1:101395/150‑1 (MQ=255)

                                                          cTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGcc                                                                                         >  2:707722/1‑150 (MQ=255)

                                                           tCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCa                                                                                        >  2:911742/1‑150 (MQ=255)

                                                           tCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCa                                                                                        >  2:155673/1‑150 (MQ=255)

                                                                  tCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTg                                                                                  <  1:153608/149‑1 (MQ=255)

                                                                  tCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTg                                                                                  >  2:153608/1‑149 (MQ=255)

                                                                      tATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTg                                                                             <  2:889460/150‑1 (MQ=255)

                                                                       aTGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTACAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGa                                                                            <  2:808942/150‑1 (MQ=255)

                                                                        tGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCAAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGaa                                                                           >  1:704921/1‑150 (MQ=255)

                                                                          tttGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTGCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGc                                                                         >  2:1304038/1‑150 (MQ=255)

                                                                               ggcggcAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                    <  1:676272/150‑1 (MQ=255)

                                                                                gcggcAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCg                                                                   >  2:1658160/1‑150 (MQ=255)

                                                                                gcggcAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCg                                                                   >  2:2190900/1‑150 (MQ=255)

                                                                                gcggcAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATACAGTGGTCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCg                                                                   >  1:145667/1‑150 (MQ=255)

                                                                                  ggcAACAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 <  1:440270/150‑1 (MQ=255)

                                                                                    cagcagCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCt                                                               <  1:2412136/150‑1 (MQ=255)

                                                                                       cagcCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAt                                                            >  2:1635365/1‑150 (MQ=255)

                                                                                        agcCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtc                                                           >  1:2449194/1‑150 (MQ=255)

                                                                                          cccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtctt                                                         <  1:155673/150‑1 (MQ=255)

                                                                                          cccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtctt                                                         >  1:2368917/1‑150 (MQ=255)

                                                                                           ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttc                                                        >  2:1554066/1‑150 (MQ=255)

                                                                                           ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttc                                                        >  1:73388/1‑150 (MQ=255)

                                                                                           ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttc                                                        >  2:2001353/1‑150 (MQ=255)

                                                                                            cGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttct                                                       >  1:882869/1‑150 (MQ=255)

                                                                                            cGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttct                                                       >  1:815806/1‑150 (MQ=255)

                                                                                            cGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttct                                                       >  2:2148287/1‑150 (MQ=255)

                                                                                             gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                      <  1:911742/150‑1 (MQ=255)

                                                                                             gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                      <  2:2378424/150‑1 (MQ=255)

                                                                                              aaaaaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                     <  1:1658160/150‑1 (MQ=255)

                                                                                              aaaaaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                     <  2:582745/150‑1 (MQ=255)

                                                                                                  aaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACa                                                 <  1:1635365/150‑1 (MQ=255)

                                                                                                   aGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACAt                                                >  1:142254/1‑150 (MQ=255)

                                                                                                   aGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACAt                                                >  1:583504/1‑150 (MQ=255)

                                                                                                       tttGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAAt                                                                                                                          <  1:304654/72‑1 (MQ=255)

                                                                                                       tttGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAAt                                                                                                                          >  2:304654/1‑72 (MQ=255)

                                                                                                       tttGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCAGGTGACCTGATCTTCTTCAACATCAag                                            >  1:959672/1‑150 (MQ=255)

                                                                                                        ttGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAaga                                           >  2:1418469/1‑150 (MQ=255)

                                                                                                           aCTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGcc                                        >  2:183809/1‑150 (MQ=255)

                                                                                                             tGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACAGCAAGAGCCGt                                      >  1:2440671/1‑150 (MQ=255)

                                                                                                                  cGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGt                                 <  1:1304038/150‑1 (MQ=255)

                                                                                                                     gCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTAAACATCAAGAGCCGTCGGGTCaa                              >  2:1424901/1‑150 (MQ=255)

                                                                                                                       tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                <  1:1247436/98‑1 (MQ=255)

                                                                                                                       tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                >  2:1247436/1‑98 (MQ=255)

                                                                                                                       tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATc                            <  2:2449194/150‑1 (MQ=255)

                                                                                                                         gTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCGCCTCCAAt                                                                                                                          >  1:213655/1‑54 (MQ=255)

                                                                                                                         gTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCGCCTCCAAt                                                                                                                          <  2:213655/54‑1 (MQ=255)

                                                                                                                          tCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCgg                                   >  2:1919153/1‑140 (MQ=255)

                                                                                                                          tCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCgg                                   <  1:1919153/140‑1 (MQ=255)

                                                                                                                          tCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACg                         <  2:2440671/150‑1 (MQ=255)

                                                                                                                              aTACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAAAGCCATGGCCCGGGGCCCCGGGGTGAAGGCGGCCAAGGGTGGTCTGAAGCCTGGCGACCTGCTCTTCTTCAACATCAAGAGCCGTCGGGGCACTCACGCAGc                     >  2:2150490/1‑150 (MQ=255)

                                                                                                                                aCGCCTTCAACGCCGTCGCCGACGACGACCTGCCGCGCACCTCTCATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCa                   <  2:145667/150‑1 (MQ=255)

                                                                                                                                aCGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCa                   >  1:1119169/1‑150 (MQ=255)

                                                                                                                                  gCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAAt                             >  1:2169217/1‑138 (MQ=255)

                                                                                                                                  gCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAAt                             <  2:2169217/138‑1 (MQ=255)

                                                                                                                                    cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                     <  2:948490/112‑1 (MQ=255)

                                                                                                                                    cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                     >  1:948490/1‑112 (MQ=255)

                                                                                                                                     ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAATGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc              >  2:1370095/1‑150 (MQ=255)

                                                                                                                                     ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCt                                                               >  1:1490509/1‑101 (MQ=255)

                                                                                                                                     ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc              <  2:882869/150‑1 (MQ=255)

                                                                                                                                       ggtcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCt            >  1:678035/4‑150 (MQ=255)

                                                                                                                                        aaCGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                    >  2:2074808/1‑93 (MQ=255)

                                                                                                                                        aaCGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                    <  1:2074808/93‑1 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                      >  2:284626/2‑106 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                      <  1:284626/105‑1 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCaa                                             <  1:977294/114‑1 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          >  2:833623/2‑150 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTG‑AGGTGGCCAAGGGTGATCTGGAGCCTGGCGACCTGAt                                                            >  1:1731670/2‑99 (MQ=255)

                                                                                                                                         tcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTG‑AGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAt                                                            <  2:1731670/98‑1 (MQ=255)

                                                                                                                                         acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  1:971292/150‑1 (MQ=255)

                                                                                                                                         acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  1:1370095/150‑1 (MQ=255)

                                                                                                                                         acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  1:1803679/150‑1 (MQ=255)

                                                                                                                                         acgacgTCGCCGACGTCGACCCGCCGCGCCCCTACAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg          <  2:1947657/150‑1 (MQ=255)

                                                                                                                                          cgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCTATCTACCTggg         <  1:2423099/150‑1 (MQ=255)

                                                                                                                                          cgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTggg         >  1:2045690/1‑150 (MQ=255)

                                                                                                                                                 gCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAc  >  1:208339/1‑150 (MQ=255)

TGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816441‑5816735

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