BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,000,732	T→C	6.3%	E131G (GAG→GGG) ‡	ycaC‑IV ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,000,732	0	T	C	6.3%	999.5 / 34.8	316	E131G (GAG→GGG) ‡	ycaC‑IV	putative hydrolase
Reads supporting (aligned to +/- strand): ref base T (160/135); new base C (9/11); total (170/146)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 342 total aligned reads displayed.

TGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAC  >  NC_002947/6000583‑6000881

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

  tCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAActct                                                                                                                                                     <  1:7637676/150‑1 (MQ=255)

   cGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACtctt                                                                                                                                                    <  2:5545171/150‑1 (MQ=255)

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

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

        cgcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTcc                                                                                                                                               <  1:8829856/150‑1 (MQ=255)

        cgcCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTcc                                                                                                                                               <  2:725212/150‑1 (MQ=255)

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

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

             gTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGc                                                                                                                                          >  2:4102445/1‑150 (MQ=255)

              tCGCGGTGCAGCTCTCAGGCCACGCCGAACCAGGTCATCTGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCTGGTCATGGCGTTGATGGTACCGGAGGCATCAGTGACCACGAACTCGTCGAACTCTTCTTCCAGGGcc                                                                                                                                         <  1:6773005/150‑1 (MQ=21)

                gcgGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGa                                                                                                                                       <  1:3063243/150‑1 (MQ=255)

                  ggTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACa                                                                                                                                     <  1:5626470/150‑1 (MQ=255)

                         cTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCgg                                                                                                                              <  1:2310147/150‑1 (MQ=255)

                         cTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCgg                                                                                                                              <  1:5304384/150‑1 (MQ=255)

                         cTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCgg                                                                                                                              <  1:3524109/150‑1 (MQ=255)

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

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

                           cacaGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGa                                                                                                                            <  1:4163241/150‑1 (MQ=255)

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

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

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

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

                               ggCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGc                                                                                                                        >  2:1647278/1‑150 (MQ=255)

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

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

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

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

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

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

                                          aCCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACtt                                                                                                             >  1:5092448/1‑150 (MQ=255)

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

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

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

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

                                                  atcaGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTcaca                                                                                                     >  1:3671946/1‑150 (MQ=255)

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

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

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

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

                                                             aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgat                                                                                          >  1:4734334/1‑150 (MQ=255)

                                                               ctgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatga                                                                                        >  1:4655068/1‑150 (MQ=255)

                                                                 gcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGa                                                                                                                                       <  2:5295776/101‑1 (MQ=255)

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

                                                                  cctgTCTCATGCGGTCGTGGGCGGCGTCGCGGGTCAGGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCTGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCAGAACGCCGGCgatgatga                                                                                     <  2:8945551/150‑3 (MQ=255)

                                                                  cctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCAACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCa                                                                                     <  2:3361600/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                              ggTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGc                                                                         >  2:2209353/1‑150 (MQ=255)

                                                                              ggTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGc                                                                         <  1:9073370/150‑1 (MQ=255)

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

                                                                                            tCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACt                                                           <  1:2504142/150‑1 (MQ=255)

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

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

                                                                                                    cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                   <  1:8406592/150‑1 (MQ=255)

                                                                                                    cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGACTTCACTGCCTTGAc                                                   <  1:145834/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                          gTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTc                                             >  2:9143266/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                 gTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTg                                      <  2:6328218/150‑1 (MQ=255)

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

                                                                                                                     cGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTccc                                  >  2:7095481/1‑150 (MQ=255)

                                                                                                                        aGGCATCAGTGACCACGAACACGTCGAACACTTCTTCCAGGGCCGACAGGGCCGGGAAAGCAACGCACACATCAGTCACAACGCCGGCGAGGATCAGTAGCTTCTAGCCAGTGGCCAACACTGCCTAGACGAAGTCTTCGATGTCCCAgg                               >  2:4946191/1‑150 (MQ=255)

                                                                                                                         ggCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGc                              <  2:407729/150‑1 (MQ=255)

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

                                                                                                                            aTCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                           <  2:4188119/150‑1 (MQ=255)

                                                                                                                            aTCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                           >  1:7303127/1‑150 (MQ=255)

                                                                                                                              cAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGa                         <  1:5783965/150‑1 (MQ=255)

                                                                                                                                 tGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCt                      <  2:2510846/150‑1 (MQ=255)

                                                                                                                                      aCGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCa                 >  1:1516623/1‑150 (MQ=255)

                                                                                                                                      aCGAACACGTCGAACTCTTCTTCCAGGGCCGAAAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCa                 >  2:1360099/1‑150 (MQ=255)

                                                                                                                                       cGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAg                >  1:5922846/1‑150 (MQ=255)

                                                                                                                                       cGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAg                >  1:5922946/1‑150 (MQ=255)

                                                                                                                                        gcacacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgg               <  2:6758949/148‑1 (MQ=255)

                                                                                                                                         aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg              <  1:8071427/150‑1 (MQ=255)

                                                                                                                                         aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg              <  1:1794437/150‑1 (MQ=255)

                                                                                                                                          acacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc             >  2:370380/1‑150 (MQ=255)

                                                                                                                                               tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa        >  2:569680/1‑150 (MQ=255)

                                                                                                                                                 gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt      >  1:4038177/1‑150 (MQ=255)

                                                                                                                                                 gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt      <  2:2489128/150‑1 (MQ=255)

                                                                                                                                                 gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt      >  1:6514165/1‑150 (MQ=255)

                                                                                                                                                  aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg     >  2:6238579/1‑150 (MQ=255)

                                                                                                                                                  aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg     >  1:3767940/1‑150 (MQ=255)

                                                                                                                                                   aCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGt    <  2:2069347/150‑1 (MQ=255)

                                                                                                                                                    ctctTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTa   <  1:9133913/150‑1 (MQ=255)

                                                                                                                                                     tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAg                                               >  2:5842641/1‑105 (MQ=255)

                                                                                                                                                     tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc  >  1:1198801/1‑150 (MQ=255)

TGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAC  >  NC_002947/6000583‑6000881

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