BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,446,670	A→C	5.6%	L43R (CTG→CGG)	nicF ←	maleamate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,446,670	0	A	C	5.6%	552.8 / 13.3	173	L43R (CTG→CGG)	nicF	maleamate amidohydrolase
Reads supporting (aligned to +/- strand): ref base A (73/89); new base C (5/6); total (78/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

CCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGATTTCAGTAG  >  NC_002947/4446521‑4446819

ccGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACa                                                                                                                                                       <  1:859443/150‑1 (MQ=255)

 cGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAg                                                                                                                                                      >  2:5405424/1‑150 (MQ=255)

     gcCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGt                                                                                                                                                  <  2:549981/150‑1 (MQ=255)

     gcCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGt                                                                                                                                                  >  2:2382493/1‑150 (MQ=255)

       cTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTgc                                                                                                                                                <  1:5414282/150‑1 (MQ=255)

       cTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTgc                                                                                                                                                <  1:4258222/150‑1 (MQ=255)

         tGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTgcgc                                                                                                                                              >  1:3747832/1‑150 (MQ=255)

             cacccaccCGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgg                                                                                                                                          <  2:2428933/150‑1 (MQ=255)

             cACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgg                                                                                                                                          <  1:4873606/150‑1 (MQ=255)

             cACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgg                                                                                                                                          <  1:69901/150‑1 (MQ=255)

              aCCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCggg                                                                                                                                         <  2:4988458/150‑1 (MQ=255)

               cccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGGATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCC‑GGGCGCACCGCTGTACGCCCcgggc                                                                                                                                       <  2:4307718/150‑2 (MQ=255)

                ccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGTATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggt                                                                                                                                       <  1:4197393/150‑1 (MQ=255)

                ccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTCTGCACCACCAAGGTCCCGCCGTCGCGGGCCCGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACCGCGGTGCGCCCggggg                                                                                                                                       >  2:3904590/1‑149 (MQ=255)

                ccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggt                                                                                                                                       <  2:3039818/150‑1 (MQ=255)

                ccACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCCGGGCGCACCGCTGTACGCCCcgggc                                                                                                                                       <  1:4840034/150‑2 (MQ=255)

                 cacaCGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCAC‑CCC‑GGGCGCACCGCTGTACGCCCcgggcgtg                                                                                                                                    <  2:381300/150‑5 (MQ=38)

                 cacaCGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAACAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtg                                                                                                                                      <  1:4237450/150‑1 (MQ=255)

                  acacGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtgg                                                                                                                                     >  1:842491/1‑150 (MQ=255)

                  acacGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtgg                                                                                                                                     >  2:2012432/1‑150 (MQ=255)

                     cgccgcCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGgtgt                                                                                                                                  <  1:583458/150‑1 (MQ=255)

                      gccgccGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTa                                                                                                                                 >  2:2078651/1‑150 (MQ=255)

                      gccgccGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGCACAGCGGTGCGCCCGGGGTGGTGTa                                                                                                                                 <  1:4444651/150‑1 (MQ=255)

                         gccgTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGcc                                                                                                                              >  1:3443304/1‑150 (MQ=255)

                               aGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCAt                                                                                                                        <  1:5226669/150‑1 (MQ=255)

                                   aaaTGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAg                                                                                                                    >  1:1957619/1‑150 (MQ=255)

                                    aaTGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGt                                                                                                                   <  1:5404934/150‑1 (MQ=255)

                                      tGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCg                                                                                                                 <  1:1449255/150‑1 (MQ=255)

                                        cGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGat                                                                                                               <  1:300557/150‑1 (MQ=255)

                                           aGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcat                                                                                                            >  1:742783/1‑150 (MQ=255)

                                            ggCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcatc                                                                                                           <  1:2955175/150‑1 (MQ=255)

                                                    gCGGATGTTGCTATGCACCACCAAGGTACCGCAGACGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGCTGCGCCCGGGGTGGTGTAGCCCTGCTTGAAGTCGATCAACAGCAGGGc                                                                                                   >  2:2841132/1‑150 (MQ=255)

                                                       gATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCgg                                                                                                <  1:3914336/150‑1 (MQ=255)

                                                        aTGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGc                                                                                               <  1:5326220/150‑1 (MQ=255)

                                                            tGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGc                                                                                           >  1:1648633/1‑150 (MQ=255)

                                                             ggTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCg                                                                                          <  2:2276731/150‑1 (MQ=255)

                                                              ggaTGCCCCCCCAAGGGACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGCTCATCAGCAGGGCCGGCTTGCgg                                                                                         <  2:4765536/148‑1 (MQ=255)

                                                              gTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCgg                                                                                         <  1:1008057/150‑1 (MQ=255)

                                                                   caccacCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCaaa                                                                                    <  2:5333713/150‑1 (MQ=255)

                                                                   caccacCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCaaa                                                                                    <  2:2575451/150‑1 (MQ=255)

                                                                    accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAg                                                                                   >  2:1924939/1‑150 (MQ=255)

                                                                    accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAg                                                                                   >  2:3485740/1‑150 (MQ=255)

                                                                     ccaccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGc                                                                                  <  1:2238507/150‑1 (MQ=255)

                                                                       accaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtggt                                                                                                                                    <  2:3100888/98‑1 (MQ=255)

                                                                       accaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGggtggt                                                                                                                                    >  1:3100888/1‑98 (MQ=255)

                                                                                ccGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACCGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTgg                                                                       >  2:4222938/1‑150 (MQ=255)

                                                                                  gCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCGCGACTGCCGCGCCCACGCCCGGGGCGTACCGCGGTGCGCCCGGGGTGGGGTAGCCCTGCATGAAGTCGATCATCAGCCGGGCCGGCGTGCGGCCAAAGCCGATGCGCTGGcc                                                                     >  1:1830539/1‑150 (MQ=255)

                                                                                     gTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCca                                                                  <  2:635899/150‑1 (MQ=255)

                                                                                      tCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCcac                                                                 >  1:1780636/1‑150 (MQ=255)

                                                                                          gggCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCgg                                                                                         <  2:4793725/122‑1 (MQ=255)

                                                                                          gggCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCgg                                                                                         >  1:4793725/1‑122 (MQ=255)

                                                                                          gggCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCacacac                                                             <  1:866161/150‑1 (MQ=255)

                                                                                           ggCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAcc                                                            <  1:14403/150‑1 (MQ=255)

                                                                                            gccagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAccc                                                           <  1:323141/150‑1 (MQ=255)

                                                                                             ccagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCt                                                          >  2:2183864/1‑150 (MQ=255)

                                                                                             acagcgccagcAGGCCGGCGGCCTGCTCGACTGACGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGTTGGTGTAGACCTGCATGAAGTCTATCATCAGCAGGGCCGGCTAGCGGCAAAAGCCGATGCGCTGGCCCCACACACCCt                                                          >  2:3018339/2‑150 (MQ=255)

                                                                                              cagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTg                                                         >  2:41287/1‑150 (MQ=255)

                                                                                              cagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTg                                                         >  1:353590/1‑150 (MQ=255)

                                                                                              cagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTg                                                         >  2:4589105/1‑150 (MQ=255)

                                                                                                  tcgagcagGCCGCCGGCCTGCTGGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAg                                                     >  1:5289749/4‑150 (MQ=255)

                                                                                                    cagcagGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGtt                                                   >  2:2528009/1‑150 (MQ=255)

                                                                                                    cagcagGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGtt                                                   <  2:842491/150‑1 (MQ=255)

                                                                                                     agcagGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGt                                                                                                                   >  2:3213184/1‑85 (MQ=255)

                                                                                                      gcagGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGt                                                 >  1:3376179/1‑150 (MQ=255)

                                                                                                          gCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACg                                             >  1:1379603/1‑150 (MQ=255)

                                                                                                               cggcCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGgcgc                                        <  2:5014112/150‑1 (MQ=255)

                                                                                                                  ccTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCACTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCttt                                     >  1:158990/1‑150 (MQ=255)

                                                                                                                   cTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtg                                    >  1:1880656/1‑150 (MQ=255)

                                                                                                                    tGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAccc                                                           >  2:2277600/1‑126 (MQ=255)

                                                                                                                    tGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAccc                                                           <  1:2277600/126‑1 (MQ=255)

                                                                                                                    tGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgt                                   >  2:5015591/1‑150 (MQ=255)

                                                                                                                    tGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgt                                   <  1:3585155/150‑1 (MQ=255)

                                                                                                                    tGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgt                                   <  2:1725720/150‑1 (MQ=255)

                                                                                                                      cTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGc                                 >  2:1206336/1‑150 (MQ=255)

                                                                                                                        cGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGt                               <  1:2727816/150‑1 (MQ=255)

                                                                                                                         gACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTc                              >  1:5140260/1‑150 (MQ=255)

                                                                                                                          aCTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatc                                                                                                              <  1:3721651/69‑1 (MQ=255)

                                                                                                                          aCTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatc                                                                                                              >  2:3721651/1‑69 (MQ=255)

                                                                                                                           cTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGc                            <  1:1800459/150‑1 (MQ=255)

                                                                                                                            tGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGGTGCGCTGGCCCCACACACCCTGGGAGTTGTCACGGGCGCTTTGTGCGTCGCt                           >  2:4932235/1‑150 (MQ=255)

                                                                                                                            tGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCt                           >  1:5027608/1‑150 (MQ=255)

                                                                                                                              ccGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGGGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTGTGGGCGTCGCTCa                         >  2:4455938/1‑150 (MQ=255)

                                                                                                                               cgcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCAt                        >  2:2003873/1‑150 (MQ=255)

                                                                                                                               cgcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCAt                        >  1:593557/1‑150 (MQ=255)

                                                                                                                                gcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATg                       >  2:4874031/1‑150 (MQ=255)

                                                                                                                                 cgACCACGCCCGGGGTGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGt                      >  1:26695/1‑150 (MQ=255)

                                                                                                                                 cgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGt                      >  2:3796096/1‑150 (MQ=255)

                                                                                                                                 cgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGt                      >  2:2405119/1‑150 (MQ=255)

                                                                                                                                     cACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGt                  >  1:1054453/1‑150 (MQ=255)

                                                                                                                                      aCGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGtc                 >  2:5284599/1‑150 (MQ=255)

                                                                                                                                       cacccGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGtct                >  1:4017252/3‑150 (MQ=255)

                                                                                                                                        gCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGtctc               >  1:3467654/1‑150 (MQ=255)

                                                                                                                                        gCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGtctc               <  2:683940/150‑1 (MQ=255)

                                                                                                                                        acccGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGtctc               <  2:4017252/149‑1 (MQ=255)

                                                                                                                                            ggggCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGa           <  1:1206336/150‑1 (MQ=255)

                                                                                                                                             gggCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGAt          <  1:2573024/150‑1 (MQ=255)

                                                                                                                                             gggCGTACAGCGGTGCGCCCGGGGAGGTGGAGCCCTGCATGACGTCGCTCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCGGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGAt          <  2:2989593/150‑1 (MQ=255)

                                                                                                                                               gCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGAttt        <  2:5140260/150‑1 (MQ=255)

                                                                                                                                                cGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTc                                                <  1:2137598/109‑1 (MQ=255)

                                                                                                                                                     aGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGATTTCAGTAg  <  1:4589105/150‑1 (MQ=255)

CCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGTCTCCCGATTTCAGTAG  >  NC_002947/4446521‑4446819

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