BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,254,654	A→C	5.2%	K428Q (AAG→CAG)	ureC →	urease subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,254,654	0	A	C	5.2%	1382.7 / 44.7	403	K428Q (AAG→CAG)	ureC	urease subunit alpha
Reads supporting (aligned to +/- strand): ref base A (183/199); new base C (16/5); total (199/204)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 410 total aligned reads displayed.

GCACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTTC  >  NC_002947/3254506‑3254803

gcACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCaa                                                                                                                                                      >  1:5159826/1‑150 (MQ=255)

 cACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAg                                                                                                                                                     >  1:4086556/1‑150 (MQ=255)

     aGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTgg                                                                                                                                                 <  1:6398018/150‑1 (MQ=255)

      gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                                <  2:2979953/150‑1 (MQ=255)

      gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                                <  1:302609/150‑1 (MQ=255)

       aTGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGcc                                                                                                                                               <  2:4534058/150‑1 (MQ=255)

          aaGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAc                                                                                                                                            <  1:6949233/150‑1 (MQ=255)

             cTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGACCtc                                                                                                                                         <  2:8667936/150‑2 (MQ=255)

             cTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTg                                                                                                                                         >  1:7897051/1‑150 (MQ=255)

                  tCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCt                                                                                                                                    <  2:5595926/150‑1 (MQ=255)

                  tCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCt                                                                                                                                    <  1:7453920/150‑1 (MQ=255)

                    gcgGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgt                                                                                                                                  <  1:8790433/150‑1 (MQ=255)

                    gcgGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgt                                                                                                                                  >  1:5473416/1‑150 (MQ=255)

                    gcgGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgt                                                                                                                                  >  1:1356706/1‑150 (MQ=255)

                         ccGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCg                                                                                                                             >  1:1581003/1‑150 (MQ=255)

                            cTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCa                                                                                                                          <  1:8501349/150‑1 (MQ=255)

                            cTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGACa                                                                                                                          <  2:1630075/150‑1 (MQ=255)

                             tGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAg                                                                                                                         >  2:2006152/1‑150 (MQ=255)

                               gCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCg                                                                                                                       <  2:1857315/150‑1 (MQ=255)

                                 aCCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGtg                                                                                                                     >  2:6056507/1‑149 (MQ=255)

                                  ccGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTc                                                                                                                    >  2:1745315/1‑150 (MQ=255)

                                   cGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCt                                                                                                                   >  1:6648692/1‑150 (MQ=255)

                                     gACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCttt                                                                                                                 <  1:3864007/150‑1 (MQ=255)

                                     gACACCCCGGACCGCGACACCTTCCGGGTCAAGCGCCCCATCGCCAAGGACCCCAGCCACCCGGCACTCCCCCATGGCATCGGTCACGACGGGGGCTCGGAGGAGGTCGGCAAGCGGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCttt                                                                                                                 <  2:4014515/150‑1 (MQ=255)

                                      acacCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTg                                                                                                                <  2:5170979/150‑1 (MQ=255)

                                          cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGt                                                                                                            >  2:625359/1‑150 (MQ=255)

                                          cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGt                                                                                                            >  1:6499647/1‑150 (MQ=255)

                                          cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGt                                                                                                            >  2:6490904/1‑150 (MQ=255)

                                            cGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCa                                                                                                          >  2:8093769/1‑150 (MQ=255)

                                                   cGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGc                                                                                                   >  2:7558109/1‑150 (MQ=255)

                                                    gACAACTTCCGGGTCAAGCGCTACATCGCCAAGTCCACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGcc                                                                                                  >  2:4393181/1‑150 (MQ=255)

                                                      cAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCt                                                                                                <  2:8021848/150‑1 (MQ=255)

                                                       aaCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTg                                                                                               <  2:5268933/150‑1 (MQ=255)

                                                       aaCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTg                                                                                               >  2:7869728/1‑150 (MQ=255)

                                                          ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                            >  1:2056962/1‑150 (MQ=255)

                                                             cGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTc                                                                                         <  2:473187/150‑1 (MQ=255)

                                                                   aaGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTTGAGCACCAGGGc                                                                                   <  1:5329473/150‑1 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  >  1:184536/1‑150 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  >  2:1530213/1‑150 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  >  2:1242284/1‑150 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  <  1:7105106/150‑1 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  <  1:5810399/150‑1 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  >  1:4744651/1‑150 (MQ=255)

                                                                    aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                  <  1:1571491/150‑1 (MQ=255)

                                                                     gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCgg                                                                                 >  1:8221500/1‑150 (MQ=255)

                                                                      cgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCggg                                                                                <  1:2239442/150‑1 (MQ=255)

                                                                         tACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatg                                                                             >  1:7304589/1‑150 (MQ=255)

                                                                         ggctgcGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatg                                                                             >  2:1383661/6‑150 (MQ=255)

                                                                          aCATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatga                                                                            <  1:7260422/150‑1 (MQ=255)

                                                                           cATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatgat                                                                           <  2:2418085/150‑1 (MQ=255)

                                                                            ctCGCCACGTACACCCTCCACCCGGCACTCCCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                          <  1:6168780/149‑1 (MQ=255)

                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                          >  1:2297780/1‑150 (MQ=255)

                                                                                 cAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCAc                                                                     >  1:7656043/1‑150 (MQ=255)

                                                                                     tACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACcgcg                                                                 >  2:5736366/1‑150 (MQ=255)

                                                                                        acCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGccc                                                              <  1:2866968/150‑1 (MQ=255)

                                                                                        acCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGccc                                                              >  2:752065/1‑150 (MQ=255)

                                                                                         ccATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCa                                                             >  1:8084945/1‑150 (MQ=255)

                                                                                           aTCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATg                                                           <  1:457251/150‑1 (MQ=255)

                                                                                           aTCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATg                                                           >  2:2836937/1‑150 (MQ=255)

                                                                                            tCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                          <  1:2509510/150‑1 (MQ=255)

                                                                                            tCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                          >  2:4185792/1‑150 (MQ=255)

                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGc                                                        <  1:1455093/150‑1 (MQ=255)

                                                                                               aCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCg                                                       >  1:4979843/1‑150 (MQ=255)

                                                                                               aCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCg                                                       <  1:368016/150‑1 (MQ=255)

                                                                                                 ccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCGGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGAc                                                     <  2:3250651/150‑1 (MQ=255)

                                                                                                      aCTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCaa                                                >  1:3484730/1‑150 (MQ=255)

                                                                                                       cTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAAc                                               <  2:1090503/150‑1 (MQ=255)

                                                                                                        tCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACg                                              >  2:3440240/1‑150 (MQ=255)

                                                                                                        tCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACg                                              >  1:753799/1‑150 (MQ=255)

                                                                                                          aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGc                                            <  1:5915278/150‑1 (MQ=255)

                                                                                                           cccATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCt                                           <  1:2654026/150‑1 (MQ=255)

                                                                                                            ccATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTc                                          >  1:6806436/1‑150 (MQ=255)

                                                                                                            ccATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTc                                          >  2:2280255/1‑150 (MQ=255)

                                                                                                            ccATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTc                                          <  2:2161126/150‑1 (MQ=255)

                                                                                                             cATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACTGCTcc                                         <  2:7270367/150‑1 (MQ=255)

                                                                                                               tGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCAt                                       >  2:3719567/1‑150 (MQ=255)

                                                                                                                       gTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATacccaccc                               <  2:6648692/150‑1 (MQ=255)

                                                                                                                           cGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCa                           <  1:776779/150‑1 (MQ=255)

                                                                                                                             aaGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGc                         <  2:8505741/150‑1 (MQ=255)

                                                                                                                             aaGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGc                         >  1:7141557/1‑150 (MQ=255)

                                                                                                                               gTGGGCGCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCg                       <  1:6318889/150‑1 (MQ=255)

                                                                                                                                tGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCgg                      >  2:836425/1‑150 (MQ=255)

                                                                                                                                tGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCgg                      >  2:3804258/1‑150 (MQ=255)

                                                                                                                                 gggCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGt                     <  1:5615087/150‑1 (MQ=255)

                                                                                                                                  ggCTCGGTGGAGGGCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTg                    <  2:7123279/150‑1 (MQ=255)

                                                                                                                                   gCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc                   <  2:6531092/150‑1 (MQ=255)

                                                                                                                                   gCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc                   <  1:115917/150‑1 (MQ=255)

                                                                                                                                   gCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc                   <  1:8001877/150‑1 (MQ=255)

                                                                                                                                     tCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCAc                 >  2:4869358/1‑150 (MQ=255)

                                                                                                                                           gAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCgg           <  2:5473416/150‑1 (MQ=255)

                                                                                                                                           caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCgg           >  2:6365937/2‑150 (MQ=255)

                                                                                                                                           caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCgg           >  2:5895205/2‑150 (MQ=255)

                                                                                                                                           caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCgg           >  1:1050/2‑150 (MQ=255)

                                                                                                                                            aGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc          >  1:3813935/1‑150 (MQ=255)

                                                                                                                                              gTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCg        >  2:5671319/1‑150 (MQ=255)

                                                                                                                                              gTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACAATCGGctc        >  1:3598110/1‑148 (MQ=255)

                                                                                                                                                 ggCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATg     >  2:2586815/1‑150 (MQ=255)

                                                                                                                                                  gCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGt    >  1:843213/1‑150 (MQ=255)

                                                                                                                                                  gCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGt    >  2:3540402/1‑150 (MQ=255)

                                                                                                                                                   cAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGtt   <  2:6853052/150‑1 (MQ=255)

GCACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTTC  >  NC_002947/3254506‑3254803

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