BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,438,254	A→T	5.5%	K457M (AAG→ATG) ‡	PP_3932 →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,438,254	0	A	T	5.5%	387.6 / 11.0	127	K457M (AAG→ATG) ‡	PP_3932	diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base A (57/63); new base T (3/4); total (60/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

Only 100 of 130 total aligned reads displayed.

CCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTAC  >  NC_002947/4438105‑4438403

ccgcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCaa                                                                                                                                                        <  1:1643738/150‑1 (MQ=255)

 cgcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                                                       <  1:1368494/150‑1 (MQ=255)

 cgcgccATACAGGACCGGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                                                       <  1:229100/150‑1 (MQ=255)

     ccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGccc                                                                                                                                                   >  1:339082/1‑150 (MQ=255)

      cATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCg                                                                                                                                                  >  2:255060/1‑150 (MQ=255)

         aCAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGTCCAAGGCCCgggg                                                                                                                                               <  1:1014152/150‑1 (MQ=255)

         aCAGGACCTGCACGTGACACAGGCCGACCAGCGAATGCCGATTAAAACGTGCAGAAGTTTGGCAGTCTGCAAGTGAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGGGCCGACCAGGCGCTGTACCGGGCCAAGGCCCgggg                                                                                                                                               >  1:1863182/1‑150 (MQ=255)

             gACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGc                                                                                                                                           >  2:360585/1‑150 (MQ=255)

               ccTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCaa                                                                                                                                         <  1:529937/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGTGCCAAGGCCCGGGGCCGCAACCGCTTc                                                                                                                                  <  2:1172044/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCg                                                                                                                                                  >  2:417216/1‑134 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                  <  2:1035156/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                  <  1:360585/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                  <  1:150811/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                  <  1:2239857/150‑1 (MQ=255)

                         acgcaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGaa                                                                                                                               <  2:377650/147‑1 (MQ=255)

                         acacAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGaa                                                                                                                               >  1:1536975/1‑150 (MQ=255)

                            caGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGAGACGCTGACATGCCGGTCACCAGCCCTGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAg                                                                                                                            <  2:1147303/150‑1 (MQ=255)

                                ccGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCg                                                                                                                        >  1:679145/1‑150 (MQ=255)

                                 cGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGa                                                                                                                       >  1:1726747/1‑150 (MQ=255)

                                  gACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtac                                                                                                                                              >  1:592172/1‑123 (MQ=255)

                                   aCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                                     <  1:800622/150‑1 (MQ=255)

                                     cAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGca                                                                                                                   >  2:847935/1‑150 (MQ=255)

                                      aGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGcac                                                                                                                  <  2:2316209/150‑1 (MQ=255)

                                         gcCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAgc                                                                                                               <  1:1680152/150‑1 (MQ=255)

                                             gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggt                                                                                                                                                >  2:442774/1‑112 (MQ=255)

                                             gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggt                                                                                                                                                <  1:442774/113‑2 (MQ=255)

                                                gATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGt                                                                                                       <  2:1649542/150‑1 (MQ=255)

                                                 aTTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtac                                                                                                                                              <  2:2321324/111‑4 (MQ=255)

                                                 aTTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTg                                                                                                      <  1:579801/150‑1 (MQ=255)

                                                   tACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGAt                                                                                                    <  2:1704552/150‑1 (MQ=255)

                                                   tACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGAt                                                                                                    <  2:281124/150‑1 (MQ=255)

                                                   tACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGAt                                                                                                    <  2:1609584/150‑1 (MQ=255)

                                                     ccGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCt                                                                                                  >  1:2105374/1‑150 (MQ=255)

                                                        tGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGaa                                                                                               <  2:710923/150‑1 (MQ=255)

                                                            tATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAg                                                                                           >  1:779010/1‑150 (MQ=255)

                                                              tCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGcc                                                                                         >  2:474505/1‑150 (MQ=255)

                                                                 gtgtGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGc                                                                                      >  2:1977823/1‑150 (MQ=255)

                                                                 gtgtGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGc                                                                                      >  2:806944/1‑150 (MQ=255)

                                                                   gtgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCca                                                                                    <  2:1755451/150‑1 (MQ=255)

                                                                   gtgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCca                                                                                    <  2:2105374/150‑1 (MQ=255)

                                                                   gtgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCca                                                                                    <  1:255060/150‑1 (MQ=255)

                                                                    tgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCcag                                                                                   >  2:1477751/1‑150 (MQ=255)

                                                                      gCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCcagca                                                                                 <  1:812025/150‑1 (MQ=255)

                                                                      gCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCCGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGACGGCcagca                                                                                 >  2:1002431/1‑150 (MQ=255)

                                                                       cAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACA‑CGCGCTGGTGATCTGAACTAGCCGGCcagcag                                                                                <  1:62675/150‑1 (MQ=255)

                                                                            tGCGACGCAGACATGCCGGCCCCCAGCCTCTACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTccc                                                                           >  1:1147083/1‑150 (MQ=255)

                                                                            tGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTccc                                                                           <  2:1726747/150‑1 (MQ=255)

                                                                            tGCGACGCAGACATGCAGGCACATCGCATTGACGAAGTGCTGGACATTTCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTccc                                                                           >  2:25642/1‑150 (MQ=255)

                                                                              cGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCTTATGCTCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCAc                                                                         <  2:790386/150‑1 (MQ=255)

                                                                                aCGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                       >  1:2269542/1‑150 (MQ=255)

                                                                                  gCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACa                                                                     >  2:1457107/1‑150 (MQ=255)

                                                                                       cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                                <  2:2176589/150‑1 (MQ=255)

                                                                                        aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGa                                                                                                     <  1:1685713/112‑1 (MQ=255)

                                                                                            cGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGAt                                                           <  1:1558231/150‑1 (MQ=255)

                                                                                             ggCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATc                                                          >  1:570093/1‑150 (MQ=255)

                                                                                              gCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCa                                                         <  1:1477751/150‑1 (MQ=255)

                                                                                                 cccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGca                                                      >  1:226451/1‑150 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                             <  1:156715/79‑1 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                             >  2:156715/1‑79 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcac                                                     >  2:387213/1‑150 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcac                                                     >  2:640602/1‑150 (MQ=255)

                                                                                                     gCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcacaca                                                  <  1:847935/150‑1 (MQ=255)

                                                                                                        tCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGACAGCCGGTCCCACGACCGCTTGGCGATCAGCACACacaa                                               <  1:2322848/150‑1 (MQ=255)

                                                                                                           acgaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAATCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGc                                                       <  1:2111828/139‑1 (MQ=255)

                                                                                                             gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                                     <  2:2446700/76‑1 (MQ=255)

                                                                                                             gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                                     >  1:2446700/1‑76 (MQ=255)

                                                                                                                ccgctgGCCAGTGCCGAGCAGGCGACGTCCCGCGAAAACACTTGAGATGGTAACTGCTTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCa                                       <  2:2231421/148‑1 (MQ=255)

                                                                                                                 tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                      >  1:128773/1‑150 (MQ=255)

                                                                                                                 tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                      >  1:377481/1‑150 (MQ=255)

                                                                                                                    tgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGaa                                   <  1:1457107/150‑1 (MQ=255)

                                                                                                                     ggCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGAACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAAc                                  >  1:1280784/1‑150 (MQ=255)

                                                                                                                      gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACa                                 <  1:1231526/150‑1 (MQ=255)

                                                                                                                      gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACa                                 >  1:1187068/1‑150 (MQ=255)

                                                                                                                       ccAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAg                                >  1:525285/1‑150 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGc                                                                                                                          >  1:104247/1‑60 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGc                                                                                                                          <  2:104247/60‑1 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGc                               >  1:302574/1‑150 (MQ=255)

                                                                                                                           tGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcg                            >  1:794893/1‑150 (MQ=255)

                                                                                                                            gCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                           <  2:779010/150‑1 (MQ=255)

                                                                                                                            gCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                           <  1:25642/150‑1 (MQ=255)

                                                                                                                               gACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGa                        <  2:1351099/150‑1 (MQ=255)

                                                                                                                                aCCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGTGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGAc                       >  1:988145/1‑150 (MQ=255)

                                                                                                                                      cgcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                 <  2:570093/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                 <  2:2390257/149‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                 >  1:2390257/1‑149 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg                <  2:1536975/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg                <  2:525285/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg                <  2:226451/150‑1 (MQ=255)

                                                                                                                                           tACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAc            <  1:1977823/150‑1 (MQ=255)

                                                                                                                                            aCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgcg                          >  2:242164/1‑135 (MQ=255)

                                                                                                                                            aCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAcc           >  2:474567/1‑150 (MQ=255)

                                                                                                                                             ccGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAccc          <  2:2000128/150‑1 (MQ=255)

                                                                                                                                                  ccAAGGCCCGGGGCCGCAACCGCGGCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCg     >  2:1623934/1‑150 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                                >  2:2350891/1‑90 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                                <  1:2350891/90‑1 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGt    >  2:409981/1‑150 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGt    >  1:2134737/1‑150 (MQ=255)

                                                                                                                                                    aaGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTa   >  2:2262495/1‑150 (MQ=255)

                                                                                                                                                     aGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTAc  <  2:338368/150‑1 (MQ=255)

CCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTAC  >  NC_002947/4438105‑4438403

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