BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,000,710	A→T	6.3%	T138T (ACT→ACA)	ycaC‑IV ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,000,710	0	A	T	6.3%	1056.3 / 37.1	319	T138T (ACT→ACA)	ycaC‑IV	putative hydrolase
Reads supporting (aligned to +/- strand): ref base A (165/132); new base T (9/12); total (175/144)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

CACAGCGCGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAT  >  NC_002947/6000561‑6000859

cacaGCGCGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCa                                                                                                                                                       <  1:5998976/150‑1 (MQ=255)

cacaGCGCGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCa                                                                                                                                                       <  1:4275854/150‑1 (MQ=255)

  caGCGCGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGt                                                                                                                                                     >  1:6758949/1‑150 (MQ=255)

    gcgcGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGa                                                                                                                                                   <  1:2296295/150‑1 (MQ=255)

           cAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGaa                                                                                                                                            >  1:5161732/1‑150 (MQ=255)

           cAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGaa                                                                                                                                            >  1:3064948/1‑150 (MQ=255)

                 ttCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTc                                                                                                                                      >  1:7022847/1‑150 (MQ=255)

                 ttCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTc                                                                                                                                      >  1:4188119/1‑150 (MQ=255)

                 ttCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTc                                                                                                                                      >  1:3715692/1‑150 (MQ=255)

                 ttCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTc                                                                                                                                      >  2:6992849/1‑150 (MQ=255)

                 ttCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTc                                                                                                                                      <  1:193283/150‑1 (MQ=255)

                    gATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGaa                                                                                                                                   >  2:1570470/1‑150 (MQ=255)

                      tGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAAct                                                                                                                                 <  1:6620655/150‑1 (MQ=255)

                          gTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCa                                                                                                                                                       >  2:3058786/1‑124 (MQ=255)

                          gTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACtcttc                                                                                                                             >  1:2142048/1‑150 (MQ=255)

                          gTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACtcttc                                                                                                                             >  1:6323821/1‑150 (MQ=255)

                               gcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCGCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCa                                                                                                                        >  1:9107055/1‑150 (MQ=255)

                               gcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCa                                                                                                                        <  2:415178/150‑1 (MQ=255)

                               gcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCa                                                                                                                        >  1:6443136/1‑150 (MQ=255)

                               gcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCa                                                                                                                        >  1:2510846/1‑150 (MQ=255)

                                  aGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAggg                                                                                                                     >  2:1407267/1‑150 (MQ=255)

                                            cAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGc                                                                                                           <  1:3158977/150‑1 (MQ=255)

                                             aGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGcc                                                                                                          >  1:1599524/1‑150 (MQ=255)

                                              gCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCg                                                                                                         <  1:4046098/150‑1 (MQ=255)

                                              gCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCg                                                                                                         <  2:6216793/150‑1 (MQ=255)

                                               cTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCgg                                                                                                        >  1:8841550/1‑150 (MQ=255)

                                                tCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCggg                                                                                                       >  1:8466812/1‑150 (MQ=255)

                                                tCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCggg                                                                                                       >  1:7303552/1‑150 (MQ=255)

                                                  acaGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGaa                                                                                                     >  1:4352272/1‑150 (MQ=255)

                                                   caGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGaaa                                                                                                    >  1:1511681/1‑150 (MQ=255)

                                                   caGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGaaa                                                                                                    <  1:8365500/150‑1 (MQ=255)

                                                    aGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAg                                                                                                   >  2:5755970/1‑150 (MQ=255)

                                                    aGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAg                                                                                                   >  1:2739549/1‑150 (MQ=255)

                                                     ggCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGc                                                                                                  <  2:3356690/150‑1 (MQ=255)

                                                       ccACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTa                                                                                                >  2:4081094/1‑150 (MQ=255)

                                                        cACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCa                                                                                                                                                       <  2:2682976/94‑1 (MQ=255)

                                                        cACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTAc                                                                                               >  1:3660960/1‑150 (MQ=255)

                                                         aCGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACg                                                                                              >  2:7491201/1‑150 (MQ=255)

                                                           gCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACGCGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGCAAGCTACGca                                                                                            >  1:8001025/1‑150 (MQ=255)

                                                           gCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGca                                                                                            >  1:1121117/1‑150 (MQ=255)

                                                                  cAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCa                                                                                     <  2:5024448/150‑1 (MQ=255)

                                                                  cAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCa                                                                                     >  2:1919834/1‑150 (MQ=255)

                                                                    ggTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGt                                                                                   <  2:5560617/150‑1 (MQ=255)

                                                                     gTCATCAGTTGCTCACCTGCCTGGCTCTTGCGGTCGTGGGCTGCGTCGCGGGTCATTGCGTTGAAGTTACCTGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTc                                                                                  <  2:8552802/150‑1 (MQ=255)

                                                                     gTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTc                                                                                  <  1:3764298/150‑1 (MQ=255)

                                                                     gTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTc                                                                                  <  2:6599228/150‑1 (MQ=255)

                                                                      tcatcaGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTca                                                                                 >  1:2091948/1‑150 (MQ=255)

                                                                      tcatcaGTTGCGCACCTGCCTGGCTCAGGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGGGGCAGGGACCGGAGGCCGCCGTGACCACGAACACGGCGAACTCTTCTTCCCGAGCCGACAGGGCCGGGCCAGCTCCGCCCACTTCAGGca                                                                                 >  2:939894/1‑150 (MQ=255)

                                                                       catcaGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTcac                                                                                <  1:7567577/150‑1 (MQ=255)

                                                                          caGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAAc                                                                             <  1:3342544/150‑1 (MQ=255)

                                                                          caGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAAc                                                                             >  1:3072687/1‑150 (MQ=255)

                                                                              tGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCg                                                                         <  1:4531790/150‑1 (MQ=255)

                                                                                cgcACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGc                                                                       <  2:8241430/150‑1 (MQ=255)

                                                                                 gcACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCg                                                                      >  1:4737530/1‑150 (MQ=255)

                                                                                  cACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCga                                                                     <  2:4362417/150‑1 (MQ=255)

                                                                                   aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgag                                                                    >  2:5664138/1‑149 (MQ=255)

                                                                                     ctgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatga                                                                  <  2:7212839/150‑1 (MQ=255)

                                                                                      tgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatgat                                                                 <  1:4689583/150‑1 (MQ=255)

                                                                                       gcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATc                                                                >  2:1037219/1‑150 (MQ=255)

                                                                                         ctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAg                                                              <  1:5891228/150‑1 (MQ=255)

                                                                                              tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                         >  2:6837281/1‑150 (MQ=255)

                                                                                              tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                         >  2:4586243/1‑150 (MQ=255)

                                                                                              tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                         >  1:6016204/1‑150 (MQ=255)

                                                                                              tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                         >  1:8784514/1‑150 (MQ=255)

                                                                                                   cGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTg                                                    <  2:697428/150‑1 (MQ=255)

                                                                                                   cGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTg                                                    >  2:4485728/1‑150 (MQ=255)

                                                                                                         tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTg                                              <  2:912569/150‑1 (MQ=255)

                                                                                                         tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTg                                              <  2:2152548/150‑1 (MQ=255)

                                                                                                         tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGAACAGTTGCTTCTTGCCAGTg                                              <  2:3180350/150‑1 (MQ=255)

                                                                                                                     cgGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGcc                                  <  1:8912180/150‑1 (MQ=255)

                                                                                                                       ggTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCtt                                <  2:847225/150‑1 (MQ=255)

                                                                                                                       ggTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCtt                                <  2:396176/150‑1 (MQ=255)

                                                                                                                         tCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGa                              <  2:8898191/150‑1 (MQ=255)

                                                                                                                          cacttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGAc                             <  2:1786679/143‑1 (MQ=255)

                                                                                                                            tGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGa                           >  1:5004311/1‑150 (MQ=255)

                                                                                                                             ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                          >  1:1030459/5‑150 (MQ=255)

                                                                                                                             ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                          <  2:1511681/150‑1 (MQ=255)

                                                                                                                             ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                          <  1:2050370/150‑1 (MQ=255)

                                                                                                                                gTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTc                       <  2:3477271/150‑1 (MQ=255)

                                                                                                                                 ttGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTcc                      >  1:2323280/1‑149 (MQ=255)

                                                                                                                                  tGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCtt                     <  2:5458944/150‑1 (MQ=255)

                                                                                                                                   gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                    >  2:4220440/1‑150 (MQ=255)

                                                                                                                                   gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                    <  2:3757191/150‑1 (MQ=255)

                                                                                                                                    aaGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCg                   >  1:2972637/1‑150 (MQ=255)

                                                                                                                                     aGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGt                  >  1:6256065/1‑150 (MQ=255)

                                                                                                                                     aGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGt                  >  2:3572904/1‑150 (MQ=255)

                                                                                                                                      ggTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGATTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCAtt                 <  1:7648907/150‑1 (MQ=255)

                                                                                                                                      ggTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGtt                 <  2:24068/150‑1 (MQ=255)

                                                                                                                                        tACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGt               <  2:2796374/150‑1 (MQ=255)

                                                                                                                                         aCCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTc              <  2:2879724/150‑1 (MQ=255)

                                                                                                                                          ccGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTcc             >  2:8955348/1‑150 (MQ=255)

                                                                                                                                           ccgacgCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCAACGCAGACTTCGGTCACAACACATGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTcac            >  2:4744878/6‑148 (MQ=255)

                                                                                                                                             gAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAg          >  2:7980534/1‑150 (MQ=255)

                                                                                                                                              aGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAgg         >  2:4946160/1‑150 (MQ=255)

                                                                                                                                                 cATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGt      >  1:5745484/1‑150 (MQ=255)

                                                                                                                                                  aTCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGtt     <  1:5530245/150‑1 (MQ=255)

                                                                                                                                                  aTCAGTGACCACGAACACGTCGAACTCTTGTTCGAGGGCCGAGAGGGACTTTAAAGCTACGCTCCCTTGCGTCACAACTCCGGCGATGATGAGGTGGGTCTTGCGACGTGACTTCTCTTCCTTGACGTAGTCTTCGTTGTCCCAGGCGtt     <  2:7644962/150‑1 (MQ=255)

                                                                                                                                                  aTCAGGGACCAAGAACACGTCGAACTCTTCTTCCAGGGCCGACCGGGCCGGGACAGCTAAGCACACATCAGTCACAACGCCGGCGATGATCAGTGGCTTCTTGCCAGTGGCCTGCACTGCCTGGACGACGTCTTCGTTGTCCCAGGCGtt     >  2:6366672/1‑150 (MQ=255)

                                                                                                                                                   tCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTg    <  2:813177/150‑1 (MQ=255)

                                                                                                                                                     aGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt  <  1:1429672/150‑1 (MQ=255)

CACAGCGCGGCCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAT  >  NC_002947/6000561‑6000859

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