BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,882,257	Δ1 bp	5.9%	coding (1295/1356 nt)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,882,257	0	T	.	5.9%	319.7 / 12.7	101	coding (1295/1356 nt)	uacT	uric acid permease
Reads supporting (aligned to +/- strand): ref base T (48/47); new base . (2/4); total (50/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 101 total aligned reads displayed.

CAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCG  >  NC_002947/4882112‑4882401

cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGt                                                                                                                                                >  1:1964857/1‑148 (MQ=255)

cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGt                                                                                                                                                <  2:1964857/148‑1 (MQ=255)

cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGGGGATGGAGCCCATCACCCCAGGCGCCTTCTCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTc                                                                                                                                              >  1:2378574/1‑150 (MQ=255)

   cTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAACAGTAC‑GTTCAAc                                                                                                                                           <  1:2116682/150‑1 (MQ=255)

   cTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAACAGTAC‑GTTCAAc                                                                                                                                           <  2:1778291/150‑1 (MQ=255)

    tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACa                                                                                                                                          >  2:446078/1‑150 (MQ=255)

    tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACa                                                                                                                                          >  1:2053318/1‑150 (MQ=255)

    tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTAGACGAAGCCCATCCCCCACAGTTGAAT‑GCTGTGTTCCACGGTTTTTGCCCTGGTGTTGAACAGTAC‑GTTCAACa                                                                                                                                          <  1:414869/150‑1 (MQ=255)

         gttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCtt                                                                                                                                     <  2:530517/150‑1 (MQ=255)

            gCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGc                                                                                                                                  <  1:446078/150‑1 (MQ=255)

                     gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCgg                                                                                                                                >  1:1317361/1‑143 (MQ=255)

                     gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCgg                                                                                                                                <  2:1317361/143‑1 (MQ=255)

                     gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGc                                                                                                                             <  1:1347143/146‑1 (MQ=255)

                     gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGc                                                                                                                             >  2:1347143/1‑146 (MQ=255)

                      tCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGAcc                                                                                                                        >  1:706876/1‑150 (MQ=255)

                          cTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCgcgc                                                                                                                    >  2:2452281/1‑150 (MQ=255)

                           ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAACAGTAC‑GTTCAACa                                                                                                                                          >  2:20318/1‑127 (MQ=255)

                           ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAACAGTAC‑GTTCAACa                                                                                                                                          <  1:20318/127‑1 (MQ=255)

                           ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCa                                                                                                                   >  1:754463/1‑150 (MQ=255)

                            tCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCaa                                                                                                                  <  1:2232292/150‑1 (MQ=255)

                             cGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTAAACATCCTTGGCGGTGCCGACCGCGCAAt                                                                                                                 <  1:918145/150‑1 (MQ=255)

                              gggTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATg                                                                                                                >  1:266602/1‑150 (MQ=255)

                                gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGca                                                                                                              <  1:1050712/150‑1 (MQ=255)

                                   gATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  1:1650689/150‑1 (MQ=255)

                                   gATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  2:1745698/150‑1 (MQ=255)

                                   gATCCCCGTGGTGCGTCCGGAGTTCCTTGCACAGACGCCACGGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  2:123423/150‑1 (MQ=255)

                                    aTCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAac                                                                                                          <  2:2378859/150‑1 (MQ=255)

                                      ccccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                          <  1:496659/132‑1 (MQ=255)

                                      ccccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                          >  2:496659/1‑132 (MQ=255)

                                        ccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgacg                                                                                                      >  2:1478319/1‑150 (MQ=255)

                                        ccGTGGTACGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgacg                                                                                                      >  1:1311987/1‑150 (MQ=255)

                                         cGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGt                                                                                                     >  1:2298572/1‑150 (MQ=255)

                                           tgctgCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCt                                                                                                   <  1:1416191/147‑1 (MQ=255)

                                           gggtgCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCt                                                                                                   <  2:1475925/149‑1 (MQ=255)

                                                  tCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCa                                                                                            >  2:2383008/1‑150 (MQ=255)

                                                     ggAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCa                                                                                         >  1:335694/1‑150 (MQ=255)

                                                     ggAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGAATCGCCATG‑GCCACGGTCAAACCCCGATTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGGGCCGACCGCGCAATGCACAACGACGTGTGTCACCAGcg                                                                                         >  1:535098/1‑149 (MQ=255)

                                                      gAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCAt                                                                                        <  2:266602/150‑1 (MQ=255)

                                                       aGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCCGCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCAtt                                                                                       <  2:535098/150‑1 (MQ=255)

                                                          tCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAg                                                                                    <  1:1898570/150‑1 (MQ=255)

                                                            tttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                                  >  2:760296/1‑150 (MQ=255)

                                                            tttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                                  >  2:137794/1‑150 (MQ=255)

                                                             ttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCg                                                                                                                                 >  1:2156858/1‑102 (MQ=255)

                                                             ttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCg                                                                                                                                 <  2:2156858/102‑1 (MQ=255)

                                                                cccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGTGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgg                                                                              <  2:1167724/150‑1 (MQ=255)

                                                                cccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgg                                                                              <  2:1284825/150‑1 (MQ=255)

                                                                 ccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAggg                                                                             >  1:1587276/1‑150 (MQ=255)

                                                                 ccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAggg                                                                             >  2:252260/1‑150 (MQ=255)

                                                                  cAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggg                                                                            >  1:1703980/1‑150 (MQ=255)

                                                                    gaggCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggggg                                                                          <  2:285808/147‑1 (MQ=255)

                                                                    gATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggggg                                                                          <  1:1478319/150‑1 (MQ=255)

                                                                        ccACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGt                                                                      <  1:1558884/150‑1 (MQ=255)

                                                                         cacaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTTAACATCCTTGGCGGTGCCGACCGCGCCATGCACAACGACGTCTGTCACGAGCATTGAGCCAAGGGGGGCGGTg                                                                     >  1:1343318/1‑150 (MQ=255)

                                                                         cacaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTg                                                                     >  1:1491991/1‑150 (MQ=255)

                                                                          acaGTGGATGGAGCCAATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGc                                                                    >  2:394834/1‑150 (MQ=255)

                                                                            aGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCAt                                                                  <  1:137794/150‑1 (MQ=255)

                                                                               ggatggaGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGc                                                               <  2:754463/150‑1 (MQ=255)

                                                                                  tggaGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCAGTAGACGGTTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGc                                                            <  1:258090/150‑1 (MQ=255)

                                                                                    gaGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGccc                                                          <  2:2298572/150‑1 (MQ=255)

                                                                                       cccATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCg                                                                        >  1:1790927/1‑133 (MQ=255)

                                                                                       cccATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCg                                                                        <  2:1790927/133‑1 (MQ=255)

                                                                                       cccATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTc                                                       <  2:2031023/150‑1 (MQ=255)

                                                                                       cccATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTc                                                       <  2:706876/150‑1 (MQ=255)

                                                                                         cATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTACATGGCTGCCCCTCGt                                                     >  1:1304709/1‑150 (MQ=255)

                                                                                             aCCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGAGCATGGCTGCCCCTCGtgccc                                                 >  1:97532/1‑146 (MQ=255)

                                                                                                 acaGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCaa                                                                                                                                            <  2:1800768/55‑1 (MQ=255)

                                                                                                 acaGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCaa                                                                                                                                            >  1:1800768/1‑55 (MQ=255)

                                                                                                    gCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCGCCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          <  2:743090/150‑1 (MQ=255)

                                                                                                        tATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTaa                                      >  1:2217405/1‑150 (MQ=255)

                                                                                                         aTCGCCATG‑GCCATGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                     <  2:710441/150‑1 (MQ=255)

                                                                                                         aTCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                     <  2:1587276/150‑1 (MQ=255)

                                                                                                          tCGCCATG‑GCCACGGTCTGAGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                      >  2:101226/1‑132 (MQ=255)

                                                                                                          tCGCCATG‑GCCACGGTCTGAGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                      <  1:101226/132‑1 (MQ=255)

                                                                                                            gCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                >  1:233589/1‑120 (MQ=255)

                                                                                                            gCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                <  2:233589/120‑1 (MQ=255)

                                                                                                              cATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     >  1:2371670/1‑97 (MQ=255)

                                                                                                              cATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     <  2:2371670/97‑1 (MQ=255)

                                                                                                                     cACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGt                          >  2:880052/1‑150 (MQ=255)

                                                                                                                           cTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     <  2:863807/85‑1 (MQ=255)

                                                                                                                           cTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     >  1:863807/1‑85 (MQ=255)

                                                                                                                            tgtgCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccg                   >  1:1199182/1‑150 (MQ=255)

                                                                                                                              tgCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgcc                 >  1:874866/1‑150 (MQ=255)

                                                                                                                                cccTGGTGTTGAAC‑GTACTGTGCAACAGCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCgg               >  2:1403240/1‑150 (MQ=255)

                                                                                                                                caaggatgtTGAACAGTAC‑GTTCAACACCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGCGCGGTGCATGGCTGCCCCTCGCTCACCTCTGCAGTAACCGCTACCCCGTCCACCCGCCgg               >  2:518758/7‑150 (MQ=255)

                                                                                                                                 ccTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGa              >  2:2021026/1‑150 (MQ=255)

                                                                                                                                 ccTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGa              >  1:1939550/1‑150 (MQ=255)

                                                                                                                                  cTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGaa             <  2:2053318/150‑1 (MQ=255)

                                                                                                                                   tGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          >  2:539255/1‑120 (MQ=255)

                                                                                                                                   tGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          <  1:539255/120‑1 (MQ=255)

                                                                                                                                   tGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAAt            <  2:874866/150‑1 (MQ=255)

                                                                                                                                   tGGNGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAAt            >  2:1516874/1‑150 (MQ=255)

                                                                                                                                     gtgtTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                <  1:2450863/96‑1 (MQ=255)

                                                                                                                                     gtgtTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                >  2:2450863/1‑96 (MQ=255)

                                                                                                                                       gtTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGACTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGt        <  2:1311987/150‑1 (MQ=255)

                                                                                                                                          gAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGc     >  2:2246585/1‑150 (MQ=255)

                                                                                                                                             c‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCg  >  2:2311491/1‑150 (MQ=255)

                                                                                                                                                tACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAcc                                    >  1:1875107/1‑114 (MQ=255)

                                                                                                                                                tACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAcc                                    <  2:1875107/114‑1 (MQ=255)

                                                                                                                                                  cTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGcac                                                                                                             <  2:2132671/39‑1 (MQ=255)

                                                                                                                                                  cTGTTCAACATCCTTGGCGGTGCAGACCGCGCAATGcac                                                                                                             >  1:2132671/1‑39 (MQ=255)

CAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATG‑GCCACGGTCTGTGCCCTGGTGTTGAAC‑GTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCG  >  NC_002947/4882112‑4882401

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