BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,249,705	A→C	8.1%	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,249,705	0	A	C	8.1%	363.4 / 20.9	133	K428Q (AAG→CAG)	ureC	urease subunit alpha
Reads supporting (aligned to +/- strand): ref base A (58/63); new base C (6/5); total (65/68)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 135 total aligned reads displayed.

CAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTT  >  NC_002947/3249561‑3249853

cAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTg                                                                                                                                                 >  1:1804345/1‑150 (MQ=255)

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

  gATGAAGCTGCGTCGCGGTCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                               >  2:1736687/1‑150 (MQ=255)

   aTGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGcc                                                                                                                                              >  2:2693642/1‑150 (MQ=255)

    tGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCg                                                                                                                                             >  2:2198498/1‑150 (MQ=255)

     gAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGa                                                                                                                                            <  2:1595472/150‑1 (MQ=255)

       aGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCCTGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAcc                                                                                                                                          <  1:1428085/150‑1 (MQ=255)

        gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCCCCCATTGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGACCt                                                                                                                                         <  1:1467444/150‑1 (MQ=255)

        gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCt                                                                                                                                         >  1:1731602/1‑150 (MQ=255)

         cTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTg                                                                                                                                        <  2:600623/150‑1 (MQ=255)

              tCGCGGGCGGCTGGCACCGCACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCt                                                                                                                                   <  1:504426/150‑1 (MQ=255)

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

                    gccacgggCACCGGACACCCCAGACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                             <  1:432678/144‑1 (MQ=255)

                    gCCGCTGGCACCGGACACCCCGTACAGCCGCTCCTTCCAGGTCAATCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                             <  2:2803800/150‑1 (MQ=255)

                             aCCGGACAACCCGTACAGCGACAACTTCCGCGTCAACCGCTACATCGCCAAGTACACCATCAACCCGCAACTCACCCATGGCATCGGTCACGAAGTGGGATCGATGGAGGTCGGCTACCTGGACCACCTGCTGCTATGGTCGACAGCgcc                                                                                                                    >  2:267356/1‑148 (MQ=255)

                              ccGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTc                                                                                                                   <  1:1250967/150‑1 (MQ=255)

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

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

                                   cacCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGc                                                                                                              <  1:1137635/150‑1 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                >  2:2565235/1‑129 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                <  1:2621417/129‑1 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                <  1:2565235/129‑1 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                >  2:2621417/1‑129 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGTCa                                                                                                                         >  2:1193126/1‑136 (MQ=255)

                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGTCa                                                                                                                         <  1:1193126/136‑1 (MQ=255)

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

                                       ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTgg                                                                                                                                                >  2:243827/1‑112 (MQ=255)

                                       ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTgg                                                                                                                                                <  1:243827/112‑1 (MQ=255)

                                         gTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCaa                                                                                                        >  2:1514694/1‑150 (MQ=255)

                                           acacCGACAACTGCCGGCTCAAGCGCTACATCGCCAAGTACCCCCTCAACCCGGCACGCACCCATGGCATCCCTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGc                                                                                                      <  2:1496360/150‑1 (MQ=255)

                                                      ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCaa                                                                                                                                                     >  2:1371864/1‑92 (MQ=255)

                                                      ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGtt                                                                                                                    >  2:541409/1‑125 (MQ=255)

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

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

                                                      ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                           <  1:1892434/150‑1 (MQ=255)

                                                       tCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGc                                                                                          >  2:2706007/1‑150 (MQ=255)

                                                        ccGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGACCtcc                                                                                                                                       <  2:2567683/104‑3 (MQ=255)

                                                         cGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTc                                                                                        >  1:3043233/1‑150 (MQ=255)

                                                            gTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTg                                                                                                                                        <  1:601628/99‑1 (MQ=255)

                                                              cAAGCGCTACATCGCCAAGTACAGCATCACCCCGGCACTCGCCCGTGGCAGGGGTCACGAAGTGGGCTCGGTGGAGGTCGGAAAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCCTTGCGGTCAAGCCGGCCCTGGTGCTCAAggg                                                                                   <  2:64165/150‑1 (MQ=255)

                                                               aaGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGc                                                                                  >  1:2022824/1‑150 (MQ=255)

                                                                cgcgcTACCTCGCCACGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCg                                                                                 <  2:174845/149‑1 (MQ=255)

                                                                aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCttt                                                                                                                <  1:2321813/119‑1 (MQ=255)

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

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

                                                                 gcgcTACATCGCCAAGGACACCATCAACCCGGCCCTCACCCATGGCCTCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTTGAGCACCAGGGccg                                                                                <  2:1731602/150‑3 (MQ=255)

                                                                   gcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGga                                                                              <  2:2128284/150‑1 (MQ=255)

                                                                    cTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgat                                                                             <  1:1179006/150‑1 (MQ=255)

                                                                     tACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTg                                                                                                                                                 <  1:2273206/81‑1 (MQ=255)

                                                                     tACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTg                                                                                                                                                 >  2:2273206/1‑81 (MQ=255)

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

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

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

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

                                                                           gCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTc                                                                      >  1:1034222/1‑150 (MQ=255)

                                                                            ccAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCa                                                                     >  2:3040751/1‑150 (MQ=255)

                                                                             cAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCAc                                                                    <  2:195218/150‑1 (MQ=255)

                                                                                gTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCg                                                                                                                      <  1:2210426/97‑1 (MQ=255)

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

                                                                                      cATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTg                                                                                                               >  2:910501/1‑98 (MQ=255)

                                                                                      cATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTg                                                                                                               <  1:910501/98‑1 (MQ=255)

                                                                                      cATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCAt                                                           >  2:1382289/1‑150 (MQ=255)

                                                                                       aTCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGAc                                                                                                                                           <  1:99853/69‑1 (MQ=39)

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

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

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

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

                                                                                         cAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATggg                                                        <  1:2589176/150‑1 (MQ=255)

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

                                                                                           aCCCGGCCCTTACCCATGGCATCGGTCACGAAGAGGGCTCGGTCGGGTGCGGCCCGCTGGCGGACCGGGTGCTGTGGTCGCCAGCGTTCGTTCCGGTCAAGGCGGCCCTGGGACTCAGTGTCGGGCTGATCGTCACCGCGCCCATGGGCg                                                      >  2:1119289/1‑150 (MQ=255)

                                                                                            cccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCGGCGTTCttt                                                                                                                <  2:2692421/91‑1 (MQ=255)

                                                                                            cccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCGGCGTTCttt                                                                                                                >  1:2692421/1‑91 (MQ=255)

                                                                                            cccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                         <  1:2883058/130‑1 (MQ=255)

                                                                                             ccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGAc                                                    >  2:1834429/1‑150 (MQ=255)

                                                                                              cGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTACTTTGCGGTCAAGCCGACCCTGGTGCTGAAGGTCGTGATGATCGGCACCGCGTCCATGGTCGACa                                                   <  2:2806954/150‑1 (MQ=255)

                                                                                                 cACTCACCCATGGCATCGGTCACGAAGTGGGCTCGCAGGAGGTCCGCACGTTGACCGACCTGGTGCTGTGGGCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCa                                                >  1:59349/1‑150 (MQ=255)

                                                                                                  aCTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCaa                                               <  1:3016852/150‑1 (MQ=255)

                                                                                                     cACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACgg                                            <  2:2738987/150‑1 (MQ=255)

                                                                                                      aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGc                                           >  2:1296767/1‑150 (MQ=255)

                                                                                                           tGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCAt                                      <  2:2022824/150‑1 (MQ=255)

                                                                                                             gCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACGTCAACGGCTCCATac                                    <  2:3043233/150‑1 (MQ=255)

                                                                                                                     cACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCAccccc                            >  1:1851441/1‑150 (MQ=255)

                                                                                                                      aCGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATaccca                                 <  1:228845/144‑1 (MQ=255)

                                                                                                                      aCGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATaccca                                 >  2:228845/1‑144 (MQ=255)

                                                                                                                      aCGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCAccccc                            <  2:1451981/149‑1 (MQ=255)

                                                                                                                      aCGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACAGCGCCCATGGGCGACATCAACGGCTCCATACCCAccccc                            >  1:1451981/1‑149 (MQ=255)

                                                                                                                            tGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCgg                     <  1:1306539/150‑1 (MQ=255)

                                                                                                                                 tCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCAc                <  2:1755587/150‑1 (MQ=255)

                                                                                                                                     tGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATc            >  2:76934/1‑150 (MQ=255)

                                                                                                                                     accaGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                         >  1:1845620/4‑89 (MQ=39)

                                                                                                                                      ccaGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCg           >  1:1645921/3‑150 (MQ=255)

                                                                                                                                      ccaGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCg           >  1:889303/3‑150 (MQ=255)

                                                                                                                                       caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCgg          <  2:889303/149‑1 (MQ=255)

                                                                                                                                        aGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc         <  2:59349/150‑1 (MQ=255)

                                                                                                                                        aGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc         <  1:448771/150‑1 (MQ=255)

                                                                                                                                        aGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGCGACGATCTTTCTTTGCGGTCAAGCCGGCCCTGTTGCTCAAGGGCGGGATGATCGTCACCGCGCCCGTGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc         <  2:1872429/150‑1 (MQ=255)

                                                                                                                                          gTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCg       >  1:2453622/1‑150 (MQ=255)

                                                                                                                                             ggCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCa                                       >  2:2153345/1‑115 (MQ=255)

                                                                                                                                             ggCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCa                                       <  1:2153345/115‑1 (MQ=255)

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

CAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTT  >  NC_002947/3249561‑3249853

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